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RESUMEN 


Palabras clave: grafeno, aislantes topológicos, teoría de Floquet, transporte electró- 
nico, estados de borde. 


Números PACS: 73.63.-b; 73.22.Pr; 73.20.At; 72.80.Vp; 78.67.-n 


El grafeno es el material más delgado que existe y reúne muchas propiedades 
récord, incluyendo conductividad térmica y eléctrica, y resistencia mecánica. Desde 
su aislamiento y caracterización en 2004, ha comenzado una carrera de científicos 
en diversas áreas de la física, química, y la ciencia de materiales para entender y 
explotar las propiedades de este notable material. 


Esta tesis está dedicada al estudio teórico del grafeno iluminado por láseres 
intensos. Usualmente los efectos de la luz sobre la materia son perturbativos. De 
hecho la luz se utiliza como una herramienta para observar diferentes aspectos de 
un material, como por ejemplo información estructural vía espectroscopia Raman. 
En esta tesis, sin embargo, nos enfocamos en efectos no perturbativos de la luz 
sobre el material. Mostramos cómo un láser puede modificar la estructura electró- 
nica del grafeno, alterando sus propiedades de conducción, y generando también 
características topológicas ausentes en el material no iluminado. 


En los dos primeros capítulos de esta tesis discutiremos las modificaciones 
en la estructura electrónica tanto de un plano de grafeno como de distintos tipos 
de nanocintas de grafeno. En ambos casos, vemos que la radiación abre brechas 
energéticas en las bandas de grafeno. Este efecto podría ser usado para aprovechar 
las notables propiedades de conducción del grafeno, a la vez que un láser sería 
el mecanismo de encendido y apagado de la corriente, abriendo posibilidades al 
diseño de nanodispositivos optoelectrónicos basados en grafeno. 


En el tercer capítulo analizamos en detalle las brechas electrónicas inducidas por 
láser en una nanocinta de grafeno y veremos como, debido a la radiación, se forman 
estados localizados en los bordes de la nanocinta. Estos estados se propagan en 
direcciones opuestas en cada borde (son quirales), y debido a esta separación 
espacial de los estados disponibles, no hay retrodispersión y el transporte es robusto 
frente a defectos del material. Estos resultados junto con un cálculo del invariante 
topológico para el grafeno irradiado nos permite concluir que estos estados son 
topológicos, de manera que el grafeno irradiado constituye un aislante topológico de 
Floquet. 


vil 
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En el cuarto capítulo estudiamos la conductancia en una configuración de dos y 
seis terminales. Esta última configuración nos permite calcular la resistencia Hall, 
que podemos comparar con los estados disponibles en la relación de dispersión del 
grafeno irradiado. Vemos que la cuantización de la conductancia en los aislantes 
topológicos de Floquet no sigue las mismas reglas que en el efecto Hall cuántico 
entero, donde la conductancia está cuantizada en el número de estados quirales, 
que a su vez es equivalente al invariante topológico que caracteriza las bandas 
electrónicas del material. En cambio, la conductancia de los aislantes topológicos de 
Floquet está aproximadamente cuantizada en el número de canales activos para el 
transporte coherente, que son el número de estados que llevan un peso espectral 
significativo en la densidad de estados promediada en el tiempo. Los resultados 
obtenidos tienen como ejemplo un sistema basado en grafeno pero son fácilmente 
generalizables a cualquier sistema que albergue electrones de tipo Dirac. 

El capítulo quinto contiene los últimos resultados de esta tesis, donde se ca- 
racterizan en profundidad las propiedades topológicas de un aislante topológico de 
Floquet a bajas frecuencias. El espectro de Floquet de un sistema de tipo Dirac 
descrito por el formalismo de Floquet, en particular el grafeno irradiado, se vuelve 
más complejo y difícil de analizar, debido a la multiplicidad de los procesos fotónicos 
de absorción/emisión disponibles. En este capítulo daremos un método para des- 
cribir constructivamente la topología y el espectro de Floquet de bajas frecuencias, 
a medida que se consideran más procesos fotónicos. Develamos que existe una 
jerarquía de brechas energéticas anidadas una dentro de otra, y ordenadas por el 
número de procesos fotónicos involucrados en la apertura de la brecha. Además, 
el número de orden del proceso fotónico también se relaciona directamente con 
el número de estados topológicos de Floquet que ocurren dentro de cada brecha 
energética. Este método sirve para el régimen de bajas frecuencias, fuera de este 
rango debemos recurrir al cálculo explícito del invariante topológico que describe el 
número total de estados topológicos. Haciendo uso de este formalismo damos un 
mapa topológico, para un amplio rango de frecuencias e intensidades de radiación, 
proveyendo así un estudio detallado de la topología de los aislantes topológicos de 
Floquet en sistemas de tipo Dirac. 

En el último capítulo se presentan las conclusiones, además de una discusión 
sobre el contexto en el que se desarrolla el trabajo de tesis. 
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Abstract 


En este primer capítulo introduciremos algunas de las ideas básicas para entender 
el contexto en el que se desarrolla esta tesis. El aislamiento del grafeno, junto con la 
caracterización de sus notables propiedades han abierto nuevas ramas de la física teórica 
tanto como experimental. Una de ellas ha cobrado gran importancia recientemente: los 
aislantes topológicos. Descubriremos cómo la estructura del grafeno ha sido utilizada como 
base de un primer modelo simplificado de aislantes topológicos. Mostraremos cómo el 
grafeno y los aislantes topológicos están fuertemente relacionados ya que la radiación de un 
láser nos permite hacer del grafeno un aislante topológico de Floquet donde las propiedades 
topológicas son inducidas y controladas por los parámetros del láser. 


Esta tesis se desarrolla cerca de diez años después del descubrimiento del 
grafeno y siete años después de la realización experimental de un aislante topológico 
sin campo magnético. Desde entonces, la comunidad dedicada a la investigación 
del grafeno, sus propiedades y sus aplicaciones ha crecido enormemente dentro de 
la física. 

Paralelamente, se han incrementado los esfuerzos para desarrollar nuevos 
aislantes topológicos en dos y tres dimensiones y se ha desarrollado la teoría del 
orden topológico, necesaria para entender estos nuevos estados electrónicos. Hoy 
en día, se utiliza esta teoría para explicar el orden topológico de superconductores, 
excitaciones acústicas, plasmónicas y fotónicas, entre otras. 


1. Introducción 


1.1 Estructura electrónica del grateno 





Andre Geim también recibió el Pre- 
mio Ig Nobel a los trabajos curio- 
sos por hacer levitar una rana en el 
año 2000. Imagen tomada de «Of 
flying frogs and levitrons» Berry 
y Geim (1997) [3]. 





Fig. 1.1: Red hexagonal para el 
grafeno, con detalle de los orbita- 
les atómicos hibridizados sp? en el 
plano xy, que forman los enlaces 
covalentes o, mientras que los or- 
bitales atómicos p. forman los or- 
bitales moleculares , que son los 
que le dan las propiedades eléc- 
tricas al grafeno a bajas energías. 
Figura adaptada de [4]. 


La interacción de la luz con la materia une dos grandes áreas de la física: 
los sistemas de baja dimensionalidad como el grafeno y los aislantes topológicos, 
abriendo el camino para explorar las distintas propiedades que pueden ser inducidas 
externamente en un material. En esta tesis nos basamos en el estudio teórico de 
los efectos de la radiación láser sobre grafeno y explicamos las modificaciones que 
se producen en la estructura electrónica y el cambio no trivial en sus propiedades 
topológicas. El desarrollo de este trabajo se estructura en los siguientes objetivos: 


= Entender la estructura electrónica del grafeno y las modificaciones que produce 
la radiación, utilizando el formalismo de Floquet. 


= Obtener resultados físicos y predicciones observables, simulando mediciones de 
densidad de estados y obteniendo cálculos de transporte en una configuración 
de dos terminales. 


= Contrastar la cuantización de la conductancia de un aislante topológico de Floquet 
respecto de un aislante topológico estático. Explicar esta cuantización mediante 
las características del espectro y de las funciones de onda involucradas. 


= Caracterizar la topología de las bandas en el espectro de cuasienergías que se 
obtiene al someter al grafeno a un potencial periódico en el tiempo utilizando la 
teoría de Floquet. 


1.1. Estructura electrónica del grafeno 


El grafeno es uno de los materiales más simples que conocemos, sin embargo, 
su aislamiento y caracterización ha sido elusiva por mucho tiempo. Recién en el año 
2004, el grafeno fue aislado y sus notables propiedades fueron caracterizadas por 
Andre Geim y Konstantin Novoselov [1, 2]. Sus avances les merecieron el Premio 
Nobel de Física en 2010 ["The Nobel Prize in Physics 2010". Nobelprize.org. Nobel 
Media AB 2014. http://www.nobelprize.org/nobel_prizes/physics/laureates/2010 ]. 

Podría pensarse al grafeno como el material constituyente del grafito. Basta 
apilar muchas capas atómicas de grafeno hasta obtener un material bulk . De 
hecho, el proceso inverso fue usado para aislar en 2004 al grafeno. El método de la 
cinta Scotch consiste en «pelar» una muestra de grafito utilizando cinta adhesiva 
y decapar hasta obtener una capa de un átomo de espesor. Sin embargo, las 
propiedades del grafito y del grafeno son muy distintas, el grafito es blando (incluso 
usado como lubricante) y moderadamente conductor, mientras que el grafeno junto 
a los nanotubos de carbono constituyen los materiales más duros conocidos y de 
mayor conductividad eléctrica. 

Este nuevo material de tan sólo un átomo de espesor hecho puramente de 
carbono tiene propiedades que desafían el sentido común (clásico). Es el material 
más duro conocido: su módulo de Young es incluso mayor que para el diamante; y 
su conductividad térmica y eléctrica son mayores que las del cobre. 


1. Introducción 


Comenzamos abordando el primer objetivo de la tesis, describiendo la estructura 
electrónica del grafeno. Muchas de las peculiares propiedades de transporte en 
grafeno, tienen su origen en su estructura de bandas, que se conoce desde el 
año 1947, mucho antes que el grafeno mismo ya que fue analizada por Wallace 
como punto de partida para describir la estructura de bandas del grafito [5]. Esta 
descripción predice que el grafeno es un semi-metal o un semiconductor de gap cero, 
lo cual se puede deducir con un modelo muy simple. Como vemos en la Fig. 1.1, se 
trata de una red hexagonal, formada por átomos de carbono, cuyos orbitales están 
hibridizados en el tipo sp? en el plano xy, para formar orbitales o sus tres vecinos, 
mientras que los orbitales atómicos p. que forman enlaces 7 y dan las propiedades 
eléctricas y de transporte a bajas energías. 





El modelo más simple considera orbitales p, en una red hexagonal, donde sólo 
interaccionan con sus primeros tres vecinos [6], como se muestra en la Fig. 1.2 (a). 
La red hexagonal está formada por dos redes triangulares A y B, intercaladas, 
denotadas con círculos llenos y vacíos. Los vectores de red 


aj = a(v3/2,1/2) 
az = a(v3/2,-1/2), 


expanden la red triangular 4, siendo a = V3 ay la constante de red para grafeno y 
ay = 0,142 nm la distancia entre átomos de carbono. Cada sitio de la red B está 
desplazado en d;, ¿ =1,2,3, con 


d; = ag(1,0) 
d> = ap(-1/2,-V3/2) 
d3 = a0(1/2,V3/2), 


donde podemos verificar que d;, da = (d¡ — a¡) y d3 = (d¡ - az) conforman los 
vectores de desplazamiento hacia los tres primeros vecinos en la Fig. 1.2 (a). Los 
vectores d; conectan sitios de tipo A con sitios de tipo B, por lo que un sitio de tipo 
A tiene primeros vecinos sólo de tipo B. 


1.1 Estructura electrónica del grafeno 


Fig. 1.2: (a) Red hexagonal de 
grafeno con sitios de tipo A y B 
denotados con círculos llenos y 
vacíos respectivamente. Los vec- 
tores (a, az son los vectores de 
red. (b) Red recíproca, expandida 
por los vectores (b:,b2). La prime- 
ra zona de Brillouin (BZ) se mues- 
tra en gris. Los puntos K, y K_ se 
encuentran en vértices no equiva- 
lentes de la BZ. El punto M sobre 
la arista del hexágono, se encuen- 
tra entre medio de los puntos XK. 
y el punto I' en el centro de la BZ. 
Figura adaptada de [4]. 
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1.1 Estructura electrónica del grateno 


En la Fig. 1.2 (b), observamos los vectores de la red recíproca b, y b>, que 
cumplen la relación a,b; = 270, y están dados por 


b1 = b(1/2, 3/2) 
b) = b(1/2,-V3/2) > 


donde b= 47/(v3a). La primera zona de Brillouin (BZ), construida como la celda 
de Wigner-Seitz de la red recíproca, toma forma hexagonal con seis vértices, de los 
cuales sólo dos son no equivalentes, mientras que los otros están desplazados en 
vectores de red recíproca. Estos puntos los denotamos con K',, tales que 


K, = 411/3a(W3/2,-1/2) K-_=41/3a(V3/2,1/2), 


y son de gran importancia para describir la relación de dispersión a bajas energías. 
Para obtener la relación de dispersión describimos esta red con el modelo 
tight-binding utilizando el Hamiltoniano en segunda cuantización 


H = Y ejtjt! = y ea + h.c., 
1 (5,5) 

donde h.c. denota el Hermitiano conjugado, e, es la energía del sitio ¿, la suma sobre 
(i, j) recorre los sitios ¿ y sus primeros vecinos j y t; ¿ es el coste energético en que 
un electrón salte del sitio ¿ al :, llamada amplitud de hopping. 

En el caso de una red con invarianza traslacional, podemos aplicar el teorema 
de Bloch y obtenemos el Hamiltoniano de Bloch, cuyas autofunciones describen a 
las funciones de onda dentro de una celda unidad. En el caso en que los hoppings 
en cada dirección d;, con 1 = 1,2, 3, sean homogéneos en toda la red y las energías 
de sitio de cada subred sean constantes, el Hamiltoniano de Bloch para grafeno 
toma la forma 


Hop (k) > 0) | (1) 


S(k) ep 


donde S(k) = Yj-1 teilik es la suma de las fases que adquiere la onda de Bloch al 
desplazarse a sus vecinos y S(k)* su complejo conjugado. Fácilmente podemos 
diagonalizar este sistema y encontrar la relación de dispersión para el grafeno. Si 
suponemos que las energías de sitios A y B son idénticas, podemos tomarlas como 
la energía de referencia e igualarlas a cero. El espectro de energías del Hamiltoniano 
en la Ec. (1.1) con e 4 = e = O se muestra en la Fig. 1.3 y está descrito por la relación 
de dispersión 


e(k). = +|5(k)| 


k k 
= +t1/1+4c0s a cos =— +4cos? HE, 
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El espectro es simétrico respecto de e = O y mientras no existan hoppings a 
segundos vecinos se preserva la simetría electrón-hueco: esto quiere decir que 
si existe un estado 4, (k) con energía e(k), existe también un estado 4_(k) con 
energía —-e(k). En grafeno existen hoppings a segundos vecinos que rompen la 
simetría electrón-hueco, aunque en grafeno se preserva la simetría efectivamente 
para bajas energías. 


La banda de conducción y de valencia se unen en los dos puntos XK, no equi- 
valentes de la zona de Brillouin. Si k es el vector de onda medido desde alguno 
de estos puntos, podemos expandir la relación de dispersión de la Ec. (1.1) para 
pequeños valores de k y para el valle K, obtenemos 


es roy a (12) 


donde up = v3/2ta/h es la velocidad de Fermi para grafeno, k, = k, +ik,. La 
expansión alrededor del punto K_ cumple Hx_ = Hx... El Hamiltoniano de la 
Ec. (1.2) puede ser expresado como 


Hx,(k) = hop(07ky + 0yky) 41.3) 
= hvpo . k, 


AE 


son las matrices de Pauli y actúan sobre el grado de libertad de subred. En analogía 


donde 


con los sistemas donde o es un operador que actúa sobre las variables de espín, 


Fig. 1.3: Relación de dispersión 
para grafeno. La banda de valen- 
cia (azul) y la banda de conduc- 
ción (rojo) se unen en dos puntos 
no equivalentes a energía cero, for- 
mando los conos de Dirac. Para 
bajas energías la relación de dis- 
persión es lineal. 
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Fig. 1.4: Conos de Dirac a bajas 
energías. Sólo dos de los conos 
que aparecen en los vértices de la 
primera zona de Brillouin, son no 
equivalentes, llamados valles K,. 
El pseudoespín en cada valle tiene 
dirección paralela o antiparalela al 
vector de onda k y a la vez, se in- 
vierte según el estado electrónico 
pertenezca a la banda de valencia 
o ala banda de conducción. Figura 
adaptada de [4]. 


donde (0) nos da la dirección del espín, llamamos al grado de libertad de subred el 
pseudoespín, y el valor de espetación (o) nos da la dirección del pseudoespín 
Podemos verificar que el Hamiltoniano es proporcional al operador de helicidad 


h=0-k/|k|. (1.5) 
El operador h es la proyección del pseudoespín sobre la dirección del momento k de 
la partícula. Inmediatamente notamos que el Hamiltoniano del sistema conmuta con 
el operador de helicidad, por lo que la proyección del pseudoespín sobre la dirección 
del vector k es una cantidad bien definida y conservada. Dicho de otro modo, un 
estado con momento k, que sea autoestado del Hamiltoniano k - p, debe tener el 
pseudoespín bien definido, y en la dirección de k, ya sea paralelo o antiparalelo. 

El Hamiltoniano de la Ec. (1.2) es el Hamiltoniano k-p para grafeno, cuyo 
espectro 


e(k) = +up |k] 


es lineal en k, lo que da a la dispersión su forma cónica, característica de los 
fermiones de Dirac. 


E antiparalelo E 


paralelo 





paralelo antiparalelo 


momento P_——Y 


(medido desde K+.) 
pseudospin —> 


Se puede ver en la Fig. 1.4 esquemáticamente los conos de Dirac, ubicados en 
las seis esquinas de la zona de Brillouin. Se destacan los dos conos no equivalentes 
en los valles K,. Los pseudoespines son paralelos (antiparalelos) a la dirección 
de propagación k para un electrón en la banda de conducción del valle K, (K_), 
mientras que en la banda de valencia la relación se invierte en ambos valles. Esta 
propiedad se la conoce como la quiralidad de los fermiones de Dirac, y tiene 
importantes consecuencias en las propiedades electrónicas y de transporte. Algunos 
de los efectos de los fermiones quirales son: (i) la ausencia de retrodispersión, (ii) el 
tunelamiento de Klein, (iii) particularidades en la antilocalización, entre otros. 


1.1.1. Fermiones de Dirac 


Una de las consecuencias más importantes de la relación de dispersión de bajas 
energías en grafeno, es la linealidad y quiralidad de las excitaciones electrónicas. A 
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este tipo de fermiones se los conoce como fermiones de Dirac. El término fermiones 
de Dirac ha sido usado en la descripción relativista de electrones libres y es tomado 
por la física de materia condensada para describir excitaciones electrónicas dentro 
de la red cristalina de semiconductores de gap pequeño (gap cero en grafeno). Las 
excitaciones de bajas energías en grafeno pueden ser descritas por una relación de 
dispersión relativista, expandida en dos puntos independientes de la zona de Brillouin 
y donde el espín es remplazado por el espinor que describe la fase relativa entre las 
subredes A y B. La interacción espín-órbita (SO) que puede modificar la estructura 
de bandas cerca de la energía de Fermi desdoblando los niveles energéticos para 
los distintos espines. En grafeno la interacción SO es despreciable en grafeno y los 
distintos espines tienen prácticamente la misma relación de dispersión, a menos que 
se incremente externamente el acople SO mediante un campo magnético o bien, 
mediante acoplamiento SO inducido por proximidad. 

Las excitaciones electrónicas al interior del grafeno cumplen con las ecuaciones 
de Dirac y que se comportan como fermiones sin masa. Se propagan como los 
fotones, con una relación de dispersión lineal, es decir a velocidad constante. Esta 
velocidad está dada por los parámetros del material y se conoce como velocidad de 
Fermi uy = 10%m/s, tan solo unas 300 veces más pequeña que la velocidad de la 
luz. 


1.1.2. Tunelamiento de Klein 


Uno de los efectos más interesantes en grafeno, también es consecuencia de 
los efectos relativistas que cumplen los fermiones sin masa: La paradoja de Klein. 
Este efecto hace que las excitaciones electrónicas en grafeno y los electrones 
quirales en general, tengan probabilidad 1 de atravesar una barrera electrostática, 
a diferencia de otros materiales. La propagación de los electrones quirales se ve 
poco afectada por la presencia de perturbaciones externas y desorden que pueden 
tratarse como potenciales electrostáticos localizados. El efecto de tunelamiento de 
Klein da transmisión perfecta para un electrón que incide en dirección normal de 
una barrera. Si la incidencia no es normal, la trasmisión decrece con una función 
que depende de la altura de la barrera. Aún más paradójico resulta que una barrera 
suave es más efectiva en reflejar un electrón (cuando se propaga en una dirección 
arbitraria), mientras que una barrera más abrupta es más transparente. Incluso, si 
la barrera separa regiones de excitaciones de electrones y huecos, puede haber 
transmisión de electrón a hueco con una eficacia del 100 % [7, 8, 9]. 

Cuando el electrón se desplaza en la dirección del vector de onda k el pseu- 
doespín se conserva en ausencia de scattering entre distintos valles, por lo que la 
transmisión en incidencia normal a una barrera cuadrada, es perfecta. Como vimos 
en la Fig. 1.4, a una misma energía el pseudoespín tiene direcciones opuestas en 
el valle K, y K_, de modo que si se conserva el pseudoespín y se desprecia el 
scattering entre distintos valles, el único estado accesible al otro lado de la barrera 
es el estado que se propaga en la misma dirección k. De esta forma se obtienen 
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1.2 Aislantes topológicos 


caminos libres medios muy largos del orden del micrón, que presenta un desafío en 
cuanto a confinar los electrones en grafeno. 


1.1.3. Interés para la física y aplicaciones tecnológicas 


La familia de materiales bidimensionales, fundada con el descubrimiento del 
grafeno, continúa creciendo y siendo objeto de intensa investigación. Algunos de 
sus miembros son disulfuro de molibdeno (dicalcogenuros en general), fosforeno, 
nitruro de boro hexagonal, siliceno, entre otros. 

La investigación y el financiamiento dedicados a estos nuevos materiales ya 
encuentran muchas aplicaciones, que varían desde interruptores ópticos hasta fo- 
tocorrientes bombeadas para la recolección de energía [10]. Las aplicaciones en 
optoelectrónica también son vastas: ejemplos de dispositivos plasmónicos se en- 
cuentran en las Refs. [11, 12, 13, 14] y de fotodetectores en [15, 16]. En particular 
en grafeno, los campos en las frecuencias de los terahertz producen efectos no 
lineales [17], pueden producir fotoconductividad [18] y guiar el diseño de fototransis- 
tores [19]. 


1.2. Aislantes topológicos 


Los fermiones de Dirac son excitaciones electrónicas en las que el grado de 
espín está directamente relacionado con la dirección de propagación o con el vector 
de onda. Este tipo de excitaciones caracteriza los estados topológicos que aparecen 
en la superficie de los aislantes topológicos con simetría de reversión temporal. En 
esta sección estudiamos cómo los aislantes topológicos, los fermiones de Dirac y 
el grafeno se encuentran ligados desde su desarrollo teórico hasta su realización 
experimental. 

Los aislantes topológicos representan una nueva fase electrónica, donde el inte- 
rior del material es aislante, mientras que existen estados metálicos en la superficie 
o en los bordes. Estos estados electrónicos de borde están protegidos por las sime- 
trías que protegen al invariante topológico. A diferencia de los aislantes comunes 
(triviales), semiconductores y metales, estos materiales tienen propiedades que caen 
fuera de la descripción usual. La teoría de bandas fue un gran logro en la física de la 
materia condensada para describir la estructura electrónica de muchos materiales, 
en particular de los semiconductores. Sin embargo, desde los 80's con el descubri- 
miento del efecto Hall cuántico (QHE) [20] en gases electrónicos bidimensionales 
(2DEG) inmersos en un fuerte campo magnético, el orden topológico [21] de los 
sistemas debió incorporarse a la teoría de bandas para describir las propiedades de 
transporte de estos sistemas, como la cuantización de la conductancia Hall [22]. 


1.2.1. El nacimiento de los aislantes topológicos 


El modelo de Haldane, propuesto en 1988, es el primer modelo teórico que 
reproduce el efecto Hall cuántico (QHE). En un principio su importancia se vio 
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disminuida por tratarse de un artificio sin una relación clara con un sistema físico, 


donde los hoppings a segundos vecinos en una red hexagonal son imaginarios, Fig. 


1.5, con signos opuestos para cada subred. 

El modelo de Haldane incorpora hoppings a segundos vecinos, que corresponden 
a hoppings entre sitios de la misma subred. Estos hoppings imaginarios cambian 
de signo según vinculen sitios pertenecientes a la subred 4 o a la subred B. El 
Hamiltoniano del modelo de Haldane es 


H=tYN ae + Y mjotl+he, 
(1,5) ((6,5)) 

donde ((i, j)) denota la sumatoria a segundos vecinos y »; ¿ = —1;¿ = +i depende 
de la orientación del segundo vecino: si para ir del sitio ¿ al ¿ debemos hacer un giro 
a derecha, el signo es positivo y si el giro es a izquierda, el signo es negativo. En 
la Fig. 1.6 podemos ver la relación de dispersión del modelo de Haldane para una 
cinta de tipo zigzag, en la que se observan estados de borde que cruzan el gap e 
energía e = 0. 

Este modelo fue extendido por Kane y Mele para predecir el efecto espín-Hall 
cuántico (QSHE) en grafeno con acople espín-órbita (SO) lo suficientemente grande. 
Según este modelo, por medio acople SO intrínseco los electrones tienen acoples 
imaginarios a segundos vecinos, tal como en el modelo de Haldane, pero con signos 
opuestos para cada espín. En este caso existen dos estados de borde, que se 
propagan en direcciones opuestas y tienen espín opuesto. Esto es una consecuencia 
de que, a diferencia del QHE, el QSHE conserva la simetría de reversión temporal y 
está caracterizado por un nuevo invariante topológico del tipo Z2. A pesar de que en 
grafeno el acople SO es muy débil, poco después de la propuesta de Kane y Mele, 
el efecto QSH fue predicho y comprobado experimentalmente en pozos cuánticos 
de HgTe/CaTe. 


1.3. Sistemas forzados periódicamente 


Los sistemas cuánticos sujetos a fuerzas periódicas pueden alterar su estructura 
electrónica, lo que brinda gran control sobre las propiedades de transporte [24, 
25]. Los mecanismos para forzar estos sistemas pueden ser, por ejemplo, un gate 
periódico, para inducir bombeo cuántico [26, 27], o bien un campo electromagnético. 
En el último caso, muchas propiedades interesantes emergen, tales como efectos 
termoeléctricos [10], respuestas no lineales [28] y transiciones de fase cuánticas [29, 
30]. 

La luz, más específicamente, el campo electromagnético al que está sujeta una 
muestra de grafeno u otro material de Dirac, cuando es irradiada por un láser pueden 
producir la apertura de gaps en el espectro [30, 31, 32], donde se desarrollan fases 
del tipo QHE o QSHE. Una gran variedad de materiales bidimensionales bajo la 
incidencia de luz, exhibe, por ejemplo, estados de tipo Hall polarizados en espín o 
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Fig. 1.5: Red hexagonal con hop- 
pings imaginarios a segundos ve- 
cinos. La fase ganada en el hop- 
ping tiene el mismo signo cuando 
la dirección de hopping es la in- 
dicada por las flechas. El sentido 
de los hoppings entre sitios A es 
antihorario (lineas de puntos azu- 
les) y entre los sitios B es horario 
(líneas de trazos verdes). Figura 
adaptada de [4]. 





Fig. 1.6: Relación de dispersión 
para el modelo de Haldane, con 
los mismos parámetros que en Ka- 
ne y Mele (2005) [23], con t'/t = 
0,03, siendo t” el módulo del hop- 
ping imaginario a segundos veci- 
nos. Para esta geometría (zigzag) 
los estados de borde conectan los 
valles K.,. 
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1.4 Sobre esta tesis 


en valle en disulfuro de molibdeno (MoS>) [33], corrientes polarizadas por valle en 
grafeno bicapa [31] y estados sintonizables en grafeno [34]. 

De manera análoga a los aislantes topológicos, que pueden ser usados para 
obtener canales coherentes de transporte de electrones sin disipación, los sistemas 
cuánticos sujetos a perturbaciones periódicas, desarrollan una estructura de bandas 
que se ve modificada y puede tener propiedades topológicas [35, 36]. La teoría de 
Floquet predice que las propiedades topológicas pueden ser inducidas incluso en 
materiales triviales [37, 38]. La luz podría ser usada como el mecanismo de control 
sobre los estados topológicos ya sea activando y desactivando estos canales, o bien 
cambiando su quiralidad y su dirección de propagación. 

La actividad en este campo viene incrementando y se vio acelerada por la ob- 
servación de estados de Floquet-Bloch en la superficie de un aislante topológico 
tridimensional [39]. Los estados topológicos también son simulados en redes ópti- 
cas [40], en átomos fríos [41], donde el campo de gauge es producido artificialmente. 


1.4. Sobre esta tesis 


El objetivo de esta tesis es explorar las propiedades de las nuevas fases electró- 
nicas que la luz induce en grafeno y de manera similar en otros sistemas de Dirac. 
Por un lado, el grafeno brinda un candidato ideal para aplicaciones tecnológicas y 
nanodispositivos por sus sorprendentes propiedades de conducción, gran movili- 
dad electrónica y por supuesto, la posibilidad de tener dispositivos nanométricos 
eficientes. Por el otro, la nueva fase electrónica de aislantes topológicos también 
brinda propiedades con gran potencial aplicación, como canales de conductancia 
cuasi-unidimensionales, ausencia de retrodispersión en el caso de aislantes tipo 
QHE y polarización en espín para el caso QSHE. Más aun, inducir propiedades 
topológicas en grafeno mediante la aplicación de un campo electromagnético une 
todas estas características en un sólo sistema y brinda gran capacidad de control 
sobre los estados topológicos. 
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Abstract 


El grafeno es uno de los principales candidatos a ser la base de nuevos nanodispositivos, 
debido a sus excelentes propiedades de conducción y gran movilidad electrónica. Sin 
embargo, carece de un gap que provea un mecanismo de encendido y apagado. Un láser 
provee la radiación coherente necesaria para inducir gaps en sistemas bidimensionales, en 
particular en grafeno. Al tratar a la radiación como una perturbación periódica en el tiempo 
podemos utilizar la teoría de Floquet para obtener el espectro electrónico donde se inducen 
gaps en el promedio temporal de la densidad local de estados a energías conmensuradas 
con la mitad de la energía fotónica. 


Un láser puede producir modificaciones en la estructura electrónica del grafeno. 
Con el objetivo de abrir un gap en la densidad de estados promedio, investigamos 
los efectos de la radiación en grafeno, tanto para un sistema bulk como para distintos 
tipos de nanocintas, explorando los efectos que tiene la polarización de la radiación 
en cada caso. 

Para describir al grafeno irradiado utilizamos la teoría de Floquet, en la que 
se incluyen los efectos de un campo electromagnético coherente como una per- 
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turbación temporal. En la primera sección se presenta la teoría de Floquet, como 
herramienta para analizar el espectro de cuasienergías de un Hamiltoniano sometido 
a un potencial dependiente periódicamente del tiempo. Seguido de esto, podremos 
entender las primeras evidencias teóricas donde se observan gaps inducidos por 
radiación en grafeno lejos del punto de Dirac [32]. Estos resultados constituyen el 
estado del arte de la teoría que describe estos sistemas al comienzo de esta tesis y 
forman la base sobre la cual construimos nuevos resultados. En la tercera sección 
se muestran los primeros resultados propios, investigando cómo se desarrollan 
los gaps en nanocintas de grafeno irradiado, según la dirección de la polarización 
respecto del borde de la cinta y del tipo de borde [42, 43]. Los gaps en el espectro 
de cuasienergías se reflejan en gaps de transporte a dos terminales. El modelo 
de scattering utilizado para hacer los cálculos de transporte y los resultados se 
presentan en la cuarta sección; y finalmente, en la quinta sección, un resumen de 
las contribuciones originales. 


2.1. Teoría de Floquet 


Cuando un sistema es sometido a una perturbación con dependencia temporal, 
la evolución temporal de los estados cuánticos está descrita por un Hamiltoniano 
dependiente del tiempo, a través de la ecuación de Schródinger dependiente del 
tiempo 

Edad =H(t)V. (2.1) 

Ot 

Para resolver esta ecuación usualmente se consideran dos regímenes. El primero, 
cuando la perturbación o la dependencia temporal del Hamiltoniano varían mu- 
cho más lentamente que los tiempos característicos del sistema, llamado límite 
adiabático; y el segundo cuando esta aproximación no se cumple y los tiempos ca- 
racterísticos de la dinámica del sistema son comparables a la dependencia temporal 
del Hamiltoniano. En el caso de una perturbación periódica en el tiempo, debemos 
comparar el período T' del Hamiltoniano con los tiempos característicos 7, por ejem- 
plo el tiempo de pasaje electrónico por el sistema. Cuando estos tiempos son del 
mismo orden, la perturbación es no adiabática. 

Consideramos una muestra de grafeno sometido al potencial electromagnético 
periódico de un láser coherente que se extiende en todo el dominio temporal y 
queremos conocer la solución del problema para todo tiempo, incluso tiempos largos. 
Para esto necesitamos una solución no adiabática y no perturbativa, dada por la 
teoría de Floquet. 


2.1.1. El Hamiltoniano de Floquet 


El teorema de Floquet, originalmente formulado por el matemático francés Achille 
Marie Gaston Floquet en 1883 para resolver ecuaciones diferenciales ordinarias, da 
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las soluciones a ecuaciones diferenciales lineales con coeficientes periódicos [44] 
y ha sido extendido al tratamiento de Hamiltonianos dependientes periódicamente 
del tiempo por Shirley y posteriormente por Sambe alrededor de 1970 [45, 46]. El 
formalismo de Floquet ofrece un método para obtener las soluciones al problema 
de autovalores de un Hamiltoniano dependiente periódicamente del tiempo, que es 
no perturbativa, válida para tiempos largos y más aún exacta. Este formalismo es 
análogo al teorema de Bloch, donde el potencial es periódico en el espacio y las 
soluciones al Hamiltoniano de Bloch están indexadas por el cuasimomento, que 
puede ser restringido a la primera zona de Brillouin. El teorema de Floquet hace uso 
de la periodicidad temporal del Hamiltoniano y los autovalores del Hamiltoniano de 
Floquet g., llamados cuasienergías, también pueden restringirse a la primera zona 
de Floquet (FZ). Si imponemos la periodicidad temporal del Hamiltoniano, el sistema 
y sus soluciones tienen simetría discreta de traslación temporal * . Es decir que, si 7 
es el periodo del Hamiltoniano tal que H(t + T”) = A(t), luego cualquier función de 
onda Y.,, que satisface la Ec. (2.1) debe ser invariante ante traslaciones temporales 
en T', a menos de una fase 


UT) =P. (2.2) 


Dado que ei"7 = 1, podemos utilizar la identidad e*7/? = ¿ile hT/h para 
acotar los valores posibles de e. al intervalo (-A02/2, h12/2) y definir £an = € +nhA0. 


En nuestro caso, la fase e*a7/% debe ser constante en el tiempo, ya que la 
dependencia temporal periódica del Hamiltoniano es la misma desde t > —oo y la 
solución no puede depender del tiempo inicial + desde el cual definimos las funciones 
de onda. Por lo tanto, la función et (1) debe ser invariante ante traslaciones 
temporales en un período del potencial T. Llamamos a esta función ¿(t), que cumple 


w(t) = TIME) p(t+T)=0(0), (2.3) 


tal que la Ec. (2.2) se satisface inmediatamente. Insertando la Ec. (2.3) en la ecua- 
ción de Schródinger dependiente del tiempo, Ec. (2.1), obtenemos 


1) 
-itat/h ip (,ritat/h 
H (e dalt)) = Ab (e da(t)) 
esti pa (e) = eeatin (sa dalt) + in, 0a(0)) 
A OS 
Luego, obtenemos el siguiente problema de autovalores 


E N in | Palt) = €u alt) 
Hr dalt) = En Dalt) > 


(2.4) 


2.1 Teoría de Floquet 


* Cuando un operador tiene una 
simetría discreta o continua y 
sus condiciones de borde tam- 


bién, 


entonces la solución al 


problema respeta estas sime- 


trías. 
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2.1 Teoría de Floquet 


donde ?+y es el Hamiltoniano de Floquet 


0) 
Hp =H-ih—. 
ij Ot 
Esta forma del Hamiltoniano de Floquet es la de un operador dependiente del tiempo 
que actúa sobre las coordenadas del sistema y sobre la variable temporal. Podemos 
arribar a una expresión estática del Hamiltoniano de Floquet, como se describe en 


la subsección siguiente. 


2.1.2. Hamiltoniano de Floquet en forma matricial 


El Hamiltoniano de Floquet es un operador que actúa sobre la variable temporal. 
Al hacer un cambio de base adecuado podemos obtener una ecuación análoga a la 
ecuación de Schródinger independiente del tiempo. Si expresamos los elementos de 
la Ec. (2.4) en la base Calida que expande las funciones temporales de periodo 
T =27/M obtenemos 


dalt) E y ¿ans H(t) - FATE, 


donde ¿”. es la componente n-ésima de la serie de Fourier de ¿,(t), al igual que 
H” lo es de H(t). Luego 


En Nr cimot - y go pr ¿m'+n)0t pn Y nhQ pr pnot 


m',n n 


Ea 97 = Y In, HPA OD INAAOG (2.5) 


cada = DATAN]. 
n 
En esta deducción sólo utilizamos la periodicidad temporal de ¿(t) y H(t), por lo 
que sus componentes de Fourier tienen coordenadas espaciales, o de espín, que 
expanden el espacio de Hilbert del problema para el caso estático (sin dependencia 
temporal). El espacio de Hilbert para el caso estático 7 producto el espacio de las 
funciones periódicas en el tiempo 7 es el espacio de Sambe R e T. Expandido en el 
espacio de Sambe [46], el Hamiltoniano de Floquet toma la forma matricial estática 
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E 7 
cuyos autovectores $,, = (...,9",90,01,02,...)” asociados a autovalores £. 
diagonalizan el Hamiltoniano de Floquet, de manera análoga a un problema de 
autovalores de un Hamiltoniano estático pero en un espacio extendido. 


2.1.3. Réplicas de Floquet 


Al resolver el Hamiltoniano de Floquet en el espacio de Sambe R e T obtene- 
mos un espectro de cuasienergías £,. La Ec. (2.6) es particularmente útil, ya que 
podemos comenzar a analizar el espectro de cuasienergías de Hp. Cada bloque 
diagonal H? + nh£) representa al sistema estático con energías desplazadas en 
nh0Q. Alo largo de esta tesis llamamos a cada uno de estos subespacios réplicas de 
Floquet, indexadas por el número n. La réplica de Floquet número n, representa al 
sistema estático excitado en energía y podemos interpretar que el sistema absorbió 
n fotones (emitió n fotones si n < 0). 

La primer peculiaridad es que podemos sumar ná(Q/ al Hamiltoniano de Floquet 
en la ecuación (2.6) y obtener un problema de autovalores idéntico, ya que en esta 
representación la matriz del operador de Floquet es infinita, aunque las réplicas 
de Floquet están desplazadas: a la que antes llamábamos m ahora se encuentra 
en el lugar n + m. De esta forma se manifiesta la libertad de gauge de definir 
Eam = Ea + nhS), tal que e, pertenece a la primera zona de Floquet (FZ) en el 
intervalo de cuasienergías (-A0/2, 0/2). 

El espectro de Floquet, entonces debe ser invariante ante traslaciones en nñQ 
(aunque los autovectores tengan sus componentes desplazadas). Esto nos dice dos 
cosas: primero que el espectro de cuasienergías es no acotado y segundo que el 
espectro de cuasienergías es periódico * . 

Si consideramos un Hamiltoniano estático 4% con un espectro A, sometido a 
una perturbación o un potencial externo, periódico en el tiempo, entonces 


H(t) =H2+AV(0, (2.7) 


tal que V(t + T') = V(t) y donde A representa un parámetro que podemos modi- 
ficar externamente para apagar y encender la perturbación. Podemos identificar 
al parámetro A con alguna cantidad física de nuestro problema. Si la perturbación 
se produce por proximidad del sistema original (estático) a un potencial periódico, 
entonces A es una función de la distancia; si el potencial es un potencial electro- 
magnético, entonces A es una función de la intensidad del campo; y en general 
representa la capacidad física de controlar la amplitud de la perturbación: para 
apagar el potencial A = 0 y encenderlo A = 1. 

Si la perturbación V(t) es completamente sinusoidal Y” = O (sin componente 
estática), entonces los elementos no diagonales de Hp serán las componentes 
de Fourier de V(t), tal que H” = AY” para n + 0. Si apagamos la perturbación 
poniendo A = 0 en la Ec. (2.7), el Hamiltoniano de Floquet toma la forma de una 
matriz diagonal por bloques, con copias del sistema no perturbado H0 desplazado 


2.1 Teoría de Floquet 


TDebido a que el espectro de 
Floquet es periódico podemos 
restringirlo al intervalo -h0/2 < 
£a < h0/2. Sin embargo, dado 
que los autovectores en el es- 
pacio de Sambe son diferentes, 
es útil mantener todos los au- 
tovalores para reconocer a qué 
réplica de Floquet pertenece ca- 
da componente de los autovec- 
tores $, . 
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hQ>>A 
Fig. 2.1: Réplicas de Floquet de 


un sistema irradiado a distintas fre- 


cuencias. El ancho de banda del 
sistema es A, cuando AM < A hay 
solapamiento de las réplicas en 
e = hQ/2. Además, si AQ < A/2 
también hay solapamiento en e = 
0). 


= 


ta proyección sobre la réplica 
cero cobra especial importancia, 
ya que la densidad de estados 
promediada en el tiempo coinci- 
de con la densidad de estados 
para el caso estático. Así recu- 
peramos el sistema no pertur- 
bado cuando la intensidad de la 
perturbación tiende a cero. 
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en energía en cada réplica. En el límite A = O todas las réplicas están desacopladas 
y el espectro de cuasienergías será el espectro de energías del sistema estático 
replicado a intervalos de Al). 

Como vemos en la Fig. 2.1, si la frecuencia del potencial externo 242 es mucho 
más grande que el ancho de banda del sistema estático Á, entonces las distintas 
réplicas no se solapan y no hay degeneración en el espectro de cuasienergías. 
Cuando disminuye la frecuencia al punto que A < Á, el espectro comienza a estar 
degenerado. Dado que el sistema original tiene simetría electrón-hueco, es decir que 
por cada estado (y, . con energía e. existe otro estado con energía —e.,, entonces a 
energía h0/2 y momento k habrá al menos dos estados degenerados; un estado 
perteneciente a la réplica O con energía e, = AQ/2 y un estado de la réplica 1 con 
energía —¿a + hM = h0/2. Cuando la perturbación V(t) contiene elementos de 
matriz que acoplan estos estados, al momento que encendamos el potencial externo 
(A > 0) la degeneración se verá removida generando un cruce evitado. 


2.1.4. Densidad de estados media 


Para entender cómo se relaciona el formalismo de Floquet y el espectro de 
cuasienergías con una medición real, podemos calcular qué valor daría un instru- 
mento que mida la densidad de estados. Suponemos que la frecuencia de oscilación 
del potencial es mucho más rápida que el tiempo de medición del instrumento. En 
este caso, debemos promediar la densidad de estados sobre un período 7'. Para 
calcular la densidad de estados promedio utilizamos el formalismo de funciones de 
Floquet-Green [47, 30], a partir del cual obtenemos 


: (Salva) 
ote) J dk Y lí, A] (2.8) 

En la Ec. (2.8) se pone pone de manifiesto la relación entre el promedio temporal 
de la función de onda de Floquet ¿(t) y la proyección sobre la réplica de Floquet de 
orden cero. El peso w”, de un estado de Floquet sobre la réplica n es simplemente 
we = (9u/o7). En el límite de intensidad cero (A = 0 en la Ec. (2.7)), donde las 
réplicas están desacopladas, el estado indexado por a solo pertenecen a una 


réplica, es decir que si pertenece a la réplica m, entonces wi” =1 y wí = 0,sin+m 
+ 


2.2. Efectos de la radiación en grafeno bulk 


El grafeno es un gran candidato a ser la materia prima de transistores más 
rápidos y más eficientes debido a su gran movilidad electrónica y capacidad de 
disipar el calor. Al ser un semiconductor de gap cero es necesario un mecanismo de 
encendido y apagado para controlar sus propiedades de transporte. La posibilidad 
de inducir un gap en grafeno mediante la incidencia de luz abre las posibilidades a 
aplicaciones de este material en transistores y nanodispositivos optoelectrónicos. 


2. Gaps en grafeno 2.2 Efectos de la radiación en grafteno bulk 


En esta sección se discutirá un trabajo previo a esta tesis, donde los autores Calvo, 
Pastawski, Roche y Foa Torres en el artículo «Tuning laser-induced band gaps in 
graphene» [32], proponen utilizar un láser en el rango del infrarrojo (h0 - 140meV) 
para abrir un gap apreciable y sintonizable a energías de la mitad de la frecuencia 
fotónica. 


Los efectos de luz circularmente polarizada fueron estudiados previamente, 
donde se predice un pequeño gap en el punto de Dirac. Para frecuencias en el 
rango del visible [30], una muy alta potencia (1W/um*) es necesaria para abrir 
un gap apreciable; y para frecuencias en el infrarrojo lejano (£(02=29 meV) [24] se 
analiza una juntura p-n y las fotocorrientes allí generadas. En cambio la propuesta 
de Calvo y col. [32], presentada en esta sección, explora frecuencias y potencias 
más accesibles, por ejemplo un láser de CO» de |=32 mW/um? de potencia, que 
podría abrir un gap apreciable de 20 meV en grafeno. Además se muestran los 
efectos sobre la densidad de estados efectiva de distintas polarizaciones, de lineal 
a elíptica y circular. Trabajos más recientes estudian los efectos de la radiación en 
grafeno donde muestran la posibilidad de modificar los puntos de Dirac [4.8], inducir 
polarización de valle en grafeno bicapa [49], estudian las propiedades de transporte 
del grafeno irradiado [25, 50]. 

Para tratar el grafeno iluminado por un láser se necesita un método que de 
cuenta de los efectos no-perturbativos y no-adiabáticos, por lo que podemos utilizar 
el teorema de Floquet [45, 46]. Consideramos una muestra limpia de grafeno (sin 
defectos), donde la luz aplicada es perpendicular al plano del grafeno, por lo que 
el potencial electromagnético es homogéneo en el plano, no produce scattering 
entre valles y podemos considerar los valles XK, y K_ por separado. En general solo 


utilizamos el valle K,, que denotamos XK, excepto cuando sea necesario distinguirlos. 


A bajas energías usamos el modelo k- p para entender las transiciones en el espacio 
recíproco (espacio k de momentos) generadas por la luz. En este modelo hacemos 
uso de la sustitución canónica k > k + ¿ A(t), para dar cuenta de los efectos de la 
radiación. El Hamiltoniano k - p que describe el valle K es 


Hrplt) = up oO -[p+eA(t)], (2.9) 


donde d son las matrices de Pauli y A(t) es el potencial electromagnético que 
describe a la radiación monocromática, coherente y polarizada 


A(t) = Az cos(Qt) E + Ay cos(0t + 0) Y - (2.10) 


Podemos tomar a las componentes 4,, = A, = Ay, de modo que nuestro parámetro 


que controla la polarización es la fase relativa ¿y entre las componentes de A(t). 


Para do = 0 la polarización es lineal en la dirección (4% + Y), para dp = 7/2 la 

polarización es circular y en valores intermedios la polarización es elíptica. 
Armados del teorema de Floquet, podemos prever qué sucederá cuando el 

potencial vector acople las distintas réplicas de Floquet. A bajas energías podemos 


e = +h0 

e=0 

e = —h0) 

E 
ñRQ << A 
k 


Fig. 2.2: Réplicas de Floquet pa- 
ra frecuencias h(2 mucho meno- 
res que el ancho de banda A. A 
bajas energías los conos de Di- 
rac describen la dispersión del sis- 
tema, sus réplicas se encuentran 
en +A), denotadas en rojo y azul. 
Los círculos verdes muestran los 
cruces en h£2/2 entre las bandas 
que pertenecen a réplicas aleda- 
ñas. 
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2.2 Efectos de la radiación en grafeno bulk 


describir al grafeno no perturbado como un cono de Dirac (Hoy = vpo - p en la 
Ec. (2.9)) replicado a energías nñ?2, como se esquematiza en la Fig. 2.2. Cuando 
la radiación acopla los estados que pertenecen a distintas réplicas, los efectos 
más importantes se observan cerca de las degeneraciones que ocurren a energías 
+h0/2. 


Los resultados del modelo k - p coinciden con los del modelo tight-binding, válido 
en un rango mayor de energías y es más preciso. El modelo tight-binding describe 
la red en el espacio real y en este caso la radiación actúa modulando los hoppings 
t; ¿ entre dos sitios ubicados en r; y r;, a través de la sustitución de Peierls 


2 0; 
tj — ti ¿(t) = texp (E pd ? A(t) dr) z (40 
0 YT 


donde dy es el cuanto de flujo magnético y t = 2,7eV es la amplitud de hopping 
típica para grafeno. La fase de los hoppings está modulada periódicamente en el 
tiempo y nuestro Hamiltoniano de red 


Hi(t) = DN estidl - Y tig(t) est! + Ho., (2.12) 


también será periódicamente dependiente del tiempo. 


Construimos el Hamiltoniano de Floquet-Bloch utilizando el formalismo del Apén- 
dice A.1 y podemos obtener la densidad de estados mediante las funciones de 
Floquet-Green. En la Fig. 2.3 se muestran los resultados de la densidad de estados 
promedio para distintas polarizaciones. En (a) la polarización lineal produce efectos 
a energías +A11%/2, que son la reducción en la densidad de estados media, aunque 
sin formar un gap. En el panel (b) y (c), la fase relativa 0 < Hp < 7/2 produce polari- 
zaciones elípticas y se observa un gap tanto en el punto de Dirac como a energías 
+h0/2. El máximo ancho del gap se produce para polarización circular, cuando 


do = 7/2. 


En la Fig. 2.3 se puede observar que es la componente circular de la polarización 
la que abre los gaps, por lo que las polarizaciones elípticas, que pueden ser des- 
compuestas en una parte circular y una lineal, abren un gap que es máximo cuando 
la polarización es completamente circular. Se puede estimar el ancho de este gap 
notando que tanto el cono de Dirac como la polarización circular son isotrópicos 
en el sentido de que todas las direcciones en el plano (k,,, k,) son equivalentes 
y podemos elegir una dirección arbitraria en particular. En esta dirección el cruce 
evitado entre la réplica +1 y O abre un gap de Ay.-9/2v. = € Ao UF 


Además, también se observa el gap estudiado por Oka y Aoki en la Ref. [30], 
que se abre en el punto de Dirac. Este gap es de segundo orden (más pequeño) ya 
que su ancho está dado por 


(e Ay up)? 
Aro = 12 == 
k=0 v2 AO > 
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2.3 Efectos de la radiación en nanocintas 








y depende cuadráticamente de la intensidad del potencial vector. 

Se puede observar que cuando la intensidad es grande aumenta la densidad de 
estados mínima dentro de los gaps, generando «sub-gaps» de orden superior. Estos 
efectos serán tratados con mayor profundidad en los Capítulos 3, 4 y 5, aunque 
rápidamente se puede entender que los efectos provienen de los cruces entre las 
bandas de las réplicas con n + O que se cruzan lejos de k = 0 en la Fig. 2.2. 


2.3. Efectos de la radiación en nanocintas de grafeno 


En la sección anterior se analizan los efectos de la radiación de distintas pola- 
rizaciones en una muestra bulk, es decir, infinita y sin tener en cuenta los bordes. 
En este caso se observan gaps cuando la polarización es circular y no se observan 
gaps cuando la polarización es lineal. En esta sección analizamos los efectos de la 
radiación de polarización lineal en nanocintas de grafeno (GNR's) y mostramos que 
es posible abrir gaps o cerrarlos si cambiamos la dirección de polarización lineal. La 
dirección de la polarización necesaria para abrir un gap depende fuertemente del 
tipo de borde de la nanocinta, lo que se puede explicar mediante las transiciones 
intra-modo e inter-modo que induce la radiación. Estos son los primeros resultados 
originales durante el desarrollo de esta tesis. 

El grafeno presenta naturalmente bordes preferencialmente de dos tipos: zigzag 
y armchair. Utilizando como variables el tipo de borde y la dirección de la polarización 
veremos cómo se pueden abrir gaps en nanocintas de grafeno. La posibilidad de abrir 


Fig. 2.3: DOS para una monoca- 
pa de grafeno iluminada por un lá- 
ser de A? = 140meV, polarizada 
(a) linealmente, (b) elípticamente 
con do = m1, (c) elípticamente con 
do = TÍ y (d) circularmente. La 
línea sólida negra es para inten- 
sidades de / = 32mW/um? mien- 
tras que la línea roja punteada pa- 
ra I = 130mW/pum?. Para polari- 
zación lineal se comparan los re- 
sultados con los del modelo tight- 
binding (en azul a trazos y violeta a 
trazos y puntos). Figura adaptada 
de [32]. 
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Fig. 2.4: (a) Red hexagonal de 
grafeno, en gris se muestran las 
direcciones de corte para obtener 


bordes de tipo zigzag y armchair. 


(b) aGNR de N filas de sitos de 


ancho. (c) zZGNR de N filas de si- 


tos de ancho. (d) Ejemplo de un 
borde arbitrario en una muestra de 
grafeno. Figura adaptada de [4]. 
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gaps incluso con polarización lineal diferencia las nanocintas del comportamiento 
para grafeno bulk. Estos gaps muestran conductancia cero, incluso con polarización 
lineal en los cálculos de transporte, lo que hace a las nanocintas un candidato ideal 
para diseñar dispositivos optoelectrónicos que hagan uso de un mecanismo de 
encendido y apagado del transporte electrónico. 


2.3.1. Nanocintas de grafeno 


Consideraremos dos tipos de nanocintas de grafeno, las de bordes de tipo zigzag 
(ZGNR) y las nanocintas con bordes de tipo armchair (aGNR). En la Fig. 2.4 (a) 
se muestran las direcciones de corte (a) de ambos tipos de cintas. Un esquema 
de una cinta N-aGNR se muestra en el panel (b), con su celda unidad en líneas 
de puntos. El ancho W,, de la cinta es proporcional al número de filas de dímeros 
W,¿(NW) = V3/2(N - 1)a. En el panel (c) se muestra una cinta N-zGNR, con su celda 
unidad y su ancho W., depende del número de filas zigzag W.(N) = 3/2(N — 1)a. 
Una dirección de corte arbitraria, panel (d), es una combinación de tramos zigzag y 
armchair que ocurre cuando no se puede controlar el tipo de borde. 


(a) dirección armchair (b) 


d 




















Al confinar la función de onda electrónica dentro de una nanocinta, esperamos 
que los modos electrónicos dependan del ancho de la cinta y también de las con- 
diciones de borde. Las condiciones de borde son de masa infinita, pedimos que la 
función de onda se anule en sitios virtuales justo por fuera del borde. La solución 
analítica se muestra en el Apéndice A.2, siguiendo el formalismo presentado por 
Marconcini y Macucci en la Ref. [51]. Como consecuencia del confinamiento el 
vector de onda transversal está cuantizado en k, = k,, y puede ser función de la 
componente paralela k, = k, (k). 

En el caso de los bordes zigzag, k, sólo depende cerca de k¡ - 0 y para 
ky > 1/W(V) se hace independiente de manera asintótica. Para el caso de la 
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cinta armchair, k, es independiente de kj en todo su rango. En ambos casos 
kn(ky) define superficies de corte en la zona de Brillouin, cuya intersección con la 
relación de dispersión para el grafeno bulk e(k; , k, (ky )) debemos proyectar sobre 
la dirección ky y tener en cuenta todos los modos n para obtener la dispersión de 
una nanocinta de grafeno. 

Además de los modos extendidos, en las cintas zigzag existen modos que no es- 
tán en la relación de dispersión bulk. Estos modos no pertenecen completamente a 
la zona de Brillouin ya que son modos evanescentes, caracterizados por Im(k,) + 0 
y están localizados en cada borde de la cinta zigzag. Cerca de k; =- 0, la componente 
real de k, es distinta de cero y los modos son metálicos (dispersivos). Sin embargo 
rápidamente al crecer k > O la componente real de k, decae exponencialmente, 
junto con ella la velocidad de este modo electrónico y crece la componente imagina- 
ria, es decir que se localizan completamente y son prácticamente no dispersivos. 
Estos modos localizados reproducen a los modos que conectan los valles XK, del 
modelo tight-binding, como se muestra en la Fig. 2.5. 

En las cintas armchair todos los modos permitidos son una combinación de 
modos con k,, y —k, del grafeno bulk y cada k,, define planos rectos de corte en 
la zona de Brillouin. Cuando un plano de corte pasa por los puntos K, la cinta es 
metálica, esto sucede siempre que la cinta sea de ancho N = 3 * p- 1, con p un 
número entero. Cuando el número de filas de dímeros que definen el ancho de la 
cinta, no es NV =2 mód (3), la cinta es semiconductora, como se muestra en la 
Fig. 2.6. 

Por convención, tomaremos las direcciones de corte mostradas en la Fig. 2.7 
(a), siendo a = a; + az y z = a; — as. Estas direcciones son las de simetría de 
desplazamiento espacial y debemos plegar la zona de Brillouin como muestra la 
zona gris en (b). En la Figs. 2.5 y 2.6 se muestran las relaciones de dispersión 
para cada tipo de cinta, que son la relación de dispersión para el grafeno bulk, 
pero plegada como la zona de Brillouin. En ambos casos las modificaciones más 
significativas que produce el tipo de borde se dan a bajas energías: (i) las zigzag las 
cintas son siempre metálicas, con estados localizados que conectan los valles K,; 
(ii) las cintas armchair no tienen estados localizados y según el ancho de la cinta 
son metálicas o semiconductoras. 


2.3.2. Cambios en la densidad de estados 


En las nanocintas de grafeno los efectos son más particulares y pueden ser más 
pronunciados, producto de un interjuego entre las transiciones propiciadas por la 
radiación y los modos transversales electrónicos permitidos. Los modos permitidos 
estarán dados por el tipo de borde y el ancho de la nanocinta y las transiciones 
permitidas o no permitidas tienen sus casos extremos cuando la polarización lineal 
es paralela o perpendicular al borde de la nanocinta. 

Para describir el sistema, utilizamos el Hamiltoniano tight-binding donde los 
hoppings están modificados por la sustitución de Peierls, Ec. (2.12), para obtener 





Fig. 2.5: Relación de dispersión 
para una cinta 20-zGNR. 





Fig. 2.6: Relación de dispersión 
para una cinta 19-aGNR semicon- 
ductora. 
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Fig. 2.7: (a) Red hexagonal de 
grafeno con sitios de tipo A y B 
denotados con círculos llenos y va- 
cíos respectivamente. Los vecto- 
res á y z son los vectores de red 
y a la vez las direcciones de cor- 
te para obtener un borde armchair 
o zigzag. (b) Primera zona de Bri- 
llouin para una red hexagonal bulk 
y en gris para una red con borde 
armchair o zigzag. (c) y (d), relacio- 
nes de dispersión para una cinta 
con bordes armchair y zigzag, res- 
pectivamente. Figura adaptada de 


[4]. 


Fig. 2.8: En el panel de la izquier- 
da se muestra la relación de dis- 
persión para una cinta armchair 
muy angosta. Los modos de las 
réplicas +1 (en líneas punteadas 
en rojo y azul) se cruzan con los 
modos de la réplica O (líneas con- 
tinuas en negro). A la derecha un 
esquema que muestra las transi- 
ciones posibles. Si la transición es 
entre un modo q y uno distinto q”, 
la transición es inter-modo, cuan- 
do ambos modos tienen el mismo 
q, la transición es intra-modo [42]. 
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la DOS media como muestra la Ec. (2.8). La celda unidad crece con el ancho de 
la cinta y por lo tanto el costo computacional de aplicar el método recursivo de las 
funciones de Floquet-Green como se detalla en [47, 43]. Además, el número de 
réplicas de Floquet consideradas en el cálculo es un factor que multiplica el tamaño 
de la matriz celda unidad. Para intensidades bajas basta considerar +2 canales 
fotónicos para dar cuenta de las modificaciones en la DOS media. 

Según el tipo de borde y la dirección de la polarización o bien las transiciones 
inter-modo o las intra-modo están permitidas (ver Apéndice A.2), como se muestra en 
la Fig. 2.8, mientras que para polarizaciones intermedias los dos tipos de transiciones 
ocurren. 





Se observan reducciones en la densidad de estados media, que ocurren exacta- 
mente a h£2/2 para las transiciones intra-modo o bien en posiciones aledañas para 
las transiciones inter-modo. Esto se debe a la simetría electrón-hueco dado que 
cuando existe un modo ¿(e, q) con energía e, existe también ¿(—e, q) con energía 
—e. Esto sucede para todas las réplicas de Floquet, por lo tanto las transiciones 
intra-modo se producen cuando el modo q de la réplica n = O se cruza con el modo 
-q de la réplica n = 1 y el cruce (evitado) entre ellos se produce exactamente en 
hQ0/2. En cambio si las transiciones son inter-modo, cuando un modo q de la réplica 
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n = 0 se cruza con un modo —q/ de la réplica n = 1. La energía a la que se cruzan 
estos modos q y -q' está desplazada de h£?/2. 

Las cintas armchair tienen transiciones intra-modo cuando la polarización se 
encuentra en la dirección £ paralela al borde. En la Fig. 2.9 se indica con una 
flecha la energía a la que ocurren las degeneraciones (4/2 = 70 meV) y se ven las 
reducciones en la DOS de una cinta 83-aGNR semiconductora para polarización 
(a) en la dirección %, (b) en la dirección y y (c) circular, donde ambos tipos de 
transiciones ocurren. En el panel (a) se puede observar que el la cinta sin radiación 
tiene dos modos en electrónicos h(2/2. Como ambos modos se cruzan con su par 
electrón-hueco de la réplica de Floquet siguiente, los dos modos desarrollan un cruce 
evitado y se forman dos gaps concéntricos. Los anchos de los gaps son distintos 
porque las transiciones son sensitivas a la componente transversal k, del modo 
electrónico que acoplan. En el panel (b) la polarización en y inhibe las transiciones 
intra-modo y propicia las transiciones inter-modo, por lo que se observan reducciones 
en la DOS «satélites» a h02/2. Según la intensidad de la radiación puede que se 
desarrollen verdaderos gaps o no. El gap a izquierda (energía menor que h(?/2) cruza 
una singularidad de van Hove —que indican el comienzo de un modo electrónico— 
y se abre un gap en el primer modo electrónico. En el panel (c) podemos ver los 
efectos combinados de la polarización en £ junto con los efectos de la polarización 
en y. Para cintas pequeñas e intensidades bajas las modificaciones producidas por 
la polarización circular no difieren cualitativamente de las de una polarización lineal 
en la dirección E + y. Sin embargo, como vimos en la Sección 2.2, en el límite bulk 
ambas polarizaciones producen efectos muy distintos. Mientras que la polarización 
circular abre un gap, la polarización lineal no abre un gap. 


Fa) N=83, polarización en "x" 
[——] sin radiación 

p 1mW/um? li 

| 47 10 mW /um? 


Fig. 2.9: Densidad de estados 
media para una cinta aGNR con 
N = 83 para polarización (a) lineal 
en la dirección 2, (b) en la direc- 
ción y, y (c) circular [43]. 


DOS (u.a.) 











Todas las cintas aGNR semiconductoras, desarrollan gaps a energías +h12/2 
cuando la polarización incidente tiene una componente en k.,. Sin embargo, cuando 
la cinta es aGNR metálica, como se muestra en la Fig.2.10, el modo metálico muestra 
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2.3 Efectos de la radiación en nanocintas 





Fig. 2.10: Densidad de estados 
media para una cinta aGNR con 
N = 52 metálica irradiada (AM = 
140meV, P = 1mW/um?) con po- 
larización lineal en la dirección E 
(línea azul punteada) y (b) en la di- 
rección y (línea roja continua) [43]. 


Fig. 2.11: Densidad de estados 
media para una cinta zGNR con 
N = 250 metálica para polariza- 
ción lineal en la dirección £ (lí- 
nea azul punteada) y (b) en la 
dirección y (línea roja continua) 
para una radiación de frecuencia 
hQ = 140meV y una potencia de 
P = 1mW/¡m? [43]. 
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gran robustez frente a la radiación de polarización lineal en cualquier dirección. Si 
este es el único modo presente a energía e = h(02/2 no se observa un gap en la 
DOS y sólo ocurren modificaciones leves. En el caso en que más de un modo 
esté permitido a energía h(2/2 en el sistema no perturbado, entonces habrá una 
reducción en la DOS, excepto para el modo metálico de la cinta. 
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polarización en "x" 
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Las cintas zigzag pueden ser analizadas de manera análoga a las armchair y en 
la Fig. 2.11 se observan efectos similares, con la gran diferencia de que ahora la 
polarización en y es la que abre gaps en h£2/2 y al polarización en £ abre los gaps 
satélites. A diferencia de las cintas armchair, las cintas zigzag son siempre metálicas 
y el estado metálico no es robusto frente a polarización lineal. Se observa un gap 
para polarización perpendicular a la dirección del borde, en la Fig. 2.11 para una 
cinta 250-zGNR. Cuando la polarización está en la dirección £ paralela al borde las 
reducciones en la DOS media son a energías aledañas a e = 70meV. 
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En el análisis presentado hasta aquí, las nanocintas aGNR y zZGNR muestran 
propiedades bien diferenciadas en cuanto a la modificación de la DOS efectiva, 


2. Gaps en grafeno 2.4 Propiedades de transporte 





y un comportamiento prácticamente opuesto respecto a la dirección de la polari- 
zación. Sin embargo, sabemos de los resultados bulk que en el límite de ancho 
muy grande, una cinta zigzag y una armchair deben aproximarse. Comprobamos 
numéricamente este límite cuando el número de bandas que se cruzan en h(2/2 es 
grande. Para mantener un costo computacional accesible, aumentamos la frecuencia 
de la radiación, manteniendo el parámetro de acople 4p7/A£ constante, para obtener 
cualitativamente el límite bulk en una cinta 300-aGNR. Lo que vemos en la Fig. 
2.12 es que la polarización lineal, tanto en la dirección (a) E, (b) y y (c) t + Y, 
produce reducciones en la DOS que se acumulan alrededor de 2(2/2. En cambio, la 
polarización circular (d) produce un gap en todos los modos semiconductores de 
esta cinta armchair y cabe notar que el estado metálico sigue mostrando robustez 
incluso para una polarización no lineal e intensidad muy grande. Como veremos en 
el capítulo siguiente, este estado es mucho más robusto de lo esperado y ocurre en 
cualquier tipo de borde. 








2.4. Propiedades de transporte 


Las modificaciones en la DOS promedio para las GNR's irradiadas, se reflejan 
en sus propiedades de transporte. Para describir los efectos del potencial periódico 
en el tiempo podemos utilizar las funciones de Green de no equilibrio y el formalismo 
de Keldysh [52] o equivalentemente el formalismo de scattering de Floquet. En la 
Fig. 2.13 podemos ver los resultados del cálculo de transporte a dos terminales para 
una cinta aGNR. En el caso de la cinta metálica, paneles (a) y (b), la polarización 


Fig. 2.12: DOS de una cinta 
aGNR con N = 300 metálica, 
en gris para sistema sin pertur- 
bar y en naranja la DOS media 
con radiación de frecuencia A = 
2700meV y polarización lineal (a) 
en la dirección 2, (b) en la direc- 
ción y, (c) en la dirección £ + y y 
(d) polarización circular [43]. 
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2.5 Contribuciones originales 





Fig. 2.13: Transporte a dos termi- 
nales para un cinta armchair me- 
tálica, de longitud N = 130 para 
(a) y (b) y una cinta armchair se- 
miconductora de longitud N = 129 
para (c) y (d). En gris se muestra 
el perfil de conductancia para la 
muestra no irradiada, en naranja 
para una radiación de 1mW/um? y 
las líneas de puntos rojos para una 
radiación diez veces mayor. En los 
paneles (a) y (c) la polarización de 
la luz está en la dirección £ y en 
los paneles (b) y (d), la polariza- 
ción está en la dirección y [42]. 











k(a7) ”- DOS(u.a.) 


Fig. 2.14: Esquema donde se 
muestran a izquierda las bandas 
de una cinta armchair muy peque- 
ña y los cruces evitados origina- 
dos de las transiciones permitidas. 
A la derecha, la DOS de color azul 
(oscuro) se obtiene cuando la po- 
larización es paralela a la cinta y 
la DOS naranja (claro) cuando es 
perpendicular. 
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lineal produce modificaciones en energías h(2/2 o aledañas pero no logra anular 
la conductancia. En el caso de una cinta semiconductora, la radiación lineal (c) en 
la dirección de E, produce un gap llevando la conductancia a cero; y la radiación 
lineal perpendicular (d) no logra reducir a cero la conductancia. Así podemos en- 
cender y apagar la conductancia de una cinta semimetálica a energía h£2/2 o bien 
encendiendo y apagando la radiación o cambiando la dirección de la polarización. 





(a) N=130, polarización en "x" 


| sin radiación — [EW 74 
FS 1 mw /um? 


asa 10 MW lu? 







"mn 


b) N=130, polarización en "y 





(c) N=129, polarización en "x" 


G (2e*/h) 








(c) N=129, polarización en "y" 
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2.5. Contribuciones originales 


= Desarrollo de códigos propios para calcular las funciones de Floquet-Green con 
el método recursivo, utilizadas para obtener la DOS promedio y también para 
realizar los cálculos de transporte a dos terminales. Los resultados obtenidos 
fueron comparados con códigos desarrollados paralelamente por Hernán Calvo y 
Luis Foá Torres para, de esta manera, verificar la corrección de los datos. Los 
códigos propios mostraron ser comparables en precisión y eficiencia. 


= Estudiamos cómo la radiación de polarización lineal induce transiciones intra- 
modo o inter-modo en GNR's según el cómo se acople la radiación a los modos 
permitidos por las condiciones de contorno de la cinta, como se muestra en la 
Fig. 2.14. 


= Obtenemos gaps en cintas armchair irradiadas con polarización lineal. Según la 
dirección de la polarización los gaps (o reducciones) en la DOS media ocurren 
a la mitad de la energía fotónica o en energías aledañas. Según el ancho de la 


2. Gaps en grafeno 2.5 Contribuciones originales 





cinta encontraremos gaps o reducciones en la DOS media. En cintas metálicas 
el estado metálico es robusto a polarización lineal. 


= Obtenemos gaps en cintas zigzag irradiadas con polarización lineal. Según la 
dirección de la polarización los gaps y reducciones en la DOS media ocurren a 
la mitad de la energía fotónica o en energías aledañas, para cintas de todos los 
anchos. 


= Demostramos que podemos utilizar la dirección de la polarización como meca- 
nismo de encendido y apagado en nanodispositivos basados en nanocintas de 
grafeno, para los dos tipos de cintas más abundantes, armchair GNR y zigzag 
GNR. 
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Abstract 


La radiación de polarización circular puede inducir gaps en la densidad de estados 
promedio de una muestra de grafeno bulk. Sin embargo, cuando la muestra tiene bordes, 
existen estados que persisten dentro del gap. Estos estados prueban ser quirales, es decir 
que los que se propagan en una dirección están localizados en un borde de la muestra, 
mientras que en el otro borde se propagan en la dirección opuesta. La quiralidad hace que la 
conductancia sea muy robusta frente a defectos estructurales y es una cualidad distintiva de 
los aislantes topológicos. Las propiedades topológicas no triviales pueden ser inducidas por 
la radiación y están caracterizadas con un invariante topológico, que muestra que el grafeno 
irradiado con luz circularmente polarizada es un aislante topológico de Floquet. 


En este capítulo se muestran los principales resultados de esta tesis: la emer- 
gencia de estados quirales en grafeno irradiado [34], su caracterización analítica 
y sus propiedades topológicas [53]. Otros antecedentes del estudio de estados de 
borde en grafeno irradiado pueden encontrarse en [54, 55, 56] En la Sección 3.1 
analizamos los resultados preliminares, acerca de la densidad de estados media 
y relación de dispersión efectiva para cintas de grafeno al incrementar el ancho de 
la cinta y disminuir los efectos de tamaño finito. Luego, en la Sección 3.2, demos- 
tramos analíticamente la existencia estados localizados en el borde zigzag de una 
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Fig. 3.1: Relación de dispersión 
efectiva para nanocintas de gra- 
feno de borde zigzag. En color se 
muestra la contribución de los es- 
tados a la DOS media. En los pa- 
neles se muestran los resultados 
para (a) una cinta de 21,3nm, (b) 
una cinta de ancho 42,6nm y (c) 
una cinta de ancho 426nm. El gap 
dinámico aparece más claramente 
para cuando el ancho de la cinta 
crece. Figura adaptada de [34] 
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3.1 Estados de borde 


muestra de grafeno irradiado. En la Sección 3.3 mostramos que la presencia de los 
estados no dependen del tipo de borde y además, sus propiedades de transporte 
son robustas frente a defectos. 


3.1. Estados de borde de Floquet inducidos por radiación 


El grafeno bajo la incidencia de radiación desarrolla gaps en el interior del 
material a energías conmensuradas con la mitad de la frecuencia fotónica nh(2/2. 
Cuando el material está confinado en una dirección, los modos normales que se 
propagan en esta dirección sufren modificaciones particulares según el tipo de borde 
y la dirección de la polarización, como se describió en la Sección 2.3. Muchos de 
estos efectos provenían de la interacción entre los efectos de tamaño finito y las 
transiciones entre distintos modos inducidas por la radiación. Además, mostramos 
que en las cintas armchair metálicas existen estados robustos frente a la radiación 
lineal y circular, incluso en el límite de una cinta muy ancha. En cambio en las cintas 
zigzag podíamos abrir gaps con polarización lineal y no se observaban estados 
robustos. ¿Son los estados robustos una particularidad de las cintas armchair 
metálicas, o existen para cintas de cualquier terminación? Vemos a continuación 
qué sucede con una cinta zigzag irradiada con polarización circular cuando el ancho 
de la cinta aumenta. 


3.1.1. Primeras evidencias 


Para entender las transiciones entre modos es útil graficar la relación de disper- 
sión efectiva. Ésta consiste en obtener la densidad de estados media p"(k, e) de los 
modos f, en la dirección de la cinta. Utilizamos las funciones de Floquet-Green [47, 
43] para obtener p"(k, e), esta densidad de estados resuelta en k.. nos muestra la 


relación de dispersión efectiva. 


0.5 





En la Fig. 3.1 exploramos los efectos de la radiación en una nanocinta de grafeno 
zigzag, para distintos anchos. Cuando el ancho es pequeño, paneles (a) y (b), 
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los efectos de tamaño finito, junto a las transiciones entre modos inducidas por 
la radiación (ver Sección 2.3), generan modificaciones en la densidad de estados 
media y en la relación de dispersión efectiva, donde los cruces evitados forman un 
gap en el espectro. Sin embargo, cuando el ancho de la cinta aumenta, panel (c), los 
efectos sobre la DOS media se acercan al límite bulk, donde claramente se observa 
un gap en las bandas de la nanocinta de grafeno, con la excepción de los estados 
que cruzan el gap. Éstos comienzan a insinuarse en el panel (b) y convergen a dos 
estados (por valle K,) dispersivos cuando el ancho de la cinta aumenta, panel (c). 
Los estados permitidos están dados por las transiciones permitidas por la radiación 
en una cinta zigzag, aunque la presencia de los estados de borde no es accidental, 
y como veremos más adelante, las razones para la existencia de estos estados son 
más firmes y tienen sus bases en la topología del sistema irradiado. 

Podemos evidenciar en la Fig. 3.1 (c) la primera diferencia entre los estados 
dentro del gap y los estados que forman las bandas en grafeno: El número de modos 
extendidos crece con el ancho de la cinta (estos estados forman las bandas en 
grafeno y corresponden con los modos permitidos en una muestra bulk), mientras 
que el número de estados dentro del gap rápidamente converge a un número 
finito (dos estados por valle K,) que no crece con el ancho de la cinta. Podemos 
preguntarnos qué otras diferencias hay y en la Sección siguiente probamos si los 
estados dentro del gap se tratan de modos extendidos o de estados localizados en 
los bordes de la muestra. 


3.2. Cálculos analíticos 


Queremos entender cuáles son los modos permitidos tanto en el bulk del grafeno 
irradiado como en una muestra terminada en bordes zigzag. Para hacer esto co- 
menzamos analizando el Hamiltoniano de Floquet-Bloch para bajas energías, luego 
obtenemos las propiedades de los estados de Floquet en en el bulk del material y 
luego deducimos las soluciones del Hamiltoniano de Floquet imponiendo condicio- 
nes de borde zigzag. La comparación entre los resultados analíticos del modelo k - p 
con los resultados numéricos del modelo tight-binding muestran muy buen acuerdo, 
tanto para cintas de grafeno como para una muestra semi-infinita. 

La radiación que actúa sobre la muestra de grafeno puede ser descrita con el 
potencial vector A(t) = Re(Age'*) que describe el campo electromagnético que 
incide perpendicularmente a la muestra de grafeno, por lo que no tiene dependencia 
espacial en el plano (x, y). La dirección de la polarización está dada por la dirección 
del vector Ay, que para polarización circular es Ay = Ap(£ + i9). Podemos utilizar 
la sustitución canónica p > p+eA en el modelo k - p y obtenemos el Hamiltoniano 
dependiente del tiempo 


H(t) = hop (s me -A(0) 0. (3.1) 


3.2 Cálculos analíticos 
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Fig. 3.2: Esquema donde se apre- 
cian los cruces entre réplicas más 
relevantes. Los círculos denotan 
los cruces entre réplicas O y +1 
y los triángulos los cruces entre 
las réplicas 1 y —1.Los rectángu- 
los en k. = 0 indican el gap en 
e = nh? que proviene de la renor- 
malización de las bandas [53]. 
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Construimos a partir del Hamiltoniano de la Ec. (3.1) el Hamiltoniano de Bloch, que 
en el espacio de Sambe toma la forma matricial 


hupk-o+h91I “E Apo_ 0 


20 
Hp = ds ZE Ayo, hvpk-0 ES Apo “sl, (3.2) 
Ze 


0 “ES Ando, hvpk-o-hQI 


donde o, = (0, +i0,). De modo que cambiamos la representación temporal de 
la perturbación, Ec. (3.1), por una representación estática en el Hamiltoniano de 
Floquet infinito, Ec. (3.2). 

El Hamiltoniano de Floquet es una matriz tridiagonal por bloques, donde los 
bloques diagonales representan las réplicas de Floquet y los bloques fuera de la 


diagonal representan el acople entre réplicas. Utilizaremos el parámetro de acople 
_ eupAo 

n= o 

están representadas las réplicas n = -1,0,+1. Los cruces entre réplicas de primer 





y trabajaremos en el régimen de acople débil (y << 1). En la Fig. 3.2 


orden en el acople están denotados por círculos, los cruces de segundo orden por 
triángulos y los cuadrados representan la renormalización de las bandas debido a 
los procesos de absorción y reemisión de fotones, por lo que se abre un gap también 
de segundo orden (O(1?)). 


3.2.1. Gaps y estados de grafeno bulk 


En el régimen de acople débil, nos centraremos en los procesos de primer orden 
en 7, que son los más significativos. Queremos describir el espectro de Floquet y los 
estados para energías cercanas a hA0/2 (círculos superiores en la Fig. 3.2) donde 
se cruzan las réplicas n = 0 y n = 1, por lo que podemos truncar el Hamiltoniano de 
Floquet y restringirlo a estos dos subespacios 


hopk-o+hQ1 “EL Apo. 
=P pa 2 a (3.3) 


“ES Ayo hupk-0 


El bulk del grafeno presenta gaps en la DOS media, inducidos por radiación de 
polarización circular (Sección 2.2). En la Fig. 3.3 están representados los estados 
en el espacio k que pertenecen al espectro de Floquet y cuya contribución a la DOS 
media se indica por la transparencia. 

Para describir los estados del espectro de cuasienergías reescribimos el Hamil- 
toniano de Floquet truncado a las réplicas n =0y n= 1, 


22 > hork 0 0 

% hupk 40 BREA 0 AO 

a me > de (3.4) 
0 e A => hupk_ 2 
0 0 hupk, -—2 
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(b) 


AD 





cuyas energías gy, tienen autofunciones 


ik- k, k, k, ka T 
ra (5) =e "dt a a nr > (3.5) 


donde e es la amplitud de la función de onda en la réplica ¿ en cada subred. 

Cuando la amplitud del campo electromagnético es cero 4y = 0 (apagamos la 
radiación), encontramos que las energías son e; , = +hupk y e] , = hQ+hupk. Dos 
de estas ramas están degeneradas en ky = 0/2vp, donde ef y, = €] 4, - Al encender 
la radiación Ay > 0 la degeneración en h(/2 se transforma en un cruce evitado, 
abriendo un gap en el espectro de Floquet. Podemos encontrar todas las soluciones 
al Hamiltoniano de la Ec. (3.4) introduciendo el vector de la Ec. (3.5) y resolviendo el 
sistema, aunque en nuestro caso es útil restringir aun más el problema, considerando 
sólo las ramas que se cruzan en 20/2. Considerando sólo este subespacio podemos 
ver que las energías de las ramas cuasi-degeneradas son 





hQ 
Es = (+ 12); pa = (1 k/ko)? +07, (3.6) 
esquematizadas en la Fig. 3.3 
Las autofunciones del Hamiltoniano completo son 
- sin (5/2) ch 
a 1 sin (0/2) 
ik-r k 
alr)=e"" —= X ; 3.7 
cos (+j/2) 


donde A es el área de la muestra de grafeno, 0, es el ángulo entre k y el eje x y el 
ángulo p;, es una medida de cuán cerca está h del gap, 


(3.8) 





3.2 Cálculos analíticos 


Fig. 3.3: (a) Grafeno bulk irra- 
diado. El campo electromagnéti- 
co es homogéneo en el plano del 
grafeno dada la incidencia normal 
de la radiación. (b) Esquema que 
muestra la relación de dispersión 
tridimensional efectiva para k cer- 
ca de K'x+ a bajas energías. La 
transparencia de la superficie in- 
dica el peso sobre la réplica 0. Fi- 
gura adaptada de [53] 
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siendo pj, - 7] lejos de k = ko y pz, = +7/2 en el límite k > k;. La velocidad instantá- 
nea de los estados es el valor de expectación del operador v = vpo que podemos 
descomponer en componentes paralelas (v,) y perpendiculares (v,) respecto del 
vector de onda k, y en nuestro caso 


(Wy)h+ = UF COS Pi 
1=k/lo 
y HAN 


dl (3.9) 
(01): = UF sin p; sin 0t 


= A sin Qt. 
Mz 


Podemos ver que v, oscila con la frecuencia de la radiación, mientras que uy es 
constante. Los valores de la velocidad promediadas en el tiempo 


((01))x2 = UF COS pj, 





(3.10) 


donde (()) denota el promedio temporal del valor de expectación. La Ec. (3.10) 
indica que las autofunciones se propagan (en promedio) en la dirección de k, es 
decir, paralelos al vector de onda. Además observamos que el valor de vj tiene el 
signo de , para k < ko y el signo opuesto para k > ko, llegando a velocidad cero 
cuando k = kg. También podemos calcular la componente (0,)x. = sin pj cos (2, 
lo que nos dice que el pseudoespín precesa alrededor de la dirección k£ con la 
frecuencia de la radiación y la amplitud de oscilación es máxima en k = ky 


07. = hecos pi + (2 x kesin y; sin 0t + Z cos Mt). (8.11) 
k k 


3.2.2. Estados quirales de borde en grafeno con terminación zigzag 


Los bordes de grafeno de tipo zigzag se caracterizan por que la primera fila de 
sitios pertenece a una misma subred mientras que la última fila de sitios es del tipo 
opuesto. 

Supongamos que cortamos nuestro plano de grafeno por el eje y y analizamos 
las soluciones para el semiplano y > 0. Las condiciones de borde de masa infinita 
cuando la primera fila en y = ay es de sitios de tipo A exigen que la función de onda 
se anule en el borde y = 0, donde están los sitios virtuales que pertenecen a la 
subred B. Esta condición de borde es particularmente útil ya que para satisfacer 
estas condiciones de contorno, no es necesario acoplar estados pertenecientes 
a los distintos valles y podemos tratar los valles K, por separado. La invarianza 
traslacional en la dirección x, nos dice que esperamos soluciones proporcionales 
a e*»", mientras que los estados no se propagan más allá del borde y esperamos 
soluciones tengan velocidad v, = O. Una forma de obtener esto es acoplar el modo 
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(bulk) correspondiente al vector de onda (k.,, k,) con el modo que se propaga en 
(k., —ky). Las soluciones de este tipo existen siempre que existan los modos bulk, 
por lo que no hay soluciones de este tipo dentro del gap. 


Buscamos soluciones a la ecuación diferencial con condiciones de borde, cuyas 
energías estén dentro del gap (bulk). Como no existen modos extendidos dentro 
del gap, las soluciones tienen que ser exponenciales de la forma exp((ik, + q)y), 
donde (k,, q) e R. Debemos imponer q < 0 a las soluciones localizadas en el borde 
del plano semi-infinito (aunque conservamos ambas soluciones para utilizarlas 
más adelante), siendo —1/q la longitud de localización. Además, para satisfacer las 
condiciones de contorno es necesario superponer los modos con +k,, y finalmente 
podemos reescribir la Ec. (3.5) para las soluciones dentro del gap como 


war) = Ce revO, my ly). (3.12) 


Utilizamos la Ec. (3.12) para obtener las soluciones al Hamiltoniano de la Ec. (3.4) 
que toman la forma 


1 +47? - ¡2 
Orly) ==) — sin (doy +04), 
Oim(y) =iy1- psin(kyy), (3.13) 


OoA(y) = +iy 1 + sin(kyy + pr), 
Ooñ(y) =+y 1 + psin(k,y), 


donde C es la constante de normalización y 


cier — kz — q+iky 0 kg + + iKky 


—|k— q+ikyl? A ke + q+ ik] ada) 


Podemos ver que la función de onda oscila y decae hacia el interior de la muestra 


A 
al E (3.15) 


donde £ = h0/eEs es la longitud característica, siendo Ej la amplitud del campo 
eléctrico Ey = (24p/c* . La parte oscilante de la función de onda k;,, es 


E? 2 l+p y? 
=al 5 - (kz - = kox / Ul==2 .1 
ky (houp y? (k q) ko 1- y A2 (3 6) 


donde ANA = hQn/y1 + 1? es el ancho del gap calculado para la matriz completa 
en la Ec. (3.4). 


La relación de dispersión <(k..) se obtiene de las raíces de la ecuación 


con un decaimiento dado por 








(1 +) - 7 (1 - 1)  nko/koy1- 12 =0, (3.17) 


3.2 Cálculos analíticos 


* Es notable ver que la longitud 
de decaimiento sólo está dada 
por los parámetros de la radia- 
ción y no por los detalles del ma- 
terial. En particular vemos que 
como consecuencia de la rela- 
ción de dispersión lineal a ba- 
jas energías, la longitud de de- 
caimiento es independiente de 
la velocidad de Fermi para gra- 
feno. 
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que toma dos valores reales para k,, dentro del gap dinámico. Las soluciones las 
denotamos con Pero enfatizando que corresponden a las soluciones para el valle 
K, —que para simplificar la notación llamamos K— y para una cinta zigzag terminada 
en sitios 4. Las soluciones toman una forma más simple cuando la energía yace 
cercana al centro del gap dinámico (2/2, donde podemos aproximar 


(+29) , 7 Ka 
2(1+12) 2(1+m2) ko 





E kg) = (kg) /dhO (3.18) 
Esta solución es lineal en k,. y corresponde a estados de borde sin masa y cuya 
velocidad es constante. La solución para e = h£2/2 cumple que k, = ko, y k» = q = 
-1/2£€ y la función de onda en el dominio temporal es 


vé _(r,t) = PEU (2EL, y 72 


o es Koy + 21 sin ji A hi o] , (3.19) 
isin koy — sin koy 


donde el superíndice indica el tipo de borde, los subíndices el valle K y el signo 
de q y donde L.,, es la longitud de la muestra asumiendo condiciones periódicas de 
contorno. Las oscilaciones perpendiculares al borde dependen de ky = Q/vp, por 
lo que no dependen de la intensidad de la radiación pero sí de la frecuencia y la 
velocidad de Fermi para grafeno. El período de las oscilaciones transversales por 
la longitud de decaimiento es 2ky£ = 1/7 y nos da un estimativo de la cantidad de 
oscilaciones que tendrán los estados de borde, por lo que para nuestra elección de 
parámetros (acople débil y << 1) se observan muchas oscilaciones antes de que el 
estado decaiga hacia el bulk. 


Cintas de grafeno 


En el caso de una cinta de grafeno zigzag, también debemos obtener las solu- 
ciones para el borde opuesto en y = W. Si suponemos que la cinta es muy ancha 
(W > £), entonces las soluciones localizadas cada borde son independientes entre 
ellas. Podemos utilizar la solución deducida anteriormente (Ec. (3.13)) para describir 
el borde en y = 0, mientras que para el borde y = W definimos y = y —- W y segui- 
mos las mismas deducciones, esta vez recordando que la solución con q > 0 es 
válida (aunque conservamos ambas soluciones), ya que queremos un decaimiento 
exponencial de la forma exp(44) para y < 0. 


am (er) = Ce a (a), (3.20) 


3. Estados topológicos de Floquet 


con 


O14(9) = Y1 - psin(kyS), 
O1B(5) = -iy1 - psin (kyG -0x), 





Ooa() = +iy/1 + usin(kyD), (3.21) 
ZA y 1+4n- y? . 
OoB(9) =*+ E sin(kyZ — pr), 
y longitud de decaimiento 
1 l+p 
qe : 3,29 
¡Ta (3.22) 


Comparamos el decaimiento exponencial q con q y notamos que a una misma 
energía difieren entre ellos y se cumple 


y le 
e po 





(3.23) 





l-u p? 
ka =k 3.24 
WVil+p ¡AR il 
La relación de dispersión se obtiene de las soluciones de 
(1 =p) (+ po) Fky ko 1 12 =0, (3.25) 


que denotamos ¿4 , para las funciones de onda definidas en el valle K' para un 
borde zigzag terminado en sitios B. La relación de dispersión cerca del centro del 
gap dinámico toma la forma 


1 hOn kz 
2(1 +1?) “A 


siendo la dependencia de la energía lineal en k.,, por lo que describe excitaciones 


Elk) = (ko) /AO e AO 





(3.26) 


electrónicas sin masa y dispersivas, localizadas en el borde. 


En la Fig. 3.4 vemos que la relación de dispersión analítica está en muy buen 
acuerdo con los resultados numéricos para el modelo ee ado La solución 
exacta en el dominio temporal para e = h02/2 cumple con k., = -G = -1/2£, ky = ko y 
toma la forma 

1 


B _ ¡miy/2€ 3/2€ 
r,t)=e e E 
DE sl ) DEL, 


isinkoZ Y ; sin kot 
s ( 0 Jen + 0 
cos kyT - 


-)| - (8.27) 


cos kyí — 29 sin kyZ 


3.2 Cálculos analíticos 














Fig. 3.4: Comparación entre las 
relaciones de dispersión utilizan- 
do los resultados analíticos (círcu- 
los azules) y los numéricos (pun- 
tos en escala de color) utilizando 
un modelo tight-binding para una 
cinta de N = 1000, AQ = 0,3t y 
z = 1073. La escala de color indi- 
ca el peso sobre la réplica n = 0, 
mientras que los cálculos están 
realizados con 4 réplicas de Flo- 
quet. Figura adaptada de [53]. 
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Fig. 3.5: Relación de dispersión 
(a) para los estados de borde en el 
valle K, y (b) para los estados en 
el valle K_. La línea continua roja 
corresponde a la dispersión los es- 
tados localizados en el borde y = 0 
(terminado en sitios 4) y la línea 
de puntos azules a la dispersión 
de los estados localizados en el 
borde opuesto y = W (terminado 
en sitios A) de una cinta zigzag 
muy ancha (W > £). Para un mis- 
mo borde, los estados en los valles 
K, y K_ se propagan en la mis- 
ma dirección, siendo la velocidad 
opuesta en el otro borde [53]. 
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En la Fig. 3.5 (a) se grafican las dispersiones que son soluciones de las 
Ecs. (3.17) y (3.25) (más allá de la aproximación a e = h(2/2). Vemos que las 
velocidades de los estados que pertenecen a un borde 0 (1) y al otro Pr. (1) 


(calculadas como 





o bien como el valor medio de ((v)) = vp((0.,))), son 


opuestas y cumplen 


A B _ 3 
UK 4 2 Tp 


(3.28) 
para energías cerca del centro del gap dinámico. 

0.65 

0.6 | 

0.55 


ly 0.5 


0.45 





ky/Ko 


ky/Ko 


Las soluciones para el valle K_ pueden ser fácilmente obtenidas si realizamos 
un reordenamiento de la base tal que el Hamiltoniano tenga la misma forma para 
ambos valles. En la base reordenada, las funciones de Floquet son ¿x_(r) = 
([-u18(r),04(r)], [-208(r), 204(r)])*, que equivale a aplicar el operador ¡0,,7p 
al vector que define las funciones de onda en el valle K,. Por lo tanto, 


dx (1) = -i0yTbk, (1) 


3.29 
Pr (1) E -ioyTodj, ¿(1), ' E 


son las soluciones en K_ para bordes terminados en la familia de sitios 4 o B y 
dependen de las soluciones en K, para la familia de sitios opuesta. En particular se 
cumple que 


UR = UR, - (3.30) 


(3.31) 


por lo que la velocidad del estado de borde en cada valle tiene el mismo signo, y 
signos opuestos al cambiar de borde. En la Fig 3.5 (b) se muestra la dispersión 
para el valle K_, donde vemos que las velocidades de los estados en K, y K_ 
para un mismo borde de la muestra tienen el mismo signo y en el borde opuesto el 
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signo opuesto, lo que hace que los estados de borde inducidos por la radiación sean 
quirales. 























Finalmente, podemos comparar tanto la dispersión de los estados de borde, 
como el decaimiento y las oscilaciones obtenidas analíticamente, con los resultados 
numéricos del modelo tight-binding. Recordemos que en este modelo los hoppings 
entre sitios dependen del tiempo de la forma 


27 fi 
t¡¿(t) = texp (E / A(t) dr) Ñ (3.32) 
Do ri 
Para construir el Hamiltoniano de Floquet obtenemos las componentes de Fourier de 
estos elementos de matriz y haciendo uso de la identidad de Anger-Jacobi llegamos 
a 


st, (3.33) 


Las componentes de Fourier de los hoppings +, ¿(t) pueden interpretarse como 
los hoppings entre distintas réplicas, tal que el elemento de matriz entre un sitio de 
la réplica m y un sitio vecino de la réplica m + n están acoplados por £;'.. Puede 
probarse que estos hoppings son proporcionales a z” para intensidades bajas, donde 
2= Ae es un parámetro adimensional que controla la intensidad de la radiación. 
Para cualquier valor de la intensidad z puede demostrarse (ver Apéndice A.1) que 
son proporcionales a J,, (2), que es la función de Bessel de orden n. La fase que 





Fig. 3.6: Espectro de cuasiener- 
gías de una cinta zigzag de an- 
cho W = 213nm, irradiada con 
luz circularmente polarizada de fre- 
cuencia Ahí? = 0,3t, e intensidad 
z = 0,05, en el panel (a) para el 
valle K, y en el panel (b) para 
el valle K_. En color se muestra 
el peso sobre la réplica n = O de 
blanco (w = 0) a negro(w” = 1), 
mientras que los estados de bor- 
de se destacan en color azul (iz- 
quierda) y rojo (derecha). La distri- 
bución espacial de la densidad de 
probabilidad de los estados de bor- 
de sobre la subred A se muestra 
en (c) para K, y en (d) para K_. 
El panel (e) muestra un corte en 
el mapa espacial (c) para sitios de 
tipo A y B (cuadrados y círculos) 
mientras que los resultados analí- 
ticos se muestran en líneas a tra- 
mos y continuas. Figura adaptada 
de [34]. 
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adquiere el hopping entre réplicas está dada por 





tij= toa Mz) pij= tan 


(rj=r09 (3.34) 


T 
var. Pida — [rj— ri) 


donde «p;,¿ un factor geométrico que da información de la red y toma esa forma para 
polarización circular. 


Con el Hamiltoniano de Floquet en forma matricial del modelo tight-binding 
podemos resolver espacialmente a lo largo del ancho de la cinta la densidad espacial 
de los estados de borde. En la Fig. 3.6 se muestran los estados de borde para el valle 
K, y K-_ para una cinta zigzag de ancho W = 213nm, equivalente a N = 2000. En 
los paneles (a) y (b) vemos la relación de dispersión donde se distinguen las bandas 
del grafeno irradiado, en escala de color se indica el peso sobre la réplica n = 0, con 
puntos negros para w” = 1 y blancos para 1w/ = 0. El peso de los estados de borde es 
cercano a w/ = 0,5 y se distinguen con los mismos colores (azul y rojo) utilizado para 
mostrar la densidad espacial de los estados en (c) y (d) sólo se muestra la densidad 
de probabilidad de los sitios de tipo A para destacar las oscilaciones. El estado en 
color azul está localizado a izquierda de la muestra, mientras que el estado de color 
rojo a derecha, tanto para el valle K, como el K_. Las pequeñas líneas horizontales 
claras de los paneles (c) y (d) de la densidad espacial se corresponden a cruces 
evitados entre los estados de borde y las bandas de réplicas de Floquet de orden 
superior n = 2,—1. En la Fig. 3.6 (e) se muestra un corte en el mapa espacial de 
estados en (c) para e = h02/2 pero mostrando las densidades para sitios de tipo A y 
B. Podemos comparar el decaimiento y las oscilaciones de las funciones de onda 
localizadas en los bordes. Las líneas indican las soluciones analíticas: para los sitios 
de tipo A líneas oscuras continuas y B líneas claras a trazos. Los puntos indican 
el peso sobre la réplica n = O de las funciones de onda obtenidas con el modelo 
tight-binding resuelto espacialmente, para sitios de tipo A círculos rellenos y sitios 
de tipo B cuadrados vacíos. Tanto el decaimiento como las oscilaciones muestran 
un perfecto acuerdo para energías cercanas a h02/2 (< = 0,15 t en la figura). 


Resultados numéricos para una muestra semi-infinita 


Una forma de eliminar los efectos de tamaño finito y corroborar la quiralidad de 
los estados de borde es tomar la densidad local de estados (LDOS) resuelta en el 
momento transversal k., cerca del borde de un plano semi-infinito de grafeno 


1 
pijle,k) = -—Im(Gjy(e +10*,k)), (3.35) 


y sumar sobre los sitios j entre el borde y un ancho r 


PedgekEs k) a y) pjjle, k), (3.36) 
j=1 
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tal que satistaga r > £/ap, de modo que el peso de los estados de borde esté com- 
pletamente incluido. Para simular un plano semi-infinito, tomamos el Hamiltoniano 
de Floquet-Bloch de una cinta de ancho N y aplicamos el proceso de decimación 
hacia el interior de la cinta. El algoritmo de decimación aplicado a grafeno está 
explicado en detalle en Calvo y col. (2013) [43] y en general se utiliza para obtener 
las funciones de Green o espaciales. En el caso del plano semi-infinito, simulamos 
una cinta muy ancha e infinita en una dirección. La invarianza traslacional en esta 
dirección nos permite definir el momento k. El Hamiltoniano de Bloch para una 
cinta de grafeno toma una estructura tridiagonal por bloques a la que podemos 
aplicarle el mismo algoritmo de decimación. La versión rápida de este algoritmo 
nos permite simular una cinta de ancho 2Y en M iteraciones, a un muy bajo costo 
computacional. Según la resolución de energía buscada, es suficiente con tomar 
entre M = 20 y M = 40 iteraciones para que la LDOS cerca de un borde no se vea 
influenciada por los estados del borde opuesto (notar que la cinta tendrá anchos 
entre N » 10% y N » 101?) y obtenemos efectivamente la LDOS sobre el borde de 
un plano semi-infinito. 


La densidad de estados local ahora corresponde a un solo borde y podemos 
distinguir la quiralidad de los estados de borde ya que la pendiente en la relación 
de dispersión es su velocidad (v = 2 
LDOS cerca del borde de un semiplano de grafeno irradiado con luz de polarización 


). En la Fig. 3.7 vemos un mapa de la 


circular y A(0 = 3t, z = 0,5, lo que nos da una longitud característica de £ = 2ag 
(tomamos r = 50). Estos parámetros están exagerados para denotar efectivamente 
que mientras que el en gap en e = 0 hay un estado con velocidad positiva, en el 
gap dinámico e = h(2/2 vemos dos estados, uno por valle, ambos con velocidad 
negativa. En la Fig 3.7 es evidente que las quiralidades de los estados en uno y otro 
gap son opuestas. La escala de color es meramente ilustrativa ya que al tomar una 
densidad local de estados cerca del borde, obtenemos un valor parcial y las bandas 
de grafeno (zonas oscuras) y los estados de borde no pueden estar normalizados al 
mismo tiempo. 





Fig. 3.7: Densidad local de esta- 
dos cerca del borde de una mues- 
tra semi-infinita de grafeno irradia- 
da con luz de frecuencia A£) = 3t, 
intensidad z = 0,5 y circularmente 
polarizada. Se toman los primeros 
50 sitios del borde para obtener el 
mapa de color del peso de los es- 
tados sobre la réplica n = 0, mien- 
tras que se consideran 5 réplicas 
de Floquet en el cálculo. Los es- 
tados de borde cruzan los gaps 
ene =0 y e = +h0/2 con quirali- 
dades opuestas. Figura adaptada 
de [53]. 
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Fig. 3.8: Densidad local de esta- 
dos resuelta en k.. para paráme- 


tros más realistas de radiación. En 
el panel (a) se puede observar que 
el gap en e = 0 está fuertemente 


modificado por las réplicas de or- 
den superior y es mucho más chi- 


co que el gap A; en e = ñ0/2. En 
el panel (b) el gap en e = A0/2 
es más robusto y está cruzado por 
un estado de borde (en círculos 


azules se muestra la solución ana- 


lítica). Aumentando la precisión y 
la escala de la LDOS vemos en 
(c) que el gap de orden superior 


A presenta 3+1 estados quirales. 


Figura adaptada de [53]. 
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3.3 Topología de las bandas de Floquet 


En la Fig. 3.8, utilizamos parámetros más realistas de radiación (h0 = 0,05t 
(135meV), z = 2,8 x 107?) para simular un láser en el infrarrojo medio, e incluimos un 
16 de réplicas de Floquet, aunque los resultados para no difieren de considerar solo 6 
réplicas. En el panel (a) observamos tanto el gap dinámico A, como el gap en e =0, 
aunque para estos parámetros este último es mucho más pequeño y la hibridización 
de los estados con las réplicas +1 difumina la dispersión cerca de k., = 0. El gap 
dinámico es más robusto a los efectos de las réplicas de orden superior, por ejemplo 
en el panel (b) podemos ver estados de borde, en excelente acuerdo con la solución 
analítica (círculos azules). En el panel (c) las réplicas de orden superior también 
tienen un gap As y asociado a este gap más pequeño (dos órdenes más grande 
en 7), nuevos estados quirales aparecen con la misma quiralidad. La contribución 
de estos estados a la LDOS promediada en el tiempo es despreciable, aunque el 
hecho de que estos estados existan en el espectro de Floquet tiene importantes 
consecuencias sobre la topología de las bandas (ver Sección 3.3) y es tratado en 
detalle en el Capítulo 5. 
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3.3. Topología de las bandas de Floquet 


En esta sección dilucidamos si nos encontramos frente a estados de borde 
accidentales, propios de la terminación zigzag, o bien si los estados de borde son 
independientes de la terminación y ocurren para cualquier tipo borde. De manera 
similar a los aislantes topológicos estáticos, existen invariantes topológicos para los 
sistemas periódicos en el tiempo que nos permiten anticipar la existencia de estados 
de borde. Estos estados topológicos de Floquet que son dispersivos, quirales y 
robustos frente a defectos. 


3. Estados topológicos de Floquet 3.3 Topología de las bandas de Floquet 


Las propiedades topológicas de los Hamiltonianos estáticos de sistemas de dos 
y tres dimensiones es un tema maduro en la física actual [23, 57, 58, 59, 60]. En 
dos dimensiones, el número de Chern, definido como una integral de la curvatura 


de Berry sobre la primera zona de Brillouin caracteriza la topología de las bandas  tcabe aclarar que existen más 


presentes en el sistema Y . En un sistema de dos bandas separadas por un gap, el de un invariante topológico para 
número de Chern C; de la banda de menor energía caracteriza las propiedades  Caracterizar distintos aspectos 
de un sistema físico. Además, 
la relevancia de la información 
que nos brinda el invariante to- 
al número total de estados de borde n, donde n es la suma con signo —dado por la pológico está relacionada con 


quiralidad— de los estados de borde y coincide con la cuantización de la conductancia las simetrías conservadas por 
el Hamiltoniano del sistema [60] 


topológicas de los estados de borde que ocurren dentro del gap. Como consecuencia 
de la correspondencia entre volumen y superficie, el número de Chern es equivalente 


del sistema [21]. 


En el caso de Hamiltonianos periódicos en el tiempo debemos hacer uso de un 
invariante topológico distinto [61, 36, 62, 63] debido principalmente a que el espectro 
de cuasienergías del Hamiltoniano de Floquet es periódico y no acotado. Este inva- 
riante topológico llamado número de winding, deducido en Rudner y col.(2013) [36], 
describe correctamente la topología de un sistema bidimensional de Floquet sin 
pedir ninguna simetría extra conservada. Una discusión más detallada sobre los 
invariantes en sistemas forzados periódicamente se dará en el Capítulo 5. Por el 
momento es suficiente con comprender que, a pesar de que el espectro de Floquet 
es no acotado, las propiedades topológicas y el espectro de cuasienergías efectivo 
de un Hamiltoniano de Floquet puede ser bien descripto utilizando un Hamiltoniano 
de Floquet truncado a un número M de réplicas, con el cuidado de elegir un M sufi- 
cientemente grande para que todos los procesos fotónicos de interés estén incluidos. 
En el régimen de acople débil, es decir de intensidades bajas, el número de réplicas 
necesarias para dar cuenta de los efectos relevantes es M > A/AQ, donde A es el 
ancho de banda del sistema. Podemos construir entonces, un invariante topológico 
que de cuenta del número y quiralidad de los estados topológicos que se encuentren 
entre las bandas de Floquet n y la banda (n + 1). El invariante topológico es W (e), 
el número de winding, definido para energías dentro del gap como la suma de los 
números de Chern de las bandas inferiores, 


W(s) = Y G,, (8.37) 
k=1 


donde «e yace en el gap entre las bandas de Floquet n y (n +1). 


De las M bandas de Floquet que se solapan dentro del rango de energías de 
del ancho de banda del sistema original (réplica n = 0), sólo algunas producen 
modificaciones significativas en la DOS media. En particular, si el acople entre 
réplicas es de orden y << 1, la réplica de Floquet de orden n tendrá peso sobre la 
DOS media de orden 7”. Vemos entonces que las correcciones más significativas 
pueden ser consideradas incluyendo solamente las réplicas n = -1,0,+1. Las 
réplicas n = —1,+1 tienen el efecto de renormalizar los puntos de Dirac abriendo un 
gap de orden 1? a energías ey = 0, mientras que los cruces entre las réplicas n = +1 
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abren un gap de orden 7 a energías e. 1/2 = +1/2A0. Existen además solapamientos 
de réplicas de orden superior cuyas implicancias sobre la topología del sistema 
será consideradas en el Capítulo 5, pero que no relevantes para mediciones de 
cantidades medias, como la DOS media. 


A bajas frecuencias (ñRh( << t), tanto ey como e, 1/2 yacen en el régimen de bajas 
energías y podemos utilizar el modelo k - p para obtener los Hamiltonianos efectivos 
que describen al sistema a energías próximas a ep y €,1/2. Bajo la aproximación 
de acople débil (7 << 1) las bandas están cuasi- degeneradas para las energías 
de los gaps y para regiones del espacio k donde las bandas se solaparían en 
ausencia de radiación (7 = 0). En las zonas de cuasi-degeneración se encuentra 
la mayor contribución al número de Chern de una banda, ya que son las zonas 
de mayor curvatura de Berry. Las zonas del espacio k cercanas a k = O (marcado 
con rectángulos en la Fig. 3.2) son las zonas de donde provienen las principales 
contribuciones a W (<p), el número de winding del gap a energía cero, mientras que 
las contribuciones a W (s1/2) provienen de k = 0/2vp (marcado con círculos en la 
Fig. 3.2). 


Para calcular estas contribuciones al número de winding primero debemos 
identificar las zonas del espacio k donde se abre un gap en e, luego obtener el 
Hamiltoniano efectivo en esos puntos para cada valle, que es una matriz 2 x 2 de la 
forma 


HUXMk, v)=he¿,(k)-0 


El paso siguiente es calcular la contribución al número de Chern de cada cruce y 
sumar las contribuciones de ambos valles, suponiendo que la curvatura de Berry 
decae lo suficientemente rápido fuera de estas zonas. 


El Hamiltoniano efectivo que describe el gap en ey está dado por 
heov(k) = hop(kg, Vky, UN Q/Up) , 


donde y = +1 es el índice de valle. La contribución al número de Chern de cada 
valle es 1/2 y dado que las contribuciones tienen el mismo signo, obtenemos un 
total de 1 estado. Para analizar el gap en e, podemos considerar solamente 
las réplicas n = 0,+1 y expresar el Hamiltoniano de Floquet en la base de las 
autofunciones de ?y, el Hamiltoniano no perturbado. En esta base, el cruce entre 
la rama como energía positiva (+hvpk) de la réplica n = 0 y la rama con energía 
negativa (h0 — hupk) de la réplica n = 1 puede describirse con el Hamiltoniano 
efectivo 


he; 2, (k) = hup((k-— ko) k — ko nÉr), 
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donde és = sin 0:01 + cos 09,2. La contribución de este cruce al número de Chern es 
de —1 por cada valle, lo que nos da un total de -2 estados topológicos. Notoriamente 
el gap en ey, tiene el doble de estados que el gap en ep y con quiralidad opuesta. 

Las propiedades topológicas que deducimos a bajas energías prueban ser 
válidas incluso para algunos casos de altas frecuencias (ver Fig. 3.7), donde e, ¡2 
yace más allá de la singularidad de van Hove para grafeno. 

Cabe aclarar que sólo consideramos los cruces que generen estados topológicos 
con peso significativo sobre la réplica n = 0. Podemos anticipar, por la periodicidad 
del espectro de Floquet y por la simetría electrón-hueco del grafeno, que el peso de 
los estados topológicos exactamente en el centro del gap, será igual al promedio 
de los pesos de las réplicas de Floquet. Por ejemplo, el gap en ey es originado 
por una renormalización de orden 7? de las bandas de la réplica n = 0, o bien 
podríamos pensarlo como un cruce de dos bandas pertenecientes a la réplica n = 0, 
en cualquier caso podemos anticipar que el peso espectral del estado topológico 
sobre la réplica cero es w% = 1- O(7?). En cambio, en el gap en £1/2 el cruce es 
entre la réplica n = 0 y la réplica n = 1 y las simetrías del espectro nos aseguran que 
exactamente a energía h(2/2 el peso de los estados topológicos generados debe 
ser w% = 1/2- O(m). 

En la Fig. 3.1 vimos cómo los estados en e; ¡2 tienen peso 1w% = 1/2 (en color rojo), 
mientras que los estados topológicos que ocurren en ey no llegan a desarrollarse ya 
que se encuentran fuertemente afectados por los efectos de tamaño finito (el gap es 
más pequeño) y por los acoples con las réplicas de orden superior. 


3.3.1. Distintos tipos de borde 


En la sección anterior vimos que el número y la quiralidad del total de estados 
de borde está caracterizado por un invariante topológico: el número de winding. 
Al igual que el número de Chern en aislantes topológicos en equilibrio, la corres- 
pondencia volumen-superficie [21] asegura que los estados de borde aparecen 
independientemente de la terminación de la muestra de grafeno. 

Los estados topológicos de Floquet que emergen dentro del gap dinámico, 
además de ser robustos frente al desorden estructural, también deben aparecer en 
el borde de una muestra de grafeno irradiada, para cualquier dirección del borde 
dentro de la red de grafeno y sin importar los detalles de la terminación del borde. 
En la Fig. 3.9 se muestran algunos casos típicos de bordes que aparecen en las 
muestras de grafeno. En el panel (a) se muestra una cinta con bordes zigzag como 
referencia, en el panel (b) una cinta armchair, en el panel (c) una cinta zigzag con 
bordes de tipo Klein [64] y en el panel (d) una cinta zigzag con bordes calados (cove 
edges) [65]. 

Los distintos tipos de borde cambian la proyección de la relación de dispersión, 
es decir cómo se pliega la zona de Brillouin y también puede cambiar los detalles de 
la dispersión de los estados topológicos. Sin embargo su número total y quiralidad 
permanecen invariantes. En los paneles (a), (c) y (d) se pueden distinguir los valles 
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Fig. 3.9: Espectro de cuasiener- 


gías para cintas de grafeno con 
distintas terminaciones, irradiadas 
con luz circularmente polarizada 


de frecuencia Ahí? = 0,8t e intensi- 
dad z = 0,04. En color se mues- 


tra la contribución a la densidad 


de estados promediada en el tiem- 


po que es la proyección sobre la 
réplica n = O (en los cálculos se 
incluyen desde la réplica n = -2 
alan = +2. En los paneles (a), 
(c) y (d) se muestran los resulta- 
dos para cintas zigzag de ancho 
W = 31,95mm con bordes zig- 
zag, bordes de tipo Klein y bor- 
des calados, respectivamente; y 
en el panel (b) los resultados pa- 
ra una cinta armchair de ancho 
W = 31,48mm. Figura adaptada 
de [53]. 


+Los detalles de los cables no 
son relevantes en el límite de 
acople débil, aunque sí es im- 
portante que los cables no es- 
tén irradiados para que las ocu- 
paciones electrónicas en los ca- 
bles estén bien definidas[66] y 
los procesos de disipación ocu- 
rran en ellos. 
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K, y K_, que alojan un estado de borde cada uno; mientras que para una cinta 
armchair (b) ambos conos están plegados sobre el mismo punto (k., = 0) y los 
estados de borde se ubican cerca del borde de banda 

Podemos verificar que en todos los casos dos estados cruzan el gap dinámico, 
un estado por valle K,. El gap en e = 0 es mucho más pequeño para los parámetros 
elegidos y esto, junto a los efectos de tamaño finito de las cintas de grafeno, hacen 
que no se desarrolle el estado de borde. Vemos en la Fig. 3.9 que los detalles del 
borde sí cambian la relación de dispersión, pero no el número efectivo de canales de 
conducción. El número de estados de borde dispersivos junto a la quiralidad están 
dados por el invariante topológico. 


3.3.2. Robustez del transporte frente a defectos 


Para mostrar numéricamente la quiralidad y la robustez del transporte frente a 
los defectos en la red podemos calcular la probabilidad total de transmisión entre 
dos sitios (1,52) de la muestra irradiada, separados por una distancia d = 63nm a 
lo largo de la dirección x, que equivale a una distancia de 512 celdas unidad de una 
cinta zigzag. La probabilidad de transmisión está dada entre dos cables de prueba 
virtuales, débilmente acoplados a los sitios s; y s2 * y para verificar la quiralidad 
de los estados, fijamos el sitio s; en la capa £L a izquierda, cerca del borde de la 
muestra, mientras que variamos la ubicación transversal (coordenada y) del sitio sa 
en la capa R (a derecha). 

En la figura 3.10 se muestra la probabilidad de transmisión T;,-.s, en azul y 
Ts¡>s,. La marcada asimetría de la probabilidad de transmisión verifica la quiralidad 
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de los estados. En los paneles (a) y (b) la muestra de grafeno está libre de defectos 
y podemos observar cómo la probabilidad de transmisión decae exponencialmente 
a una tasa 1/1) (líneas negras a trazos), siguiendo el comportamiento de la LDOS 
media de los estados de borde. Mientras que en el panel (a) el estado localizado 
cerca del borde y = 0 ay se propaga desde la capa R hacia la L (puntos rojos), en el 
panel (b) el estado localizado cerca del borde y = 750 ag tiene velocidad opuesta y 
se propaga desde la capa L hacia la R (puntos azules). 


El desorden en un material conductor, tiende a disminuir la conductancia y a 
producir efectos de restrodispersión. En cambio en estados conductores quirales, 
los estados con direcciones de propagación opuestas se encuentran separados 
espacialmente, lo que hace que sean robustos frente a muchos tipos de defectos. En 
los paneles (c) y (d) de la Fig. 3.10, verificamos que los estados son robustos frente 
al desorden estructural. Entre las capas L y R introducimos un 0,1% de vacancias 





Fig. 3.10: Distribución espacial de 
la probabilidad de transmisión lo- 
cal entre dos sitios (s1,s2) Ubica- 
dos en las capas L y R de una 
cinta de grafeno irradiado con luz 
de polarización circular, frecuencia 
hQ = 2,7eV y acople y = 3/16, a 
un nivel de Fermi cercano al gap 
dinámico e = 0,497A%2. Ambas ca- 
pas están separadas por 63nm, co- 
mo se denota en líneas de puntos 
en los paneles. El sitio s, se ubi- 
ca en el segundo átomo desde la 
izquierda (paneles (a), (c) y (e)); 
y el sitio sz en el segundo átomo 
desde la derecha (paneles (b), (d) 
y (f)). La coordenada y entre es la 
posición del sitio s2 en la capa R y 
recorre el ancho de la cinta 750 ao, 
siendo ap la distancia interatómi- 
ca. Se muestran las probabilida- 
des de transmisión T',, ,s, en azul 
y Ts,>s, en rojo, normalizadas a 
su valor máximo. Los paneles (a) 
y (b) corresponden a una muestra 
sin defectos, los paneles (c) y (d) 
muestran los resultados cuando se 
introduce un 0,1% de vacancias 
entre las capas L y R y en los pa- 
neles (e) y (f) se cortó un pequeño 
trozo de 2 nm de ancho y 15,7 nm 
de largo del borde derecho. Figura 
adaptada de [34]. 
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que, más allá de que inducen dispersión de un borde a otro, los estados inducidos 
por radiación mantienen sus propiedades quirales de transporte en un alto grado. 
Otro tipo de defecto se muestra en los paneles (e) y (f), donde una franja de 2 nm 
de ancho y 15,7 nm de largo fue cortada del borde derecho de la muestra entre las 
capas L y R. En el panel (e), el perfil de transmisión local es prácticamente idéntico 
a la cinta sin defectos del panel (a), por lo que el estado en el borde y = Oag no se 
ve afectado por las modificaciones que se hagan en el borde opuesto. En el panel 
(f), el perfil de transmisión se ve levemente afectado pero mantiene sus propiedades 
quirales de transporte. Esto demuestra que los estados de borde se comportan 
como canales unidireccionales de transporte, quirales y robustos frente a defectos. 


3.4. Contribuciones originales 


= Desarrollo de códigos propios para obtener resultados espectrales y resolución 
espacial en el espacio real de los aislantes topológicos de Floquet, tanto para 
cintas de grafeno como para muestras semi-infinitas. 


= Describimos analíticamente la aparición de estados de borde de Floquet en una 
muestra de grafeno irradiada de terminación zigzag. Junto con la forma explícita 
de la función de onda para cada subred 4 y B y para cada valle K, y K_ pudimos 
deducir la longitud de localización, la velocidad y la relación de dispersión de los 
estados. Notoriamente vemos oscilaciones en la densidad de probabilidad de los 
estados de borde y que su velocidad es independiente de la velocidad de Fermi 
en grafeno. 


= Verificamos numéricamente la existencia de los estados topológicos de Floquet 
en un modelo de red, en muy buen acuerdo con los resultados analíticos. También 
este modelo de red nos permite comprobar la quiralidad de los estados de borde 
mediante un cálculo (local) de transporte. 


= Todos estos resultados están soportados por el cálculo de las propiedades 
topológicas de las bandas de Floquet del grafeno irradiado. Verificamos que los 
estados de borde en cintas de distintos tipos de terminaciones coinciden en 
número y quiralidad con el número de winding. 
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AISLANTES TOPOLÓGICOS DE FLOQUET 








Contenidos 
4.1. Transporte a dos terminales ..................«..... 51 
4.1.1. Mismatch de los modos electrónicos. ............ 53 
4.2. Transporte multiterminal .............. o... o... +... 54 
4.2.1. Plateaus de conductancia Hall ................. 57 
4.3. Contribuciones originales .............. 0... o... . 58 
Abstract 


La característica más distintiva de los aislantes topológicos es la cuantización de los 
Pplateaus Hall de conductancia. Esta cuantización está dada por el número y quiralidad 
de los estados de borde cuya presencia se deduce de las propiedades topológicas del 
bulk. En cambio, para los aislantes topológicos de Floquet, los plateaus de conductancia 
no están cuantizados en exactamente el total de estados de borde. A diferencia del caso 
estático, sólo los estados topológicos que suman al promedio temporal de la densidad de 
estados son los que contribuyen a la magnitud de la conductancia. Utilizando un esquema de 
scattering con una muestra de grafeno irradiada conectada a cables no irradiados mostramos 
que la cuantización de la conductancia Hall no está directamente relacionada al invariante 


topológico de las bandas de Floquet. 


Los aislantes topológicos de Floquet son sistemas en los que las propiedades 
topológicas son inducidas por un potencial periódico en el tiempo. La luz circular- 
mente polarizada irradiada sobre grafeno abre gaps en el punto de Dirac y el gap 
dinámico (ver Sección 3.3). Estos gaps son no triviales y alojan estados topoló- 
gicos de borde de no equilibrio, lo que los diferencia de los estados topológicos 
que ocurren en los aislantes topológicos en sistemas estáticos y en equilibrio. En 
los sistemas de Floquet, los estados son de no equilibrio lo cual pone dificultades 
en determinar la ocupación de los estados. Los argumentos termodinámicos que 
definen la función de ocupación f(<) de los estados en sistemas de equilibrio no 
puede aplicarse al formalismo de estados de Floquet directamente [67], ya que 
el espectro de Floquet es no acotado. Si consideramos un sistema cerrado con 
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Fig. 4.1: Esquema utilizado en 
los cálculos de transporte. (a) Un 
láser ilumina el centro de una re- 
gión de scattering conectada a ca- 
bles no irradiados. (b) La configura- 
ción multiterminal en forma de «H» 
permite corrientes de bombeo ya 
que los cables son no equivalentes. 
Las terminales 1 y 6 tienen bor- 
des de tipo zigzag, mientras que 
las terminales 2 a 5 tienen bordes 
armchair. (c) Cuando los cables 
son equivalentes entre sí el bom- 
beo de corriente es nulo. Figura 
adaptada de [80]. 
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un potencial dependiente periódicamente del tiempo, la temperatura del sistema 
puede tender a infinito [68], aunque puede probarse que el sistema se estabiliza 
si tenemos en cuenta la localización de estados de muchos cuerpos [69]. Uno de 
los caminos para determinar la ocupación de los estados de Floquet es conectar 
el sistema a reservorios térmicos [70] relevantes, que pueden ser: electrón-fonón, 
recombinación radiativa, emisión de fotones, scattering asistido por fotones, etc. 
Según qué mecanismos sean tenidos en cuenta, los baños con reservorios térmicos, 
y la configuración de acople entre el sistema forzado y el baño electrónico [71, 72], 
será el estado estacionario al que llega el sistema y este estado define las propie- 
dades de transporte del sistema. Otra punto de vista a considerar es cuáles son 
los estados de Floquet accesibles según el estado inicial del sistema y el protocolo 
de encendido del potencial periódico aunque en presencia de no-integrabilidad e 
interacciones el estado estacionario a tiempos largos es independiente del estado 
inicial [73]. Estas cuestiones siguen siendo tema activo de investigación y algunas 
propuestas pueden encontrarse en las Refs. [38, 74, 75, 76, 77, 78, 79]. 


(a) 


área 


¡iluminada 





corrientes no nulas 


(a cero voltaje) no hay corriente de bombeo 


Una situación física donde donde podemos definir las ocupaciones de un aislante 
topológico de Floquet, como el grafeno irradiado, es cuando la inyección y extracción 
de carga se hace mediante cables no irradiados y la región irradiada o de scattering 
es lo suficientemente pequeña para que no ocurran mecanismos de disipación al 
interior de la región. En este esquema, podemos usar el formalismo de scattering 
de Floquet, que ha sido claramente mostrado como equivalente al formalismo 
de Keldysh [52] para obtener las propiedades de transporte. Las ocupaciones 
están definidas en los cables y también allí ocurren los procesos de disipación [81, 


4. Transporte en aislantes topológicos de Floquet 4.1 Transporte a dos terminales 


82, 66]. En la Fig. 4.1 (a) se esquematiza una muestra de grafeno irradiado que 
está conectada a cables no irradiados en una configuración multiterminal. En una 
configuración multiterminal la corriente media sobre la terminal a, es el promedio 
temporal Z., = 2 fo Za(t)dt y toma la forma [81]: 


Te Y SAL TARO te], (4.1) 


Bra n 


donde To es la probabilidad de transmisión del cable a al P para un electrón 
que absorba n fotones (emita |[n| fotones si n < 0) y fa(<) es la función de Fermi en 
el cable a. 

La conductancia diferencial se obtiene al expandir linealmente la Ec. (4.1) en 
los voltajes de bias, aunque al estar rota la simetría de reversión temporal, no 
podemos asumir que las probabilidades de transmisión 7. ¿ y Ts, Sean iguales. 
Esta desigualdad se extiende a la suma 2 ¿Tg a(e) + Eg Tas (e) lo que significa que 
en ausencia de voltaje existe una corriente efectiva de bombeo entre las terminales 
como se muestra en el esquema de la Fig. 4.1 (b). El término de bombeo debemos 
considerarlo en los cálculos de conductancia a multiterminal, aunque para simplificar 
los cálculos utilizamos una configuración hexagonal (Fig. 4.1 (c)) para que por 
simetría el bombeo de carga se anule. 

En este Capítulo mostraremos los resultados presentados en «Multiterminal Con- 
ductance of a Floquet Topological Insulator» (Foa Torres, Perez-Piskunow, Balseiro 
y Usaj [80]), para el cual el autor de ésta tesis ha contribuido con los resultados preli- 
minares de conductancia a dos terminales, en la explicación del empalme imperfecto 
entre los modos electrónicos de los cables y de la muestra y en discusiones sobre 
todo el trabajo. 


4.1. Cálculos de transporte a dos terminales 

En la configuración de sólo dos terminales L y R, exigimos no acumulación 
de carga (o unitariedad) por lo que 7, = -Ix = 1 y reescribimos la Eg. 4.1 en 
términos de el promedio T(<) = (Tr, L(e) + Th, n(e))/2 y la diferencia 97 (e) = 
(Tr,L(e) - Tr, n(e))/2, como 

= le 
TELS 1) + IT ELO + no) de (42) 

Expandiendo a orden lineal en el voltaje Y obtenemos 


2 2 
ía (er) VE f artos, (4.3) 


donde el segundo término es la corriente de bombeo (1 p) que no depende del voltaje 
y la conductancia diferencial queda como (2e?/h)7 (e). 
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Fig. 4.2: En el panel (a) la relación 
de dispersión efectiva para una 
cinta de grafeno zigzag de ancho 


W = 125a irradiada con luz circu- 


larmente polarizada de frecuencia 


hQO = 1,5t e intensidad z = 0,15. 


En escala de grises se denota el 


peso (de blanco para w” = 0 a ne- 


gro para w% = 1) sobre la DOS 


promedio. En el panel (b) la con- 


ductancia a dos terminales para 
una muestra de grafeno irradiada 


de largo L = 420a conectada a ca- 


bles de grafeno del mismo ancho 
que en (a). Las líneas con círculos 
rellenos azules corresponden a los 
resultados para una muestra sin 
defectos y las líneas con círculos 
vacíos rojos corresponden a los 
cálculos cuando la zona irradiada 


tiene 15 defectos (vacancias) dis- 


tribuidas aleatoriamente. En este 
último caso, el rompimiento de la 


simetría de inversión induce bom- 


beo de corriente como se muestra 
en el recuadro [80]. 
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4.1 Transporte a dos terminales 


Los cálculos numéricos fueron realizados con un código construido sobre la 
base en el módulo Kwant [83], mientras que los resultados a dos terminales fueron 
comparados con los resultados de otros códigos hechos en casa en lenguaje Fortran. 


Para simular la conductancia a dos terminales utilizamos la descripción atomísti- 
ca del modelo tight-binding de una cinta zigzag de grafeno irradiado conectado a 
cables no irradiados semi-infinitos a izquierda L y a derecha R, todos del mismo 
ancho. Los resultados se muestran en la Fig. 4.2, donde en el panel (a) vemos la 
dispersión de cuasienergías en escala de color que indica la contribución de los 
estados a la DOS promedio, es decir el peso de los estados sobre la réplica n = 0, 
siendo blanco cuando la proyección es nula (1% = 0) y negro cuando la proyección 
es unidad (1% = 1). En esta figura vemos la apertura de gaps en el punto de Dirac 
(Ap) a energía e = 0 y el gap dinámico (A;¡) a energía e = h(2/2. Estos gaps se 
reflejan en el panel (b), donde la disminución de la conductancia dentro de los gaps 
hasta conductancias del orden de un cuanto de conductancia Gy = 2e?/h. 


prístino 
desordenado 





0.75 


Los estados topológicos dan las propiedades de transporte dentro de los gaps 
y podemos ver que esta conductancia no está cuantizada y muestra una fuerte 
modulación. La conductancia dentro del gap del punto de Dirac oscila entre 0,5 Go 
y 1,0Go, siendo máxima cerca de e = 0. En particular dentro del gap dinámico la 
conductancia pasa de G = 0Go hasta G = 2Gy a medida que incrementamos la 
energía de Fermi desde e» < A0/2 a ep > hQ/2. Es necesario explicar por qué 
prácticamente la conductancia se anula dentro del gap dinámico, siendo que existen 
estados conductores quirales dentro de todo el ancho A;¡. En la Subsección 4.1.1 
se mostrará que se debe a un efecto de desempalme o mismatch. 


4. Transporte en aislantes topológicos de Floquet 


4.1 Transporte a dos terminales 





En Fig. 4.2 (b) se muestran los resultados para una muestra sin defectos en 
círculos llenos azules y para una muestra con 15 vacancias distribuidas al azar en el 
interior de la muestra en círculos vacíos rojos, siendo ambos resultados cualitativa- 
mente similares. Debemos notar que si los cables y la muestra son perfectamente 
simétricos, la simetría de inversión asegura que L y R son indistinguibles por lo 
que 7 = 0 y no hay corriente de bombeo. Sin embargo, en una configuración 
realista, cualquier defecto estructural o diferencias en el empalme entre los cables 
y la muestra, romperían la simetría de inversión. En el recuadro de la Fig. 4.2 (b) 
vemos la corriente de bombeo para el sistema con defectos. 


4.1.1. Mismatch entre los modos electrónicos de los cables y de la 


muestra 


La supresión de la conductancia en el gap dinámico y su fuerte modulación 
se observa tanto para el transporte a dos terminales (Fig. 4.2 (b)), como para la 
conductancia calculada en la configuración multiterminal (ver Fig. 4.7 (a) en la 
Sección 4.2). Este efecto no se debe a la falta de estados conductores, como 
podemos observar en la relación de dispersión efectiva en la Fig. 4.2 (a). Tanto a 
energías mayores a h(2/2 = 0,75t como a energías menores existen estados quirales, 
conductores y localizados en los bordes. 


En la Fig. 4.3 se muestra la densidad local de estados a energías (a) e = 0,7t < 
hQ0/2 y (b) e = 0,8t > h(Q/2, en la configuración multiterminal, de una muestra de 
grafeno irradiado, conectado a cables no irradiados. Podemos ver que tanto en el 
panel (a) como en el (b) existen modos extendidos en los cables, mientras que en la 
muestra irradiada se observa el gap bulk y los estados de topológicos que circulan 
por los bordes. En cambio, podemos encontrar la explicación a la supresión de la 
conductancia en el empalme entre los modos propagantes disponibles en los cables 
y los modos conductores al interior de la muestra irradiada, es decir, los estados 
topológicos de borde. 


Para explicar la supresión de la conductancia, supongamos un esquema como 
el mostrado en la Fig. 4.4, donde el cable no irradiado se encuentra a izquierda 
y la muestra irradiada a derecha, por lo que en una configuración de scattering, 
estamos interesados en los modos propagantes de los cables de izquierda a derecha. 
Los modos entrantes en los cables pertenecen completamente a la réplica n = 0 
(w% = 1 ya que no son modos excitados por la radiación) y se acoplan a los modos 
topológicos dentro de la muestra. La magnitud del acople estará dado por el peso 
de los estados sobre la réplica n = O y su pseudoespín. Debido a la textura de 
pseudoespín del grafeno, podemos asegurar que los modos que se mueven hacia la 
derecha en los cables se encuentran polarizados en pseudoespín. La proyección 
de (0,) en el valle K es negativa, mientras que la proyección en el valle K” es 
positiva. Ya que la muestra y el cable están perfectamente conectados y la radiación 
es encendida suavemente a lo largo de una longitud de 30a, no se producen saltos 
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Fig. 4.3: Densidad local de es- 
tados promedio para energías (a) 
e = 0,7t y (b) e = 0,8t. La muestra 
está irradiada con luz circularmen- 
te polarizada de frecuencia hQ = 
1,5t e intensidad z = 0,15 [80]. 








electrones 
entrantes 


n=0 


Acople entre réplicas de Floquet 
(los procesos inelásticos están permitidos) 


Fig. 4.4: Esquema que represen- 
ta el formalismo de Floquet aplica- 
do a una configuración de scatteri- 
ng. A izquierda los cables no irra- 
diados (intensidad cero) impiden 
transiciones entre distintas répli- 
cas, mientras que a la derecha la 
radiación induce transiciones elec- 
trónicas. [80]. 
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Fig. 4.5: Relación de dispersión 
para una zGNR irradiada con luz 
circularmente polarizada de fre- 
cuencia há(? =, e intensidad z =. En 
escala de color se muestra la pro- 
yección sobre la réplica n = O de 
(0,) para cada estado. La escala 
está saturada a un valor máximo 
de +0,1 para resaltar el cambio en 
la dirección del pseudoespín [80]. 
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grandes en el vector de onda k., (la dirección de propagación) cuando el modo 
electrónico pasa desde el cable hacia la zona irradiada. 

A diferencia de los modos dentro del cable no irradiado, los estados conductores 
dentro de la zona irradiada siguen una textura de pseudoespín como la mostrada 
en la Fig. 4.5, donde la textura de pseudoespín ya no está ligada a la dirección de 
propagación. De hecho, el pseudoespín sigue ligado a la dirección de k proyectado 
sobre k.,, aunque ahora, el valor de k.. ya no refleja la dirección de propagación. 
A energías mayores a h(2/2 los modos que se propagan a derecha (rojos en el 
valle K y azules en el valle XK”) tienen una proyección de pseudoespín que coincide 
con la proyección sobre (0,) de los estados de borde dentro de la zona irradiada. 
En cambio, a energías menores a h(2/2, la proyección de pseudoespín los modos 
propagantes hacia la derecha del cable es opuesta a la proyección de pseudoespín 
de los estados con velocidad positiva de la muestra irradiada. Esto hace que el 
empalme entre los modos sea nulo y consecuentemente los modos de los cables 
no pueden entrar a la muestra. Esto da origen a la supresión de la conductancia 
observada dentro del gap dinámico. 





XK "GB o, 
MESS E 
0.75 
a SS E 
= INS E 
R, NS Y y 0 
w Ñ / 
0.0 
-0.1 
0 0.5 1 
k (27/a) 


4.2. Cálculos de transporte multiterminal 


Obtenemos la conductancia diferencial en un sistema multiterminal mediante 
el cálculo de las transmitancias entre dos terminales por los que pasa corriente, 
mientras que las otras terminales juegan el rol de electrodos de un voltiímetro 
ideal. Debemos asegurar que los «electrodos» que conectamos a la muestra no 
inyecten ni extraigan corriente. Los voltímetros ideales constan de electrodos que son 
reservorios de electrones que ajustan sus parámetros internos (como el potencial 
químico y) para asegurar que no haya corriente neta sobre la terminal a: —es posible 
que los voltímetros tengan que ejercer trabajo sobre el sistema para asegurar esto. 
Sin embargo, aun puede existir una corriente de bombeo en los terminales por los 
que pasa la corriente. Partiendo de esta premisa podemos reescribir la Ec. (4.1) 
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4.2 Transporte multiterminal 





para pequeñas desviaciones ÓN. = la — O 
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h 
donde la corriente de bombeo Zp,, es nula en los voltímetros. Luego aplicamos una 
pequeña diferencia de voltaje entre las terminales de los cables por donde pasa 
corriente mientras pedimos que la corriente media (o corriente dc) sea nula en los 
electrodos de los voltímetros. Las desviaciones del potencial químico du, obtenidas 
en este proceso determinan la conductancia Hall del sistema que, aunque haya 
corrientes de bombeo remanentes, no afectan la conductancia diferencial. 


En un esquema de transporte multiterminal, ya no es suficiente pedir simetría de 
inversión para asegurar que la corriente de bombeo sea nula cuando la diferencia 
de voltaje es cero. En particular en una configuración de «H», como la mostrada 
en la Fig. 4.1 (b), la simetría de inversión (respecto de los ejes de simetría de la 
muestra, es decir el eje horizontal y vertical que atraviesan el centro de la muestra), 
nos asegura que las corrientes 1; = 15, que 13 = 15 y que 13 = 14 y no significa que 
deban ser nulas. Podemos por ejemplo tener corrientes distintas de cero, entrando 





por las terminales 2 a 5 y saliendo por las terminales 1 y 6, incluso a sin diferencia de 
voltaje. En la Fig. 4.6 podemos ver las corrientes en los cables 1 en azul, 2 en rojo y 3 
en verde para una muestra irradiada conectada a cables en la configuración de forma 
de «H». La principal diferencia entre los cables conectados es que las terminales 
1 y 6 son cintas de grafeno zigzag, mientras que los cables de las terminales 2 a 5 
son cintas de tipo armchair. Esto da como resultado diferencias en los coeficientes 
de transmisión 97;(€) = Y; (T;¿(e) - Tj¡(<)), como se muestra en el recuadro de 
la Fig. 4.6 con el mismo código de color. Podemos notar que los cables 2 y 3 al ser 
similares entre ellos tienen perfiles de corriente de bombeo similares y que su suma 
es 12 + 13 = —I1, lo que asegura que no haya carga acumulada. 


Queremos asegurar que la corriente de bombeo sea cero en todas las terminales, 
para hacer esto necesitamos que todos los cables sean equivalentes entre sí. La 
configuración necesaria para que esto suceda depende de los detalles de la red, 
que para el caso de grafeno debe tener simetría de rotación en 27/3 y simetría de 
inversión. La configuración más simple que cumple con estas simetrías análogas a 
las simetrías de la red subyacente es la configuración hexagonal. 


En la Fig. 4.7 (a) vemos los resultados para la conductancia en la configuración 
hexagonal de scattering en detalle para el gap en el punto de Dirac y en el gap 
dinámico, midiendo la caída de potencial entre los terminales 1 y 6, que son los 
terminales entre los que pasa corriente. En el panel (b) podemos observar que la 
resistencia Hall R,.y (medida entre las terminales 2 y 4, cuando la corriente pasa 
entre las terminales 1 y 6) es no nula en los gaps de Floquet inducidos por la 
radiación. Vemos que su perfil es muy similar al de la Fig. 4.2 (b), en los que la 
conductancia está fuertemente modulada, llegando incluso a anularse dentro del 
gap dinámico a energías menores a A(2/2, donde la proyección de pseudoespín es 
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Fig. 4.6: Corrientes de bombeo en 
un esquema multiterminal con for- 
ma de «H». La asimetría entre los 
cables hace que los coeficientes 
de transmisión 47;(<) sean dis- 
tintos de cero (ver recuadro). La 
muestra irradiada (con los paráme- 
tros de la Fig. 4.7) es un rectán- 
gulo de ancho W = 99a y largo 
L = 512a conectada longitudinal- 
mente a cables de grafeno zigzag 
del mismo ancho (1 y 6) y trans- 
versalmente con cables de bordes 
armchair (1 a 5). [80]. 
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4.2 Transporte multiterminal 





Fig. 4.7: Configuración hexagonal 
de grafeno irradiado con luz circu- 
larmente polarizada de frecuencia 
hQO = 1,5t e intensidad z = 0,15. 
En el panel (a) se muestra la con- 
ductancia cuando la corriente pasa 
por terminales opuestas, por ejem- 
plo la 1 y 6 y la caída de poten- 
cial se mide entre las mismas. En 
el panel (b) vemos los resultados 
para la resistencia Hall medida en- 
tre las terminales 3 y 4, cuando la 
corriente pasa entre las termina- 
les 1 y 6. Los resultados para los 
cables no dopados aparecen en 
círculos azules, para dopaje de ti- 
po n en triángulos negros (energía 
de sitio corrida en 0,75t) y círculos 
rojos (energía de sito corrida en 
1,25t) [80]. 
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antiparalela a la de los modos propagantes en los cables, mientras que es paralela 
para energías mayores a hA(/2. Los círculos vacíos azules con área sombreada 
por debajo corresponden a cables no dopados, mientras que los triángulos negros 
y círculos rojos corresponden a cables altamente dopados. Al dopar los cables, la 
energía de los modos entrantes respecto del punto de Dirac (de los cables) aumenta 
considerablemente, más allá de la singularidad de van Hove donde los modos del 
grafeno ya no son de tipo Dirac y el pseudoespín no está fijado por la dirección de 
propagación o del vector de onda k. Esto evita los efectos de mismatch (explicados 
en la Sección 4.1.1) entre los cables y la muestra irradiada. 


Cuando los cables están dopados la resistencia Hall (Fig. 4.7) (b) muestra 
plateaus de cuantización aproximados a valores de R..y = 12, para el gap en 
e = 0, mientras que para el gap en e = 0/2 la resistencia Hall esté cuantizada en 


Ray = 1/22. 














A pesar de que las intensidad y frecuencias utilizadas hasta el momento son 
muy altas, reflejan las características principales de transporte de los aislantes 
topológicos de Floquet para un sistema de Dirac con dispersión lineal. Los resultados 
presentados hasta ahora son un compromiso entre frecuencia e intensidad lo más 
baja posible para obtener gaps y estados topológicos bien definidos, manteniendo 
la viabilidad de los cálculos numéricos y la validez del modelo. Simular un sistema 
multiterminal de baja frecuencia e intensidad (parámetros más realistas) con cables 
lo suficientemente anchos para que los estados topológicos se desarrollen requiere 
un muy alto costo computacional. Para mostrar la estabilidad de los plateaus Hall, 
podemos comparar la resistencia R,.y dentro del gap del punto de Dirac entre 
distintas intensidades de láser. En la Fig. 4.8 vemos la resistencia no local para un 
láser de intensidades z = 0,15 (triángulos oscuros) y z = 0,25 (línea con cuadrados y 
área sombreada). La intensidad más alta genera un rango más ancho donde Ry es 
no nula, que se acerca aun más al valor de cuantización Ry = -1L>. 
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4.2.1. Relación entre los plateaus de conductancia Hall y el número 


de canales efectivos 


En los aislantes topológicos de Floquet el número de winding es el invariante 
topológico que relaciona las propiedades del bulk con el número de estados de 
borde del sistema. A diferencia de los aislantes topológicos estáticos, el número 
de estados de borde no está directamente relacionado con la cuantización de los 
plateaus de conductancia Hall. 


Por ejemplo, la resistencia Hall del grafeno irradiado con luz circularmente 
polarizada de frecuencia h( = 1,5t e intensidad z = 0,15 o z = 0,25 es de Ry * 
-1 Go. Esta conductancia está dada por un solo estado de borde de quiralidad 
negativa. Sin embargo, podemos observar en la Fig 4.9 (a) y (b) que existe más 
de un modo propagante en el gap en e = O para ambas intensidades. En ambos 
casos existe un par de estados de borde que cruza e = O cuando k = 7/a y otro 
par de estados en k = 0,27/a. En los paneles (c) y (d), vemos la misma relación 


de dispersión de los paneles superiores, pero proyectada sobre la réplica n = 0. 


En escala de color se muestra la contribución de los estados al promedio temporal 
de la densidad de estados de blanco (w% = 0) a negro (1 = 1). Mientras que los 
estados en k = 7/a tienen peso = 1 sobre la réplica n = 0, los estados en k = 0,27/a 
tienen peso = 0. Estos estados permanecen ocultos a este cálculo de transporte, ya 
que no tener peso sobre la réplica n = 0, los modos propagantes de los cables no 
pueden acoplarse a ellos. Más aun, podemos verificar que los estados topológicos 
cerca de k = 27/a tienen quiralidad positiva, opuesta a los estados en k = 7/a, 
como puede observarse en el recuadro de la Fig 4.9 (d), donde vemos la dispersión 
proyectada sobre la réplica n = O y pesada solo en una mitad de la cinta, para 
resaltar la quiralidad de los estados topológicos. 


Aunque los estados lejos de k., = 7/a no contribuyan a la conductancia Hall, sí 
podemos ver (Fig. 4.9 (a) y (c)) que el gap efectivo, donde la resistencia Hall es 
distinta de cero, se ve reducido por zonas donde las réplicas de orden superior se 
hibridizan con modos de la réplica n = 0, generando estados extendidos que no 
contribuyen a la conductancia no local. Este tipo de modificaciones en el ancho del 
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Fig. 4.8: Resistencia Hall (medi- 
da como se muestra en el recua- 
dro) cerca del punto de Dirac pa- 
ra la configuración hexagonal de 
grafeno irradiado con AQ = 1,5t 
con cables zigzag de ancho W = 
99a dopados tipo n (energías de 
sitios desplazadas en 1,5t). Los 
triángulos son los resultados para 
un láser de intensidad z = 0,15, 
mientras que los cuadrados para 
z = 0,25. Figura adaptada de [80]. 
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Fig. 4.9: Relación de dispersión 
efectiva para una cinta zigzag de 
ancho VW = 99a, iluminada por un 
láser circularmente polarizado de 


frecuencia Ahí? = 1,5t e intensida- 
des (a) z = 0,15 y (b) z = 0,25. 
En los paneles (c) y (d) se mues- 


tran los mismos resultados pero 
destacando (en escala de color) el 


peso de los estados sobre la ré- 


plica n = 0, en blanco w% = 0 y 


en negro w” = 1. El pequeño re- 


cuadro en (d) muestra la densidad 
local de estados sobre una mitad 


de la cinta para resaltar la distin- 


ta quiralidad de los estados cerca 
de e = 0: el estado que pasa por 
k = Tr/a tiene quiralidad opuesta 
al estado en k = 27,/a. [80]. 
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4.3 Contribuciones originales 





0 





gap no se observan en el gap dinámico a e = h(2/2, ya que este gap es de primer 
orden en la intensidad z y es más robusto frente a efectos de orden superior. 

Esto muestra que para una configuración de scattering, donde los cables no 
irradiados poseen ocupaciones bien definidas, la conductancia Hall de los aislantes 
topológicos de Floquet quiebra la correspondencia entre el número de estados qui- 
rales y el valor de cuantización de la conductancia Hall. En cambio, la conductancia 
de los aislantes topológicos de Floquet alcanza plateaus en el número de canales 
activos, una vez que se ajustan los empalmes entre los modos del sistema y los 
cables. 

El número de Chern puede no ser suficiente para caracterizar la topología de 
los sistemas forzados periódicamente. En el esquema de scattering con cables no 
irradiados debemos considerar el número de estados que contribuyen al promedio 
temporal de la densidad de estados, para dar cuenta correctamente de la cuantiza- 
ción de la conductancia Hall. Algunos estados topológicos de Floquet permanecen 
ocultos a esta configuración de transporte, aunque podrían manifestarse a intensida- 
des grandes o en otro tipo de medición, como por ejemplo correlaciones de ruido o 
en mediciones que utilicen corrientes dependientes del tiempo. 


4.3. Contribuciones originales 


= Desarrollo de códigos propios para obtener resultados preliminares de las rela- 
ciones de dispersión y conductancia a dos terminales. 


= Explicación de la modulación en la conductancia dos terminales con cables no 
dopados, donde existe un desacople entre los modos normales soportados por el 
cable no irradiado y la muestra irradiada, este desacople proviene principalmente 
del hecho de que el estado topológico tiene una textura de pseudoespín distinta 
a la de los modos propagantes de los cables. 


4. Transporte en aislantes topológicos de Floquet 4.3 Contribuciones originales 





=_ Se aborda la respuesta no local de un aislante topológico de Floquet con cálculos 
de scattering multiterminal, presentado los primeros resultados tanto para el gap 
en el punto de Dirac, como el gap dinámico. Se muestra que la cuantización de la 
conductancia depende del número de canales topológicos activos, una vez que 
se reduce el desempalme entre los modos electrónicos de los cables y la región 
de scattering. 
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Abstract 


El potencial electromagnético de un láser induce gaps y estados de borde en grafeno. A 
medida que baja la frecuencia de la radiación el espectro de Floquet tiene una estructura 
cada vez más compleja. Mientras que el régimen de altas frecuencias es el más simple de 
analizar, en este Capítulo nos enfocamos en los múltiples procesos de emisión y absorción 
fotónica que se habilitan a bajas frecuencias. Para describir estos procesos utilizamos 
Hamiltonianos efectivos que dan cuenta de las propiedades topológicas de cada gap y de 
los estados topológicos asociados y revelan la jerarquía anidada que clasifica a los gaps 
según el orden del proceso fotónico asociado. Cuando las frecuencias son intermedias, 
recurrimos al cálculo explícito del invariante topológico , que nos permite alcanzar distintas 
fases topológicas al variar la intensidad. Este método junto al mapa de transiciones de fase 
topológica da cuenta de la compleja estructura de bandas del Hamiltoniano de Floquet para 


un amplio rango de frecuencias. 


Los aislantes topológicos de Floquet son propuestos como la realización diná- 
mica de los aislantes topológicos donde las propiedades topológicas son inducidas 
por un potencial periódico en el tiempo. La descripción topológica de estos sistemas 
sigue siendo tema de intensa investigación [84, 67, 68, 70, 85], tanto al respecto 
del llenado de las bandas de Floquet, como de los mecanismos de disipación que 
estabilizan el sistema [69, 71, 72] (ver discusión en el Capítulo 4). Por otro lado, 
también se busca un marco teórico para caracterizar la topología del espectro del 
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Hamiltoniano de Floquet. La principal diferencia con respecto a los materiales en 
equilibrio es que el espectro de Floquet es no acotado y periódico en A. Este espec- 
tro también puede obtenerse de los autovalores del operador evolución temporal [36, 
86, 62, 63]. Si observamos la evolución a tiempos discretos nT', con T' el periodo 
del potencial [61], el operador evolución es equivalente al de un sistema estático 
gobernado por un Hamiltoniano efectivo Hp. Entonces, analizando la topología de 
Hof, podemos obtener las propiedades topológicas del sistema de Floquet. Esta 
aproximación es válida siempre que no existan singularidades en el espectro ins- 
tantáneo a tiempos menores que T' [62]. Además, cuando exigimos la conservación 
de alguna simetría extra, puede que necesitemos de un nuevo invariante topológico 
(relacionado con el número de winding) para describir la topología del sistema. Por 
ejemplo si tanto el sistema como el potencial periódico conservan la simetría de 
reversión temporal (TRS), el invariante propuesto en [63] es un invariante 2, para 
describir la topología de estos sistemas y recientemente también se ha propuesto 
un invariante para sistemas de Floquet con simetría quiral [87] 


Cuando la frecuencia es mucho más grande que el ancho de banda del sistema, 
el Hamiltoniano efectivo que describe la evolución estroboscópica del sistema puede 
ser aproximado por una renormalización de la réplica n = 0, que a primer orden sólo 
depende de las réplicas n = +1 [30], 


Her = HSH EGUO)H HAGA. (5.1) 


Esta aproximación deja de ser válida cuando la frecuencia baja y el número de 
procesos fotónicos de absorción/emisión permitidos aumenta considerablemente; o 
bien cuando la intensidad del potencial periódico aumenta. En este caso podemos 
recurrir al método de la zona de Floquet extendida presentado en [36], que se basa 
en considerar un número suficiente de réplicas de Floquet y calcular los números de 
Chern de cada banda de Floquet, para obtener el número de winding. 


En este Capítulo nos basamos en el número de winding para describir la to- 
pología de los gaps y estados topológicos de Floquet asociados. El número de 
winding W (€) es el invariante topológico que describe correctamente el número de 
estados topológicos en cada gap en el espectro de Floquet y podemos utilizarlo 
para frecuencias intermedias y altas. Sin embargo cuando la frecuencia es muy 
baja, el cálculo exacto de los números de Chern se vuelve computacionalmente 
muy costoso, lo que imposibilita obtener los números de winding directamente. En la 
Sección 5.1 derivamos un algoritmo alternativo, a partir del cálculo del número de 
winding, que es válido cuando la frecuencia fotónica es mucho más baja que el ancho 
de banda del sistema A(2 << A y nos permite obtener de manera constructiva las 
contribuciones al número de winding, a medida que incluimos un mayor número de 
réplicas de Floquet. En la Sección 5.2 presentamos el cálculo completo del número 
de winding, analizando las transiciones de fases topológicas que se producen al 
variar la frecuencia e intensidad de la radiación. 


5. Jerarquía de los estados topológicos 5.1 Gaps del orden del número de procesos fotónicos 


5.1. Gaps del orden del número de procesos fotónicos 


En esta sección trataremos el caso de bajas frecuencias y bajas intensidades, en 
donde la aproximación k - p es válida —como veremos luego esta validez se extiende 
hasta energías cercanas a la singularidad de van Hove para grafeno ya que sólo 
necesitamos que el modelo describa bien la topología y textura de pseudoespín del 
sistema. 


Primero analizamos la topología del grafeno irradiado, calculando el número 


de Chern de cada banda de Floquet como la integral de la curvatura de Berry. 


La radiación abre gaps en los cruces (evitados) entre las distintas réplicas de 
Floquet. Las zonas en el espacio k donde ocurren los cruces representan las zonas 
de principal contribución al número de Chern (mayor curvatura de Berry). Luego, 
veremos cómo podemos describir los cruces de las distintas réplicas a través de 
Hamiltonianos efectivos. Utilizando estos Hamiltonianos podemos probar que el 
orden del gap está dado por el número de réplicas involucradas, lo que genera una 
estructura anidada de gaps y estados topológicos de Floquet. 


En el caso de grafeno, donde la relación de dispersión a bajas energías toma la 
forma de dos conos enfrentados, los cruces entre distintas réplicas se producen en 
anillos centrado en k = 0, cuyo radio está determinado las réplicas involucradas (ver 
Fig. 5.1), de forma similar a la estructura anidada de una muñeca rusa. 


5.1.1. Invariantes topológicos de las bandas de Floquet 


Como vimos en la Sección 3.3, el número de winding es la suma de los números 
de Chern de las bandas inferiores y describe correctamente el número de estados y 
quiralidad que se encuentran en un gap en el espectro de Floquet. 


Podemos simplificar aún más este cálculo si tenemos en cuenta que el grafeno 
en ausencia de radiación es un semiconductor de gap cero, cuyas bandas son 
triviales (número de Chern nulo). Recordemos que el número de Chern de la banda 
o es 

i 
Cu=— (var [Vr lUak) - dk 
277 9€ 
1 2 
= —Im 139 Ny Uat lOho Ut) Ka 


TE 


(5.2) 


donde |u.) es la autofunción de Bloch y C el contorno de la zona de Brillouin [21]. 
El número de Chern también puede escribirse como la integral de la curvatura de 
Berry en la primera zona de Brillouin 


1 


A 
“2 JBz 


Dor: dSy, (5.3) 


Y*X/ 


gap de 1* orden 


Y 





Fig. 5.1: Esquema de las réplicas 
de Floquet (sólo de las ramas posi- 
tivas), donde los cruces en el gap 
dinámico se muestran como cortes 
en la relación de dispersión. Los 
cruces de orden superior envuel- 
ven a los de orden inferior de ma- 
nera anidada lo que recuerda a las 
mamushkas. 
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*La densidad de puntos en el es- 
pacio k que debemos conside- 


rar crece notoriamente porque 


los gaps de orden superior (cu- 


yo origen es el cruce de réplicas 
de orden superior) son cada vez 
más pequeños y su curvatura 


de Berry muy pronunciada, ade- 


más de que múltiples réplicas 


deben ser incluidas en el cálcu- 


lo numérico y crece la matriz del 
Hamiltoniano de Floquet. 
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donde I'., es explícitamente la curvatura de Berry. 


Como vemos en la Ec. (5.4), los picos en la curvatura de Berry provienen de las 
zonas donde las energías están cuasi-degeneradas. Cuando sometemos al sistema 
a un potencial electromagnético periódico en el tiempo, obtenemos un espectro 
de Floquet cuyas réplicas se solapan unas con otras. Estos solapamientos son 
levantados (se remueven las degeneraciones) cuando la intensidad de la radiación 
es mayor que cero. Por lo tanto a intensidades muy bajas la principal contribución 
a la curvatura de Berry proviene de las zonas en el espacio k donde las energías 
están cuasi- degeneradas. 


El cálculo del número de winding involucra un espectro de Floquet truncado 
hasta un cierto número de réplicas, donde tenemos un número finito de bandas de 
Floquet. Si truncamos el espectro de Floquet entre las réplicas -M y M (con un 
total de 2M + 1 réplicas), obtenemos dos de bandas de Floquet por cada réplica 
considerada. Llamemos C., al número de Chern de la banda «a de Floquet, donde «a 
está definido entre -(2M +1) y (2M + 1). El esquema de zonas repetidas [36] nos 
dice que el número de winding W., del gap por encima de la banda «a está dado por 
la suma de los números de Chern de las bandas inferiores 


Q 


Wísa)= Y; Cp, (5.5) 
B=-(2M+1) 
siempre que el número de réplicas incluidos en el cálculo sea suficientemente grande 
para que el cálculo converja. 


En la Fig. 5.2 vemos la densidad de estados que se calculada cerca del borde 
de un plano semi-infinito de grafeno irradiado, donde, por ejemplo para M = 2, los 
números de Chern de las bandas de Floquet C,, varían hasta que Ca = +4 para 
las bandas centrales. Como resultado, el número de winding W¿ = 0 para el gap 
en ep =0 y W;y2 = -4 para el gap en ey, = h02/2 reflejan el número de estados de 
cada gap. En el primer caso Wo = 0 nos dice que el número total de estados es cero, 
aunque el gap está cruzado por dos estados de quiralidades opuestas, en cambio 
W1/2 = -4 nos dice que dentro del gap hay 4 estados con la misma quiralidad. 


A pesar de utilizar el formalismo de las zonas repetidas de Floquet y el Hamil- 
toniano truncado con un número finito de bandas, debemos recordar que nuestro 
objetivo es dar cuenta de la topología del espectro de Floquet completo (periódico y 
no acotado). La mejor aproximación al espectro de Floquet completo, son las bandas 
centrales del Hamiltoniano de Floquet truncado (bandas « = 0 y a = 1 en la Fig. 5.2). 


A medida que la frecuencia se hace más baja, el cálculo explícito resulta numéri- 
camente muy costoso * . En esta sección haremos uso de otra aproximación que 
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Fig. 5.2: Esquema de zonas re- 
petidas de Floquet Las bandas 
de Floquet se muestras para un 
plano semi-infinito de una red he- 
xagonal irradiada, truncado entre 
(-M, M) réplicas con M = 2 
Los números de Chern C; de ca- 
da banda se muestran dentro de 
las 10 bandas de Floquet, mien- 
tras que los números de winding 
W(eo) y W(e1/2) aparecen den- 
tro de los gaps alrededor de ey y 
€1/2, respectivamente. Los núme- 
ros de winding caracterizan el to- 
tal de estados quirales dentro del 
gap, siendo W (20) = 0, por lo que 
en £o hay dos estados de quirali- 
-10 dades opuestas y W(e1/2) = -4 
por lo que en e,y2 hay cuatro esta- 
dos con la misma quiralidad. Figu- 
ra adaptada de [88]. 
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sirve para calcular de manera constructiva las contribuciones al número de winding 
en el régimen de bajas frecuencias. 


Para calcular las contribuciones al número de winding recordemos que este 
invariante topológico es la suma de los números de Chern cuando calculamos el es- 
pectro del Hamiltoniano de Floquet truncado a un número suficiente de réplicas [36]. 
El número de Chern de la banda «a de Floquet, en el límite de intensidades bajas, 
es la suma de las contribuciones del espacio k donde el espectro se encuentra 
cuasi-degenerado debido a los cruces con otras réplicas de Floquet. Los cruces de 
la banda a: pueden ser tanto por encima, con la banda («a + 1) cuyas contribuciones 
Cp. provienen de los puntos k,, ,; como por debajo, con la banda a — 1, dando 
contribuciones esa que provienen de los puntos posiblemente distintos k,, . Luego 
obtenemos 


CELO ADO (5.6) 
Pp p 


siempre que no sea la primera o la última banda, es decir a+ + -(2M +1) 0a + 


(2M + 1), en cuyo caso no debemos incluir en la suma de E 0 


Debido a que las contribuciones dentro de un cruce evitado tienen distintos 
signos entre la banda de arriba y la de abajo 


Cra = ari) (5.7) 
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podemos deducir que el número de winding, al ser la suma parcial de los números 
de Chern, tendrá ambas contribuciones de manera que se cancelen entre ellas, 
excepto por las últimas contribuciones 


(2 008 + y dea] (5.8) 
P 


Vemos que el número de winding queda definido solamente por las contribucio- 
nes superiores (c,!,) al número de Chern de la banda inferior, y podemos descartar 
el superíndice. 


El espectro de Floquet completo (Hamiltoniano de Floquet no truncado) está 
dividido en dos bandas que se repiten a intervalos de A(2. Estas bandas tienen gaps 
en go y 1/2. Caracterizamos los números de winding con Wo y W y, respectivamente. 
Para obtenerlos debemos sumar las contribuciones que llamaremos c,, ¿, con PB = 0 
y f6=172. 

El siguiente paso es encontrar los puntos k,, g, las contribuciones están bien 
separadas en el espacio k y podemos indexarlas distinguiendo las contribuciones 
dentro del gap en ey = 0 y en eyy2 = hQ/2. Por la simetría electrón-hueco y la 
relación de dispersión cónica del grafeno podemos intuir que las zonas de cuasi- 
degeneración serán los anillos de radio k,, ¿ centrados en k = 0 para cada valle 
K, y K_. Para ver esto más claramente, escribimos el Hamiltoniano de Floquet en 
la base (1/42, +exp (i0z,) II que diagonaliza los subespacios de cada réplica 
Ho + nhQTy obtenemos 


Ad) 7) 0 
HE (k,0)=|-- VO)? E) VO) -l, (5.9) 
o V(0%)t Ax) 


donde 
ES k 0 
A (k) =vp he 0 + nup2ko1 (5.10) 
a ' 1 lA 
V (0%) = nkpe e En se ) (5.11) 
—e —. 


donde el? es una fase trivial entre los vectores de base, de ahora en más A = 0, 
6, es el ángulo entre k y el eje x y 2upky reemplaza AM para hacer evidentes los 
cruces en el espacio k. 
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Los valores del módulo de k donde se cruzan las réplicas son k,, ¿ = 2ko(p+ B), 
donde $ = 0 o 1/2, según estemos calculando las contribuciones al gap en £0 O e. 
En la Fig. 5.3 vemos esquemáticamente los cruces entre las distintas réplicas que se 
producen en los puntos k,, ¿. Además los cruces de orden más alto envuelven a los 
de orden inferior, lo que le da una jerarquía indexada por el orden de los procesos 
fotónicos que originan los gaps. 























5.1.2. Múltiples procesos fotónicos en bajas frecuencias 


Podemos deducir ahora, a partir de la Ec. 5.9 los Hamiltonianos efectivos que 
describen los cruces entre distintas réplicas. Teniendo en cuenta la analogía entre 
el teorema de Floquet y el teorema de Bloch, podemos definir una zona donde 
restringir la coordenada periódica. Para el Hamiltoniano de Bloch podemos restringir 
el cuasimomento k a la primera zona de Brillouin y para el Hamiltoniano de Floquet 
podemos definir la primera zona de Floquet (FZ) como la región entre +RA(2/2 de 
cuasi-energías. Debemos diferenciar los cruces entre réplicas que ocurren en el 
centro de la FZ (so = 0), de los que ocurren en el borde de la FZ (£1/2 = h£2/2). 


Gaps en el borde de la zona de Floquet 


El caso más simple es el que describe el gap dinámico a primer orden. Para 
describir este cruce sólo hace falta retener las réplicas n = 0 y n= 1 








k +2k0 0 nkgeió | nkgei% 
- 0 —k+2kp9 —nkpei% | —nkyeió 
HO L(k) = up E eN dl E ó a (5.12) 
0 =NKo 
nkpe | —nkpelór 0 —k 








Fig. 5.3: (a) Réplicas de Floquet 
para un sistema tipo Dirac, en ne- 
gro para la réplica n = O y en 
rojo (azul) para las réplicas con 
n > 0 (n < 0). Los cruces entre 
las distintas bandas n y m se pro- 
ducen a dos energías no equiva- 
lentes (módulo AQ), ey = O cuan- 
don+m=0 y Ex1/2 = +h0/2 
cuando n + m = +1. Los distintos 
cruces originan gaps de distintos 
órdenes según las réplicas involu- 
cradas. (b) Los gaps en ey son de 
orden par n — m (excepto en el ca- 
so que n = m = 0 donde el gap 
proviene de la renormalización de 
las bandas). (c) Los gaps en e1/2 
son de orden impar n — m. [88]. 
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tla rotación está hecha con la 
matriz 


O sin0 cos0 
R(O) = (o -cosó sin0 _ 
1.0 0 


que satisface la identidad 


R *0Rh».s =x hos. 


FE proceso de decimación con- 


siste en reducir un problema de 


autovalores Hu = eu reducien- 
do la matriz con el siguiente al- 


goritmo [42] 
1 - .e, 
(Po 2) (25), 
o vtr! 


donde 


V=VpV 
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donde el recuadro central contiene los subespacios degenerados cuando k = kg 1/2 = 
ko y la tilde nos recuerda que está definido en la base rotada. El Hamiltoniano efectivo 
Ak, p, 5) en este caso es 


(5.13) 


Ñ -k+2k0 —nkge Pr 
| o ES ) 


—nkyetr k 


y en el caso general podemos escribirlo como Á¿*(k,p, 8) = hy ¿(k) - 0 + egZ, 
donde la tilde indica la base rotada, a es el vector de las matrices de Pauli y 7 es la 
matriz 2 x 2 identidad. 


Para calcular la contribución al número de Chern debemos volver a la base no 
rotada y reescribiendo la Ec. (5.3) obtenemos 


1 A A A 
cpo=35 $ hipo (Oh, «0%, hs) (5.14) 
donde h,.s = h,.2/|h,,5| está definido en la base no rotada. Esta expresión tiene que 
ser reformulada en la base rotada ya que los vectores base dependen de 0,,, lo que 
implica una rotación * del vector efectivo hos = R(0,)h,g. Usando coordenadas 
polares obtenemos 


1 1 
cpa=>32 f hpa-(Oohpa Oh, 5) 7, ao 


1 pS y > 1 
> e Ss hp,g* (Oh x Oxhr,g) z 


hn dkd0 





(5.15) 


1 ll o dc - 1 
+— | ho: (R*9Rh, gx Oxhp 5) —dhd6. 
4 ú : : Rp, al? 
En la Ec. (5.15) hacemos uso de la rápida convergencia del integrando para extender 
la integral a todo el espacio k. 


A la misma energía e¡/2 el próximo cruce entre las réplicas n =—1 y n = 2 será 
en ki 1/2 = 3 ko. En este caso debemos decimar + las réplicas n = 0 y n= 1, de 
manera que se renormalizan los acoples entre las distintas réplicas, generando un 
gap de tercer orden (O(1?)) en el parámetro de acople electrón-fotón. Si definimos 
en función de la Ec. (5.10) la función de Green 


>, 1 
G"(e, ko.8) = ——==—, 5.16 
(E, Kp,a) <L - Í”(kp,5) ( ) 
entonces los acoples entre la réplica m =-—1 y n = 2 renormalizados serán 
Va -1(04) = V(01)G (e1/2, k11/2)V (0)GCeyjo, 1112) V (01), (5.17) 


donde el subíndice (2,-1) indica el orden de las réplicas que se acoplan. Además, 
las funciones de Green deben ser evaluadas a la energía e 2 y en el radio del vector 
de onda k; 1/2. 
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Es fácil ver que en la base en que 4? es diagonal, (Ec. (5.16)), dado que las 
funciones de Green están definidas como las inversas de matrices diagonales, 
entonces también G”(+<) será diagonal. Esto junto al prefactor de V (9) de la 
Ec. (5.10), hace que los acoples Vn.m(0x) sean proporcionales a ql Por lo tanto, 
el siguiente cruce en el espacio k será en ka 1/2 = 5 ko entre las réplicas m = -2 y 
n = 3 y abrirá un gap de orden 1? ya que debemos decimar las réplicas n = —1,0,1,2. 
Siguiendo este análisis, cada vez que observamos un nuevo cruce, el orden del gap 
aumenta en 17. 

Ahora podemos calcular explícitamente los vectores h,, ¡/2 que describen los 
cruces a energías e /, serán 


ho1/2(k) = -né1ko 1/2 + (-k + ko1/2) 2 
- 4 3 E 
hy 1/2(k) = Eskr 256 + (-k + k1.1/2) z 


E ES 80 e 
ha 1/2(k) = 1 Esk21/23965 sl (-k js ka,1/2) z 


(5.18) 


2 


donde, É,, = cos (n0y)4 + sin (n0z)4 gira n veces en el plano («, y) al rededor del 
eje z mientras k recorre la zona de Brillouin. 

Utilizando los vectores de la Ec. (5.18) podemos calcular las contribuciones 
respectivas al número de winding según la Ec. (5.15) y vemos que 


Co1/2 = 1 
1/2 =9 (5.19) 
C2,1/2=5 


..) 


para cada valle K, y K_. Los resultados numéricos del modelo tight-binding de una 
muestra semi-infinita de grafeno irradiado de la Fig. 5.4 están en perfecto acuerdo 
con los resultados analíticos. En el panel (a) vemos las réplicas diferenciadas por 
su color, que indica el peso en escala logarítmica sobre la réplica n = 0. Los cruces 
entre las réplicas generan los gaps de orden cada vez mayor (y ancho cada vez 
menor), de los que en (a) sólo se distinguen los primeros 2 gaps. En el panel (b) 
expandimos la escala de energías al rededor de e, para tener poder aumentar la $ Esta energía de corte £ está 


resolución cerca del gap, hasta una energía de corte ¿$ . limitada por la estabilidad numé- 
rica y la precisión utilizada en 


los cálculos. 
Gaps en el centro de la zona de Floquet 
Podemos aplicar el método anterior para encontrar los Hamiltonianos efectivos 
que describen los cruces a energía ey, con la única excepción del cruce en kg = 0 


ya que en este caso, la rotación de la base no está bien definida (el ángulo 0, no 
está bien definido en k = 0). 
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Fig. 5.4: (a) Densidad local de 
estados cerca del borde de un 
plano semi-infinito de una red he- 
xagonal irradiada con A(2 = 0,2t 
y z = 0,05 En escala logarítmica 
de color el peso sobre la LDOS 
promedio peage(k, e), resuelta en 
energía e y en el momento cris- 
talino k paralelo al borde. (b) La 
escala de energía está expandida 
logarítmicamente respecto del cen- 
tro del gap en e,/2, hasta una ener- 
gía mínima de corte € = 610” *". 
En (a) y (b) observamos que las 
réplicas de Floquet de orden alto 
se corresponden con un peso so- 
bre la LDOS media de orden más 
alto, mientras que los gaps de or- 
den superior se distinguen en el 
panel (b) junto con los estados de 
borde asociados. Figura adaptada 
de [88]. 
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Para el cruce en k = O podemos renormalizar la réplica central n = 


Hamiltoniano de Floquet 


H(k) V(0%) 0 
Hy(k)=|V (0) HUk) V(0x) |, 
0 V(0)? Ek) 


y mediante el algoritmo de decimación obtenemos 


HÉE(k,0,0) =H (ke) + V (9) 1G*(0,0)V (07) 
+ V(0,)G"*(0,0)V(01)*, 


que puede ser escrito como el vector 


hoo(k) = -2 kot + k2. 
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O del 


(5.20) 


(5.21) 


(5.22) 


Utilizando la Ec. (5.14) que hace uso del vector h,, ¿ no rotado, podemos encontrar 


que coo = -1/2 para cada valle K, y K_, lo que da un total de -1 estados. Esta 


es la única contribución negativa al número de Chern y es de esperar que este 
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estado de borde una los valles K,, como podemos comprobar en la Sección 3.3 
y en la Fig. 3.7, incluso para frecuencias e intensidades muy altas. Para el caso 
de frecuencias muy altas (Aá(2 > 1,5t) no hay cruces de orden superior en ey. Sin 
embargo, para frecuencias bajas, el cruce siguiente en k; y también abre un gap de 
orden 7”, por lo que las contribuciones coo y c1,0 deben considerarse en conjunto. 

En el punto ko o = O no podemos utilizar la base rotada, pero en el caso general 
de los cruces en k,,0 (p > 1) podemos tomar los Hamiltonianos efectivos de la base 
rotada. Debemos decimar la réplica n = O para el cruce en k; y entre las réplicas 
m =-1 y n= 1, por lo que el gap será de orden »?. El siguiente gap entre las réplicas 
m =-2 y n = 2 será de orden y* y así siguiendo. Evaluando los Hamiltonianos 
efectivos encontramos que 


a a 1 Ñ 

h1o(k) = Eh + (=k + k10)2 
a al 

hoo(k) = 1 Erkoo— + (=k + k20) 2 
2,0(k) 4K2.075 ( 2,0) Eos) 


a , 9 Ñ 
h30(k) = 9 Eska0z5 0 + (—k + k30) 2 


..) 


donde A(k,p,B) = h,¿(k) - 0 + egT. A partir de los vectores de la Ec. (5.23) 
podemos utilizar la Ec. (5.15) para calcular la contribución al número de Chern de 
cada cruce y encontramos que 


C1,0 = 2 
cap=4 

del (5.24) 
C3,0 =6 


En la Fig. 5.5 vemos las bandas para una muestra semi-infinita de grafeno 
irradiado. En el panel (a) notamos que (i) el estado que une los valles K, y K_ no 
tiene la quiralidad opuesta, ya que las contribuciones cy o = -1 y 2c1 0 = 4 son ambas 
contribuciones a un gap de orden 1?. Esto hace que el número de winding efectivo 
en este orden sea la suma de ambas contribuciones por lo que obtenemos un total 
de 3 estados, (ii) los gaps en ko o y k1,0 son del mismo orden, siendo el primero del 
doble del ancho que el segundo, mientras que (ii) los gaps de órdenes más altos 
son apenas distinguibles. En el panel (b) expandimos logarítmicamente la escala 
de energías al rededor de ey hasta una energía de corte £ para poder discernir los 
gaps y los estados de orden superior, que concuerdan con la Ec. (5.24). 

Tanto para los gaps en el borde de la zona de Floquet como centro de la zona 
de Floquet (excepto el gap en el punto de Dirac k = 0), podemos expresar los 
Hamiltonianos efectivos en la base rotada como 


hay, 9 (10) = 9 al ko [cos (s20J3 + sin (sP01)9] + (h+ só ho)2 (5.25) 
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Fig. 5.5: (a) Densidad local de 
estados cerca del borde de un 
plano semi-infinito de una red he- 
xagonal irradiada con 22 = 0,2% 
y z = 0,05 En escala logarítmica 
de color el peso sobre la LDOS 
promedio peage(k, e), resuelta en 
energía e y en el momento cris- 
talino k paralelo al borde. (b) La 
escala de energía está expandida 
logarítmicamente respecto del cen- 
tro del gap en eo, hasta una ener- 
gía mínima de corte £ = 3£10”*. 
En (a) y (b) observamos que las 
réplicas de Floquet de orden alto 
se corresponden con un peso so- 
bre la LDOS media de orden más 
alto, mientras que los gaps de or- 
den superior se distinguen en el 
panel (b) junto con los estados de 
borde asociados. Figura adaptada 
de [88] 


Para mantener la expresión en 
su versión más simple, descar- 
tamos las contribuciones de or- 
den O(n?) a la componente z 
del vector h,, ¿. Estas correccio- 
nes provienen de la renormali- 
zación de las réplicas y no modi- 
fican el resultado de la integral 
de la contribución al número de 
Chern. 


Por ejemplo para el gap en 
eo, los valores que toma s), son 
s=2,58% =4,s% =6 y así si- 
guiendo; mientras que los valo- 


1/2 1/2 
res que toma s;/? son s//” = 1, 


si? =3, si? =5, 
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a) o. 





donde sí = 2(p+ 8), 0 = tan” *(k,/k) y al es un factor numérico f . 


Usando la Ec. (5.15), donde el primer término es nulo (excepto para ho 0), 
obtenemos que la contribución al número de Chern es 


B 
S 1 
Cp,B = E 1 + , 


(atacó) sa 





(5.26) 


la cual representa el caso general para cualquier orden de gap tanto en ey como en 
€1/2- Si tomamos el límite en que 7 > 0, las aproximaciones hechas hasta ahora se 
vuelven exactas y las contribuciones son 


Gui Eh (5.27) 


La Ec. (5.27) es el resultado central de esta Sección y puede entenderse fácilmente 
si recordamos que la contribución al número de Chern debe ser igual al número de 
veces que el vector hos gira en el plano (x, y) en torno al eje z mientras k explora 
la primera zona de Brillouin | En general sf** = 2 + sf, como podemos comprobar 
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en las Figs. 5.4 y 5.5, donde cada gap de orden superior tiene +2 estados. 


5.2. Mapa de transiciones de fases topológicas 


El método de los Hamiltonianos efectivos de la Sección 5.1 describe correcta- 
mente las contribuciones al número de Chern, incluso más allá del rango de validez 
del método k - p en el que está basado. Esto se debe a que el método k - p describe 
correctamente la topología (su textura de pseudospin) de las bandas a energías por 
debajo de la singularidad de van Hove |:| < 1t = 2,7eV. 


5.2.1. Frecuencias intermedias e intensidades altas 


Cuando la frecuencia es más alta AQ - A, donde A es el ancho de banda, hay 
menos solapamientos de las distintas réplicas de Floquet, lo que significa que hay 
menos procesos inelásticos debido al alto costo energético. Este régimen es por lo 
tanto más simple de describir que el caso tratado en la Sección anterior (5.1), sin 
embargo surgen otras dificultades. Cuando aumentamos la frecuencia el número 
de réplicas bien descritas por la aproximación de bajas frecuencias rápidamente 
decrece y debemos recurrir al método tight-binding y hacer el cálculo explícito. Este 
cálculo se basa en el Hamiltoniano de Floquet truncado a un número de réplicas 
entre +M * . Entonces, debemos incluir al menos todas las réplicas que se solapen 
con la primera zona de Floquet, de modo que Min = int(A/A0) TT. 

A medida que aumenta la intensidad, los procesos inelásticos de orden más 
alto se incrementan y el peso de las autofunciones del Hamiltoniano de Floquet se 
distribuyen en un mayor número de réplicas. Como consecuencia debemos tomar 
un número mayor de réplicas para describir correctamente el espectro dentro de 
la primera zona de Floquet. El número de réplicas suficientes estará dado por el 
número de réplicas necesarias para que el cálculo del número de Chern de las 
bandas centrales converja. Por ejemplo en la Fig. 5.2 debemos incluir 5 réplicas de 
Floquet, a pesar de que A/A < 3. Vemos cómo cambia el número de estados por 
gap (en los gaps fuera de la primera zona de Floquet (-A£2/2,h(0/2)) a medida que 
nos acercamos a la réplica central. Las bandas centrales son la mejor aproximación 
al espectro del Hamiltoniano de Floquet no truncado, dentro de la primera zona de 
Floquet. 

El aumento de procesos inelásticos favorecido por el aumento de la intensidad 
de la radiación, conllevan a nuevas transiciones de fase (ver Figs. 5.6 y 5.7), como 
resultado de la distribución de las autofunciones en muchas réplicas de Floquet. 


5.2.2. Transiciones de fase topológicas al variar frecuencia e intensi- 
dad de la radiación 


En la sección anterior, observamos nuevas contribuciones al número de Chern 
cada vez que incluíamos un proceso fotónico de orden más alto, asociado a un gap 


5.2 Mapa de transiciones de fases topológicas 


**Como vimos en la Ec. (5.14) el 
cómputo del número de Chern 
depende de los autovectores 
del Hamiltoniano integrados so- 
bre (un muestreo) de la zona de 
Brillouin, por lo que es necesa- 
rio disminuir lo más posible la 
dimensión del Hamiltoniano. 


Debemos recordar que en 
nuestro caso, el ancho de ban- 
da A disminuye en función de 
la intensidad de radiación como 
6t Jo(z), ver Apéndice A.1. 
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Fig. 5.6: Número total de estados 
de borde dentro del gap en gp = 0, 
representados por el número de 
winding W (e0) (en escala de co- 
lor) al variar la frecuencia Aí2 e 
intensidad z de la radiación. A ba- 
jas intensidades, al disminuir la fre- 
cuencia nuevos solapamientos en- 
tre réplicas de Floquet ocurren ca- 
da vez que Af? = 3t/m con m en- 
tero (como se indicado) y marcan 
una línea de transiciones topoló- 
gicas. A intensidades mayores los 
hoppings entre réplicas crecen y 
nuevas transiciones ocurren, mien- 
tras que los hoppings dentro de 
la misma réplica disminuyen hasta 
cero en el valor crítico zp, 1, el pri- 
mer cero de Jo(z). Figura adapta- 
da de [88]. 
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5.2 Mapa de transiciones de fases topológicas 


de ese orden. En cambio en el régimen de altas frecuencias, el número de réplicas 
que se solapan está dado por la frecuencia de la radiación y podemos calcular 
explícitamente los números de Chern cuya suma es el número de winding, a medida 
que variamos la frecuencia e intensidad de la radiación. 

A intensidad nula (z = 0), el ancho de banda A = 61 de la n-ésima réplica (n > 0) 
se encuentra a energías entre nñáQ + 3t, por lo que se solapará con las demás 
réplicas de Floquet en ey = O cuando 


AO < 3t/n, (5.28) 


y en e1/2 = hQ/2 cuando 


RO <3t/(n— >) (5.29) 

A medida que bajamos la frecuencia, más réplicas se solapan y los nuevos 
solapamientos producen posibles transiciones de fase. Cuando la intensidad es 
baja, el ancho de banda de las réplicas se reduce aproximadamente como tJo(2) = 
t(1 -— 22), por lo que en las Figs. 5.6 y 5.7 vemos líneas de transición de fase como 
parábolas hacia abajo, que nacen a las frecuencias de indicadas por la Ec. (5.28) 
en la Fig. 5.6 y por la Ec. (5.29) en la Fig. 5.7. Cabe aclarar que no todos los 
solapamientos producen una transición, por ejemplo en la Fig. 5.7 a AQ = 3t(2 — 5), 
pero sí producen un cierre y reapertura del gap correspondiente a medida que 
bajamos la frecuencia. 





2.5 
2 2d 


El análisis anterior se basa en el hecho de que los acoples entre la n-ésima 
réplica y la m-ésima réplica son proporcionales a t J,,_.. (2), que a intensidades 
bajas «= t 2l"-"l. Luego, los principales acoples son entre sitios que pertenecen a la 
misma réplica. A medida que aumenta la intensidad, la dependencia de funciones 
de Bessel JJ, (2) hace que esta aproximación deje de ser válida y el acople entre 
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5.3 Contribuciones originales 





3t/(1-3)> 6 


h0/t 


3t/(2-5)> 2 
5 
3t/(8=3)1 427 
3t/(4 > 
20,1 0.5 1 1.5 2 42.5 
Z 20,1 


sitios vecinos pero de réplicas aledañas pase a ser incluso mayor que los hoppings 
dentro de la misma réplica, como vemos en la Fig. 5.8. A medida que aumenta la 
intensidad, los efectos de introducir una nueva réplica en el cálculo, se extienden 
mucho más allá del ancho de banda de esa réplica y se generan nuevas transiciones 
de fase. Estas son las líneas de transición de pendiente positiva en las Figs. 5.6 y 
5 

Un caso extremo ocurre cuando se anulan los hoppings dentro de la misma 
réplica, proporcionales a t Jo(7). El desvanecimiento del acople dentro de la misma 
réplica se da cuando la intensidad es zp,1 = 2,4048, el primer cero de la función de 
Bessel Jo(z). A la intensidad zo 1 una línea vertical en las Figs. 5.6 y 5.7 marca la 
aparición de transiciones de fase en ambos gaps para todas las frecuencias. A esta 
intensidad una pequeña desviación de la polarización circular hacia polarización 
elíptica puede generar transiciones de fase en el punto de Dirac [38], por lo que es 
una región altamente inestable del mapa de fases. Durante la finalización de esta 
tesis, un trabajo reciente [89] estudia las transiciones de fase con una aproximación 
de Hamiltonianos efectivos para frecuencias altas y logran explicar transiciones que 
se observan para intensidades extremadamente altas y frecuencias altas. Además 
presentan un mapa de transiciones para un rango más grande de parámetros que 
en esta tesis hemos descartado por consideraciones físicas de factibilidad. El ruido 
numérico que se observa para bajas intensidades en las Figs. 5.6 y 5.6 se debe a la 
inestabilidad del cálculo de la integral de la curvatura de Berry cuando la amplitud de 
los gaps es muy pequeña. Además, puede observarse ruido cerca de las transiciones 
de fase ya que cada línea de transición marca el cierre de un gap. 


5.3. Contribuciones originales 


= Revelar la jerarquía anidada de gaps y de estados topológicos asociados a estos 
gaps, que están paramétricamente separados por su peso o contribución a la 





Fig. 5.7: Número total de estados 
de borde dentro del gap en e,y2 = 
hQ/2, representados por el núme- 
ro de winding W(e/2) (en escala 
de color) al variar la frecuencia A 
e intensidad z de la radiación. A 
bajas intensidades, al disminuir la 
frecuencia nuevos solapamientos 
entre réplicas de Floquet ocurren 
cada vez que hQ = 3t/(m - 1/2) 
con m entero (como se indicado) y 
marcan una posible línea de tran- 
siciones topológicas. A intensida- 
des mayores los hoppings entre 
réplicas crecen y nuevas transicio- 
nes ocurren, mientras que los hop- 
pings dentro de la misma réplica 
disminuyen hasta cero en el va- 
lor crítico zo,1, el primer cero de 
Jo(z). Figura adaptada de [88]. 


1 — Jolz) 
0.8 —yi(z) 
0.6 — J,(Z) 





Fig. 5.8: Funciones de Bessel de 
orden 0, 1 y 2. El primer cero de 
Jo(2) se indica por zo 1 > 2,4048. 
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densidad de estados promediada en el tiempo. El método de los Hamiltonianos 
efectivos presentado en la Sección 5.1 tiene la ventaja de dar cuenta de los 
primeros órdenes de gaps y estados topológicos, que son las contribuciones más 
importantes a la densidad de estados promedio. 


El mapa de transiciones de fase topológica de Floquet en la Sección 5.2 muestra 
cómo obtener distintos números de estados tanto en el gap dinámico como en el 
gap del punto de Dirac, mostrando la riqueza del la topología de las bandas de 
Floquet cuando la frecuencia es baja. 


Los dos métodos presentados para describir la topología de las bandas de 
Floquet logran describir las múltiples transiciones topológicas de un sistema 
tipo Dirac para frecuencias de bajas a intermedias y para un amplio rango de 
intensidades. 


CAPÍTULO 6 
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Abstract 


En este capítulo se presenta un resumen de las principales contribuciones derivadas del 


trabajo de tesis. Por otra parte, presentamos también una breve discusión de algunas de los 


problemas que aún quedan por resolver y que están siendo investigados actualmente. 


6.1. Principales Contribuciones Originales 


= En el Capítulo 2 presentamos un estudio exhaustivo sobre los efectos de la 
radiación de en nanocintas de grafeno. Primero, estudiamos grafeno bulk bajo 
la incidencia de radiación de polarización circular y vemos que abre gaps en la 
densidad de estados promediada en el tiempo y que la polarización lineal no 
abre gaps aunque sí disminuye la densidad de estados media. En particular, 
demostramos que existen distintos mecanismos de transiciones electrónicas 
habilitados o inhabilitados según la dirección de la polarización y la terminación de 
la nanocinta. El interjuego entre la dirección de la polarización lineal de un 
láser y el confinamiento electrónico dentro de las nanocintas de grafeno 
nos da un parámetro externo de control sobre gaps fuera del punto de 
Dirac. La dirección de polarización respecto del borde de la nanocinta es el 
parámetro de control que permite o inhibe las transiciones intra-modo o inter- 
modo en cintas de terminación zigzag o armchair. 


La radiación de polarización circular rompe la simetría de reversión temporal y 
abre gaps tanto en el punto de Dirac como en el borde de la zona de Floquet y 
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ambos gaps son no triviales. En el Capítulo 3 presentamos las primeras simu- 
laciones atomísticas que evidencian los gaps y los estados topológicos de 
borde y presentamos la primera solución analítica a los estados topológi- 
cos de Floquet en grafeno irradiado. Estos resultados están apoyados por el 
cálculo explícito de las propiedades topológicas de las bandas de Floquet que 
muestra la existencia de estados de borde quirales. Verificamos la quiralidad de 
los estados de borde con cálculos de transporte local dentro de una nanocinta 
irradiada. 


= En el Capítulo 4 presentamos los primeros cálculos de transporte multiter- 
minal de un aislante topológico de Floquet. La conductancia Hall calculada 
en un esquema de scattering muestra que eventualmente se puede llegar a 
una cuantización donde los estados de borde relevantes son aquellos con 
peso significativo sobre la densidad de estados media. Esto distingue a los 
aislantes topológicos de Floquet de los aislantes topológicos estáticos. En ambos 
casos la topología del bulk se refleja en la cantidad y quiralidad de los estados 
de borde (correspondencia volumen/superficie), aunque no todos los estados 
participan en el transporte multiterminal con cables no irradiados. 


= Finalmente, vemos que más allá de las aproximaciones a primer orden, existen 
estados topológicos en una estructura de bandas de creciente complejidad a me- 
dida que bajamos la frecuencia de la radiación. En el Capítulo 5 analizamos en 
detalle el régimen de bajas frecuencias y damos un método constructivo 
para obtener las propiedades topológicas de las bandas de Floquet. Cuan- 
do la frecuencia es muy baja, la zona de Floquet es muy pequeña y el espectro de 
Floquet involucra pliegues del espectro original, con mayor número de pliegues 
a menor frecuencia. Por lo tanto en el límite A << A es altamente desafiante 
conseguir una caracterización completa de la topología de las bandas de Floquet. 
En el régimen de bajas frecuencias obtenemos Hamiltonianos efectivos, que 
caracterizan las contribuciones al invariante topológico y describen el cruce entre 
las distintas réplicas de Floquet. Con este método constructivo vemos que se 
agregan estados de borde a medida que consideramos un mayor número 
de réplicas. Estos gaps se ordenan en una jerarquía indexada por el núme- 
ro de procesos fotónicos que originan el gap, en una estructura anidada 
de gaps y estados topológicos asociados. El régimen de frecuencias inter- 
medias es analizado con el esquema de zonas repetidas de Floquet [36], con 
el que podemos calcular explícitamente el invariante topológico de cada banda 
de Floquet. A partir de este esquema obtenemos un mapa de transiciones de 
fase topológica para un amplio rango de frecuencias e intensidades de radiación 
de polarización circular sobre grafeno. Con los dos métodos de Hamiltonianos 
efectivos para frecuencias bajas y el cálculo directo del invariante topológico para 
frecuencias intermedias, caracterizamos de manera más precisa la topología de 
las bandas de Floquet del grafeno irradiado. 
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Los resultados obtenidos durante el desarrollo del trabajo de tesis han dado 
lugar a 7 publicaciones en revistas internacionales, que han colectado más de 100 
citas el momento de la presentación de esta tesis. Esto refleja el impacto de nuestro 
trabajo como también, el gran crecimiento del área y el interés hacia las propiedades 
de transporte en grafeno y aislantes topológicos de Floquet. 


6.2. Perspectivas en el campo de esta tesis 


6.2.1. Clasificación topológica de sistemas periódicos 


La primera aproximación para clasificar la topología de los aislantes topológicos 
de Floquet es utilizar los invariantes topológicos heredados del caso estático. Si 
observamos la evolución a estroboscópica del sistema a tiempos nT', donde T' es 
el periodo del potencial, esta evolución puede ser descrita por un Hamiltoniano 
efectivo independiente del tiempo. Podemos calcular el Hamiltoniano efectivo para 
frecuencias altas [30] y describir la topología de este Hamiltoniano al igual que en el 
caso estático. 

En un caso más general, para cualquier frecuencia, podemos definir un invariante 
que caracteriza la topología de sistemas bidimensionales forzados periódicamente 
mediante el número de winding [36]. El número de winding es equivalente a la 
suma de los números de Chern de las bandas de un espectro de Floquet truncado, 
sin embargo también puede obtenerse de los autovalores del operador evolución 
temporal. 

Cuando además imponemos la conservación de algunas simetrías, es posible 
que sea necesario un nuevo invariante topológico. Por ejemplo cuando exigimos 
la conservación de la simetría de reversión temporal, el invariante topológico que 
describe el sistema debe ser de tipo Zo. En la Ref. [63] los autores definen un 
invariante de este tipo, que se relaciona con el índice de Kane-Mele, al igual que el 
número de winding se relaciona con el número de Chern. 

Los últimos resultados en el área muestran que además de los invariantes 
definidos anteriormente, puede que debamos sumar contribuciones en el caso que 
existan singularidades en el espectro del operador evolución a tiempos menores 
que el período del potencial [62]. La caracterización completa de la topología de 
los aislantes topológicos de Floquet aun no está concluida y sigue siendo un tema 
activo de investigación. 


6.2.2. El problema de la ocupación de los estados de Floquet 


La ocupación de los estados de no equilibrio en un sistema sujeto a un potencial 
periódico en el tiempo, es una cuestión altamente desafiante. De hecho, para poder 
hablar de un aislante en el sentido usual deberíamos poder establecer primero las 
ocupaciones de manera de tener ocupados todos los estados de un conjunto de 
bandas dejando un gap para excitaciones. En el desarrollo del campo de los llamados 
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aislantes topológicos de Floquet, la analogía llega hasta la definición de estados 
de una partícula similares a los encontrados en aislantes topológicos. El problema 
surge al intentar ocuparlos ya que la determinación de estas ocupaciones depende 
en general del setup experimental y los mecanismos de disipación presentes en 
el material. Este es un tema debatido actualmente con diferentes propuestas que 
mencionamos brevemente a continuación. 

Para empezar, podemos imaginarnos la evolución temporal del sistema, partien- 
do de un estado de equilibrio y luego encendiendo el potencial periódico hasta llegar 
al régimen estacionario. La conexión entre los estados de equilibrio iniciales y los 
estados estacionarios de Floquet no es directa y se ha demostrado que algunos 
estados de Floquet pueden ser no accesibles dependiendo de las ocupaciones del 
estado inicial y del protocolo de encendido del potencial [67, 75, 90] —aunque en 
presencia de interacciones el estado estacionario a tiempos largos es independiente 
del estado inicial [73]. Por otro lado, se argumenta que ante evolución unitaria no 
podemos cambiar la topología del estado bulk inicial, aunque en presencia de bordes 
el invariante topológico correcto (el índice de Bott) sí puede cambiar [68]. 

En los sistemas forzados de Floquet el potencial periódico «inyecta» energía 
continuamente al sistema. Si el sistema está aislado la ocupación de los estados 
tiende a la máxima entropía y el sistema se calienta hasta infinito [68]. Esta idealiza- 
ción es poco alentadora para la realización de los aislantes topológicos de Floquet 
en materia condensada, aunque al incluir otros mecanismos como interacciones de 
muchos cuerpos podemos probar que el sistema se estabiliza [69]. 

Deben encontrarse nuevos argumentos termodinámicos para definir las ocupa- 
ciones de los sistemas de Floquet, ya que no podemos aplicar directamente los 
argumentos para obtener la función de ocupación [67], principalmente porque el 
espectro de Floquet es no acotado. Al conectar el sistema a reservorios térmicos [70] 
e incluir los mecanismos de disipación(interacciones electrón-fonón, recombinación 
radiativa, emisión de fotones, scattering asistido por fotones) que entren en juego 
en el setup del sistema, podemos definir la ocupación de los estados de Floquet. 
Sin embargo, en general se observan excitaciones en el sistema aún cuando está 
conectado a reservorios bosónicos o fermiónicos de temperatura cero. Algunas 
aproximaciones proponen hacer una ingeniería del baño térmico para garantizar 
un estado estacionario donde sólo una banda de Floquet esté ocupada [71, 72], 
sintonizando apropiadamente los mecanismos de disipación, las propiedades de 
los reservorios térmicos y la configuración de acople entre el sistema forzado y el 
reservorio. 


6.2.3. Realización experimental 


En fotónica ya se ha reportado la realización de los aislantes topológicos de 
Floquet [40, 41]. En átomos fríos se ha llegado a simular el modelo de Haldane [91]. 
El Hamiltoniano del modelo de Haldane es análogo al Hamiltoniano efectivo para 
grafeno a frecuencias mucho más grandes que el ancho de banda [30]. Los hoppings 
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a segundos vecinos del modelo de Haldane son análogos, en el contexto de grafeno 


irradiado, a absorciones y reemisión de fotones. 


En materia condensada el primer experimento es el de Wang y col. [39], donde 


se abren los gaps en el espectro de Floquet-Bloch en la superficie de un aislante 


topológico tridimensional. A través de la técnica de "pump and probe”, se observa 


la relación de dispersión mediante ARPES para diferentes polarizaciones. Las 


excitaciones electrónicas son de tipo Dirac al igual que en grafeno y se observan 


gaps que dependen de la polarización en los cruces evitados, mientras que la 


radiación de polarización circular induce además un gap en el punto de Dirac. 


6.3. 
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A.1. Modelo tight-binding para grafeno bulk irradiado 


Para describir el grafeno irradiado utilizamos un modelo de red, con orbitales 7 y 
acoples a primeros vecinos que forman una red hexagonal. 


El Hamiltoniano de red 


Hi = Yes ae = paa le (A.1) 
1 (5,5) 

donde la sumatoria corre sobre los primeros vecinos, t; ¿ es la amplitud de hopping 

entre los sitios ubicados en r; y r;, y dado que son primeros vecinos, cada sitio 

pertenece a una subred distinta. Los sitios de tipo B están desplazados de los sitios 

de tipo A en las direcciones 


d¡ = ay(1,0) 
-1 3 
d) = al» E 
-1 v3 
ER 


Si los hoppings en cada dirección son homogéneos, t;, j = 1,2,3, podemos utilizar 
el teorema de Bloch y reducir el Hamiltoniano de la Ec. (A.1) al Hamiltoniano de 
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Bloch 


(A.2) 


set) EA e 


S(k)* EB 


donde S(k) = D-1 tje"%* y S(k)* su complejo conjugado. 
Podemos describir la radiación electromagnética mediante el potencial vector 


A(T). El efecto del potencial electromagnético es modular la fase de los hoppings 
t; a través de la sustitución de Peierls 


2 
t¿ —ti(T) =texp (an) : d;) 4 
Do 


donde Y, es el cuanto de flujo magnético, donde hemos evaluado la integral ya que 
el potencial vector es homogéneo en el plano (x, y). De esta forma obtenemos un 
Hamiltoniano de Bloch dependiente periódicamente del tiempo. 


E 0 S(k,T) 
A EA as) 


: 2 
S(k,T) = Y te dir x exp (Ac) : d;) 
1 0 


Para obtener la componente Hp” del Hamiltoniano de Floquet-Bloch, debemos 
calcular la componente (4H”""") de Fourier del Hamiltoniano de Bloch dependiente 
del tiempo. Esto implica obtener la componente (m — n) del hopping t;(7). Esta 
integral depende explícitamente del potencial vector 4(7), que para radiación cohe- 
rente puede escribirse como A(7) = 24, sin(Q7) + YAysin(Q7 + Ho), donde dy 
es una fase relativa entre las componentes, con y = O para polarización lineal y 
ó = +11/2 para polarización circular. La componente de Fourier de t;(7) es 


ti(r) ES E 
n 
e E it x €Xp (Ac d;)ar 
_ Al a Sy cd sin(Qr) ¿iz sin(Q7+$0) q7 (A.4) 
T 
E > fra ¿e (rd cos($0)) sin(O7) e sin($o) cos(07) q, 


jj 27 A . Jj _ 2m A ] : 
donde z; = 5 Al : dj Y 2, = 5, Ay8 : dj. Para resolver esta integral hacemos uso 
de las expansiones de Anger-Jacobi 


ei? cos(9) _ y ¡n Jn(2)er? 


NO 


ei sin(0) - y Jaen, 
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Insertando esta expansión obtenemos 
¿” - 1 a inQT E y j ( y) imQr e 
=p], te : [2 Jm(zs z; cos(p0))e 
[Erro jsp) ar 
E E i(n-(m+m/!)) QT ¡m/ Ta J J 
= el te 1" Jr(<sin(90)) Jm(zz, + 23 cos(Po)))d7 
=1 Y) mem” Jr (2ysin(do)) Jm(22, + 24 cos(d0))) 


m,m! 


=t Za la sin(Hp)) Ta + 2 cos(P0))) 


Con esta última expresión ya podemos calcular los coeficientes de la matriz de 
Floquet, aunque los cálculos son más simples si utilizamos la identidad [48] 


Pelo Y Tama) Imajeaer, 





donde z = y/a?+ 82-208 cos(6), y y = tan (Psin(Q)/(a — Bcos(2))). En esta 
caso es evidente que los acoples entre sitios que pertenecen a distintas réplicas 
es proporcional a la función de Bessel del orden de la diferencia entre las réplicas 


"oe Janmi). 


A.2. Modelo k . p para GNR's irradiadas 


En esta Sección deducimos cómo son los acoples inducidos por el campo elec- 
tromagnético entre los distintos modos transversales de una nanocinta de grafeno 
armchair y zigzag. Para obtener los modos transversales de las nanocintas seguimos 
las deducciones de la Ref. [51]. 


El Hamiltoniano k - p que describe el sistema a bajas energías cerca de los valles 
K, es 


Hp = hup (0% + 10,8): k, 


donde utilizamos la notación k para distinguir el operador momento de su autovalor, 
v=+*+l es el índice de valle para K,, y las matrices de Pauli 


E) 
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Para ambos valles en conjunto el Hamiltoniano k - p toma la forma 


0% 0.0 
kb. 0 0 0 
Hp = hu ca 
ME Duo de 
0 0 kk 


done £, = kz +14. 

Bajo la incidencia perpendicular de radiación, el campo electromagnético es 
homogéneo en el plano del grafeno y para describirlo utilizamos el potencial vector 
A(t) = Re(Agel*), con £2 la frecuencia de la radiación, y el vector Ay = 4,4 + AyÚ. 
La magnitud de las componentes A,, y A, define la dirección de la polarización. 
Si estas componentes son reales la polarización es lineal en la dirección de Ao, 
mientras que si existe una fase relativa entre ellas la radiación es elíptica o circular 
cuando la fase relativa es +7/2. 


Para describir el Hamiltoniano dependiente del tiempo utilizamos la sustitución 
canónica k > k + ¿A(t) y obtenemos 


HUrp(t) = hop (0,4 + 10,9): [k + , A), 


Con este la dependencia temporal explícita podemos armar el Hamiltoniano de 
Floquet truncado a las réplicas O y 1 


HO y! 
Hr=l 1 pro > 
VA HPA 


donde 


e LR 
AP. 0 00 

0. 0 0 47m 

0 AE 


V” =eup 


E) 


donde 47” es la componente de Fourier de orden m de 4,(t) = (A, +i4,) 0%. 


Los modos que pertenecen a cada réplica están definidos por 4H“, y se aco- 
plan entre ellos a través de V*!. Deducimos a continuación cómo son los modos 
transversales para una cinta armchair. 


A.2.1. Nanocintas de grafeno armchair 


Las funciones de onda expandidas a bajas energías pueden ser descritas por 


v(r)= 2 VA(RaJjo(r - Ra) + Z Vp(Ra)óo(r - Ra), 
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donde 


dalr) =rE(r) ir FR (e) 


pa(r) ¡TE (r) + AR (r), 


son las funciones envolventes que describen la amplitud y fase de la función de 
onda sobre los sitos A y B. 


Definimos los puntos XK, con puntos equivalentes a la definición de los capítulos 
2 y 3. Para simplificar los cálculos elegimos los puntos XK, sobre el eje k,, siendo 
K_=-Kjyy K, = Ky, con K = Se. El vector que representa la función de onda es 


(60) EE) EE EE), 


expandidas alrededor de los puntos K_ y K, de la primer zona de Brillouin para 
los sitios de tipo A y B. 


Vamos a tratar el caso de nanocintas, donde una de las direcciones está con- 
finada, por lo tanto queremos las funciones envolventes definidas en función de 
la posición, y deben satisfacer la ecuación diferencial definida por el Hamiltoniano, 
junto a las condiciones de borde apropiadas. Dado que cambiamos entre la base de 
momentos y de posiciones es útil recordar que k es el operador momento tal que 
b==iyi 


Dado que el sistema tiene invarianza traslacional a lo largo del eje x pode- 
mos separar la componente longitudinal de la componente transversal debido al 
confinamiento en y, para cada valle 


Fr) = e of (y) 
FE-(r) = "9 f-(y) 
Flor) = rot (y) 


Fo (1) = "04 (y), 


(A.5) 


donde mantenemos k.. y k”, distinguiendo cada valle por el momento, y las funciones 
ri (0) definen los modos transversales. 


Condiciones de borde de una nanocinta armchair 


Utilizamos las condiciones de borde de masa infinita, de manera que la función de 
onda se anula en sitios virtuales adyacentes a los bordes de la muestra, localizados 
- - N+1 

en los puntos y = 0, W, con W =a 





. Luego exigimos que las funciones de onda 


97 


A. Modelos para describir el grafeno irradiado A.2 Modelo k - p para GNR's irradiadas 


98 


satisfagan 
valr,y=0)=0 
,y=0)=0=> 
p(z, y = 0) e 
o 
Vplx,y = W)=0, 


o equivalentemente 


0=D%-(0) -i0%*(0) 
0=iD0%-(0)+0%*(0) 
O = eEWgK- (17) -¡e EWgK+ (py) 


0=i Wok (17) + e kWok+ (7), 


Como podemos ver de la Ec. (A.6), las condiciones de borde armchair no acoplan 
las funciones envolventes de sitios de distinto tipo, pero sí acopla las funciones 
envolventes dentro de la misma familia de sitios A y B que pertenecen a diferentes 
valles K_ y K.,. 

Las condiciones de borde implican que la componente transversal de las solucio- 
nes para Phil y Ph están indexadas por un k,, y son discretas, mientras que 
la componente paralela al borde debe ser idéntica para los valles K_ y K,, por lo 
que k.. = k/, en la Ec.(A.5). 


Las soluciones toman la forma 


pa seg 7ikny 
K- ik 

F a e ny 

Fy* ise: e ny 
K; _s ¿1kny 

Fr je tn 


donde s = +1 es el índice de banda, C! = —— una constante de normalización, 


VaW 


¿sa Ote + ik) 


y k2 + k2 


a N+1 
nT W + 


a =tan Kn /Ka 





siendo N el número de filas de sitios a lo largo del ancho de la cinta. 

Dado que la energía es e = y/k2 + k2, podemos ver que el sistema será metálico 
si existe uno o más modos que se aproximen continuamente a e = 0. Esto sucede 
solo cuando existen una solución con k,, = O De la expresión anterior encontramos 
que 

nT dr N+1 3n-2(N +1) 


a E hi = 2 
W 3a a A "SN AD 
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Entonces, las cintas aGNR son metálicas sólo cuando 2(N + 1) = 3M, para algún 
M entero, que equivale a N- mód (3) = 2. 


Transiciones inducidas por la radiación 


Los modos transversales son ortogonales entre ellos, pero el potencial vector 
puede mezclarlos. Para ver esto analizamos el vector resultante de aplicar la matriz 
V”" a un modo |) = |k.,,n, s) 


0 AP 0 0 sel 2 g7ikny 
AP. 0) 0 0 e7ikny 
V"lw)=C0| + A 
Iv) 0. Dd az mi sele ei ny 
0 0 47 0 ie ny 
AMgikny 
-C AM geiQ giiny 


Am (cies) |? 
Am (i sel eliny) 
donde k., es el vector de onda paralelo al borde, n el índice del modo transversal y s 


representa la banda de valencia s = —-1 o de conducción s = 1. Vemos el caso m = 1 
y descartamos el índice para simplificar la notación. 


El elemento de matriz entre dos modos |) = [k.,, n, s) pertenecientes a la réplica 
0, y |”) = |k., n', s”) perteneciente a la réplica 1 definen la amplitud de acople entre 
los modos. Para que se crucen estos modos en e = h(2/2 debemos pedir s = 1, 
s' = —1, pero mantenemos los índices para obtener una expresión general. 


(v'1V |) W ia! ik, -ikn ik ia ik 
POr = 0 [se emYp ¿Uy end A se “e sl dy 
W do 7 y Lia 
ls me [agas en py ¿ny + eng p_gelergitns | dy 
y pWw 
= gia y [Aye (Ekn—k 1 )y + A od dy 
W 
+se" J [4,07 o=hnds + Allin] dy 
W o 
= y [2er A, cos((kn — kn1)y) — Ay sin((k - le )y) 
+ 2560 [Az cos((K, - kn)Y) + Ay sin((ky - h)4)]) dy 


Ahora podemos distinguir dos casos: el acople entre los mismos modos n = n' o 
entre distintos modos n + n'. Además, el acople entre modos distintos decae como 
1/(n — n”), por lo que los modos más cercanos a n serán los que se vean afectados 
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por un acople mayor. Los elementos de matriz son 


2 ia, ¿ia y 9sin((n-n')m) 
MA de eupO (4, (se +s'e )2 (na W 


(Vw) = Ay) 


a 
sin=n" > evpC? A, (set+ se LO ) 2W 
En una forma más compacta vemos que 


1- (y 


(n- nr 


Ay 


. j =A r 
ifn+n! => cup (sei2 seria ) 


(wV |) = 


. As : Se 
ifn=n! => cup (seio + seta ) 


donde podemos ver que sólo la componente 4. genera transiciones intra-modo, 
es decir entre los mismos modos n = n”, y que sólo la componente A, genera 
transiciones inter-modo, entre modos distintos n + n”. 

Si observamos el cruce de bandas en e = A(2/2, entonces s = 1 y s' = —1, y 
podemos ver que para las transiciones intra-modo n = n' = a = o, por lo tanto 


Ag ia -iQ 
(IV |) cup (40-219) 
Ag 
Un 2i sin(a) 


A 
eup— 2i En 


2 E 


por lo que el acople depende de la componente transversal y se anula cuando k:,, = 0, 


es decir que el modo metálico no se acopla a la radiación de polarización lineal en 
la dirección paralela al borde armchair. 


A.2.2. Nanocintas de grafeno zigzag 


En el caso de las cintas zigzag el borde es perpendicular a la dirección de las 
cintas armchair. Para mantener la coordenada «x en la dirección de la cinta rotamos 
los ejes tal que ahora K, = +K'% Entonces, la función de onda es 


p(r) = 2 VA(Ra)olr - Ra) + 2 Vp(Rp)o(r - Rp) 


donde 


Ppa(r)= KT PE(r) +" EÉ(r) 
Uple) =- KT FE (r) + KE (r) 


Las condiciones de borde de masa infinita exigen que la función de onda se 
anule en sitios virtuales adyacentes al borde de la cinta, de modo que una cinta 
zigzag terminada en sitios de tipo 4 debe anularse en los sitios adyacentes, todos 
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pertenecientes a la familia de sitios B. Los sitios virtuales se conforman por una fila 
+2 





ads Lee z 3N 
de sitios B en y = 0, y una fila de sitios A y = W = a . Luego, las funciones 


envolventes deben satisfacer 


E 9%(0)=0 
ente y=0)=0>| 9%-(0) =0 

aa Ja (0)=0 
va J50 Ed 


Estas condiciones de borde no acopla las funciones envolventes expandidas alrede- 
dor de los valles K,, aunque las funciones envolventes para los sitios A y B están 
acopladas por el Hamiltoniano. 


La solución general para cada valle K, toma la forma 


pÍ=(y) =+sC sinh(z (y — W)) 
pÍ=(y) = Csinh(z y), 


donde s es el índice de banda, y z es solución de la ecuación trascendental z = 
kg tanh(2W), indexadas por z,,. Esta componente transversal puede tener valores 
reales como imaginarios. Cuando toma un valor real representa un modo localizado 
en los sitios B cerca del borde y = 0, y un modo localizado en sitios 4 cerca del 
borde y = W. En esta sección nos centraremos en las soluciones con z,, = ik, que 
son los modos extendidos que se acoplan en e = A0/2. 


Los modos extendidos cumplen con 


04 (y) = FOssin(kn (y- W)) 
DÍ (y) = Csin(k, y). 


Transiciones inducidas por la radiación 


La amplitud de acople entre distintos modos del valle K_ está dado por el 
elemento de matriz entre |) = |k..,n, s) pertenecientes a la réplica 0, y |v”) = 
Ik., n/, s”) perteneciente a la réplica 1. Siguiendo los pasos que tomamos para las 
cintas armchair podemos ver que 

Gala 
(IV Iw) = e 


n n 


py 5 in(2k, W DA ze 
con constantes de normalización C,, = 1/1/W - A) Las transiciones permiti- 


das en las cintas zigzag son menos obvias de analizar, pero podemos ver que las 


(A, s+iA_ s')) x (er sin (k, W) — k,, sin (ki, W)) 


transiciones decaen con 1/(k? — k?,), lo que significa que los modos más cercanos 
serán los de mayor acople para cualquier polarización, es decir, sin asumir 4, = 00 
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Ay = 0. Tomando el límite para transiciones entre los mismos modos n = n” 


(UV 0) = CuCve ni y) x (1 cos(tntr) - DE) 
n 
donde vemos que la amplitud de acople para las transiciones intra-modo sólo 
depende de la componente A, de la radiación. 
La polarización perpendicular al borde produce las transiciones intra-modo para 
las cintas zigzag, mientras que el comportamiento opuesto ocurre para las cintas 


armchair, donde sólo la polarización paralela al borde produce estas transiciones. 


